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Abstract 

One of the most efficient methods for the evaluation of the vacuum expectation values for 
physical observables in the Casimir effect is based on using the Abel-Plana summation for- 
mula. This enables to derive the renormalized quantities in a manifestly cutoff independent 
way and to present them in the form of strongly convergent integrals. However, applications 
of the Abel-Plana formula, in its usual form, are restricted by simple geometries when the 
eigenmodes have a simple dependence on quantum numbers. The author generalized the 
Abel-Plana formula which essentially enlarges its application range. Based on this general- 
ization, formulae have been obtained for various types of scries over the zeros of combinations 
of Bessel functions and for integrals involving these functions. It has been shown that these 
results generalize the special cases existing in literature. Further, the derived summation 
formulae have been used to summarize series arising in the direct mode summation approach 
to the Casimir effect for spherically and cylindrically symmetric boundaries, for boundaries 
moving with uniform proper acceleration, and in various braneworld scenarios. This allows 
to extract from the vacuum expectation values of local physical observables the parts cor- 
responding to the geometry without boundaries and to present the boundary-induced parts 
in terms of integrals strongly convergent for the points away from the boundaries. As a 
result, the renormalization procedure for these observables is reduced to the corresponding 
procedure for bulks without boundaries. The present paper reviews these results. We also 
aim to collect the results on vacuum expectation values for local physical observables such as 
the field square and the energy-momentum tensor in manifolds with boundaries for various 
bulk and boundary geometries. 
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1 Introduction 



In many physical problems we need to consider the model on background of manifolds with 
boundaries on which the dynamical variables satisfy some prescribed boundary conditions. In 
quantum field theory the imposition of boundary conditions leads to the modification of the 
spectrum for the zero-point fluctuations and results in the shift in the vacuum expectation 
values for physical observables. In particular, vacuum forces arise acting on the constraining 
boundary. This is the familiar Casimir effect [T]. The Casimir effect is among the most interest- 
ing consequences of quantum field theory and is essentially the only macroscopic manifestation 
of the nontrivial properties of the physical vacuum. Since the original work by Casimir many 
theoretical and experimental works have been done on this problem, including various bulk and 
boundary geometries and boundary conditions (see, for instance, O [3l HJ EJ [H [7j and refer- 
ences therein). Many different approaches have been used: direct mode summation method. 
Green function formalism, multiple scattering expansions, heat-kernel series, zeta function reg- 
ularization technique, etc. Advanced field-theoretical methods have been developed for Casimir 
calculations during the past years [H [9l [10] . From a general theoretical point of view the main 
point here is the unique separation and subsequent removing of the divergences. Within the 
framework of the mode summation method in calculations of the expectation values for physi- 
cal observables, such as energy-momentum tensor, one often needs to sum over the values of a 
certain function at integer points, and then subtract the corresponding quantity for unbounded 
space (usually presented in terms of integrals). Practically, the sum and integral, taken sepa- 
rately, diverge and some physically motivated procedure to handle the finite result, is needed 
(for a discussion of different methods to evaluate this finite difference see [H]). For a number of 
geometries one of the most convenient methods to obtain such renormalized values of the mode 
sums is based on the use of the Abel-Plana summation formula (AFF) [12l [131 [IH [13 [13 [HI- 
The development of this formula dates back to Plana in 1820 [18] and further has been recon- 
sidered by Abel [I9j . Cauchy [20j and Kronecker [21]. The history of the APF is discussed in 
detail in [22] where interesting applications are given as well. Some formulae for the gamma and 
zeta functions are obtained as direct consequences of the APF (see, for instance, [13 [23]). This 
summation formula significantly and consistently improves the convergence and accuracy for the 
slowly convergent series (see, for instance, t24j). The summation method based on the APF is 
so efficient that was considered for the development of software for special functions [25]. In 
|26j the APF has been used for the renormalization of the scalar field energy-momentum tensor 
on backgrounds of various Friedmann cosmological models. Further applications of the APF in 
physical problems related to the Casimir effect for flat boundary geometries and topologically 
non-trivial spaces with corresponding references can be found in [21 [3]. The application of the 
APF in these problems allows (i) to extract in a cutoff independent way the Minkowski vac- 
uum part and (ii) to obtain for the renormalized part strongly convergent integrals useful, in 
particular, for numerical calculations. 

The applications of the APF in its usual form is restricted to the problems where the eigen- 
modes have simple dependence on quantum numbers and the normal modes are explicitly known. 
For the case of curved boundaries and for mixed boundary conditions the normal modes are given 
implicitly as zeroes of the corresponding eigenfunctions or their combinations. The necessity 
for a generalization of the Abel-Plana summation procedure arises already in the case of plane 
boundaries with Robin or non-local boundary conditions, where the eigenmodes are given im- 
plicitly as solutions of transcendental equation. In problems with spherically and cylindrically 
symmetric boundaries the eigenmodes are the zeroes of the cylinder functions and the expec- 
tation values of physical observables contain series over these zeroes. To include more general 
class of problems, in [27] the APF has been generalized (see also |28j). The generalized version 
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contains two meromorphic functions and the APF is obtained by specifying one of them (for 
other generahzations of the APF see O \TT\ 129]). Applying the generahzed formula to Bessel 
functions, in [271 EB] summation formulae are obtained for the series over the zeros of various 
combinations of these functions (for a review with physical applications see also Ref. |30l [3T] ) . 
In particular, formulae for the summation of the Fourier-Bessel and Dini series are derived. It 
has been shown that from the generalized formula interesting results can be obtained for infi- 
nite integrals involving Bessel functions. The summation formulae derived from the generalized 
Abel-Plana formula have been applied for the evaluation of the vacuum expectation values for 
local physical observables in the Casimir effect for plane boundaries with Robin or non-local 
boundary conditions [32l |33], for spherically [3H [35l [361 133 EH] and cylindrically symmetric 
|391 l40l [411 142] boundaries on the Minkowski bulk for both scalar and electromagnetic fields. 
As in the case of the APF, the use of the generalized formula allows to extract in a manifestly 
cutoff independent way the contribution of the unbounded space and to present the boundary- 
induced parts in terms of exponentially converging integrals. The case of cylindrical boundaries 
in topologically nontrivial background of the cosmic string is considered in [431 E] • By making 
use of the generalized Abel-Plana formula, the vacuum expectation values of the field square 
and the energy-momentum tensor in closely related but more complicated geometry of a wedge 
with cylindrical boundary are investigated in [451 l46l I47j for both scalar and electromagnetic 
fields. In [481 B9l [50] summation formulae for the series over the zeroes of the modified Bessel 
functions with an imaginary order are derived by using the generalized Abel-Plana formula. This 
type of series arise in the evaluation of the vacuum expectation values induced by plane bound- 
aries uniformly accelerated through the Fulling-Rindler vacuum. In all considered examples the 
background spacetime is flat. 

For curved backgrounds with boundaries the application of the generalized Abel-Plana for- 
mula extracts from the expectation values the parts which correspond to the polarization of the 
vacuum in the situation without boundaries and the boundary induced parts are presented in 
terms of integrals which are convergent for points away from the boundaries. As a result, the 
renormalization is necessary for the boundary-free parts only and this procedure is the same 
as that in quantum field theory without boundaries. Examples with both bulk and boundary 
contributions to the vacuum polarizations are considered in [511 [52l |53l [SH [55] . In these papers 
the background spacetime is generated by a global monopole and the one-loop quantum effects 
for both scalar and fermionic fields are investigated induced by spherical boundaries concentric 
with the global monopole. Another class of problems where the application of the generalized 
Abel-Plana formula provides an efficient way for the evaluation of the vacuum expectation values 
is considered in [561 EZ] • In these papers braneworld models with two parallel branes on anti-de 
Sitter bulk are discussed. The corresponding mode-sums for physical observables bilinear in 
the field contain series over the zeroes of cylinder functions which are summarized by using the 
generalized Abel-Plana formula. The geometry of spherical branes in Rindler-like spacetimes 
is considered in [58]. In [59] from the generalized Abel-Plana formula a summation formula is 
derived over the eigenmodes of a dielectric cylinder and this formula is applied for the evalua- 
tion of the radiation intensity from a point charge orbiting along a helical trajectory inside the 
cylinder. 

The present paper reviews these results and is organized as follows. In Section [2] the gen- 
eralized Abel-Plana formula is derived and it is shown that the APF is obtained by a special 
choice of one of the meromorphic functions entering in the generalized formula. It is indicated 
how to generalize the APF for functions having poles. The applications of the generalized for- 
mula to the Bessel functions are considered in Section [3l We derive two formulae for the sums 
over zeros of the function AJ^{z) + BzJl{z), where Ju{z) is the Bessel function. Specific ex- 
amples of applications of the general formulae are considered. In Section [H by special choice 
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of the function g{z) summation formulae are derived for the series over zeros of the function 
Jiy{z)Yiy{Xz) — J^{\z)Y^(z) and similar combinations with Bessel functions derivatives. Special 
examples are considered. In Section [S] we consider the application of the generalized Abel-Plana 
formula to the series over the zeros of the modified Bessel functions with an imaginary order. 
The applications to the integrals involving Bessel functions and their combinations are discussed 
in Section [6l A number of interesting results for these integrals are presented. Specific exam- 
ples of applying these general formulae are described in Section [71 In Section [HI by using the 
generalized Abel-Plana formula two theorems are proved for the integrals involving the function 
Ji^(z)Yfj_{Xz) — J^{Xz)Yu{z) and their applications are considered. 

In the second part of the paper we consider the applications of the generalized Abel-Plana 
formula to physical problems, mainly for the investigation of vacuum expectation values of phys- 
ical observables on manifolds with boundaries of different geometries and boundary conditions. 
First, in Section [9] we describe the general procedure for the evaluation of the vacuum expecta- 
tion values of physical observables bilinear in the field operator. Applications of the APF for the 
evaluation of the the vacuum expectation values of the field square and the energy-momentum 
tensor for a scalar field in the spacetime with topology x and for the geometry of two 
parallel Dirichlet and Neumann plates are given in Section [TOl In Section [T2l the summation 
formulae over the zeroes of combinations of the cylinder functions are applied for the investiga- 
tion of the scalar vacuum densities induced by spherical boundaries on background of the global 
monopole spacetime. The corresponding problem for the fermionic field with bag boundary 
conditions is considered in Section [131 In the latter case the eigenmodes are the zeroes of the 
more complicated combination of the Bessel function and its derivative. Summation formulae 
for series over these zeroes are derived from the generalized Abel-Plana formula. The application 
to the electromagnetic Casimir effect for spherical boundaries in background of the Minkowski 
spacetime is given in Section [TH Further we consider problems with cylindrical boundaries. In 
Section [15l we study the vacuum polarization effects for a scalar field induced by a cylindrical 
shell in the cosmic string spacetime assuming Robin boundary condition on the shell. The corre- 
sponding problem for the electromagnetic field and perfectly conducting cylindrical boundary is 
considered in Section [TBI Section [iTl is devoted to the investigation of the vacuum densities in the 
region between two coaxial cylindrical surfaces in the Minkowski spacetime for both scalar and 
electromagnetic fields. In Sections [181 and [191 we consider the Wightman function, the vacuum 
expectation values of the field square and the energy-momentum tensor for scalar and electro- 
magnetic fields induced by plane boundaries uniformly accelerated through the Fulling-Rindler 
vacuum. The application of the generalized Abel-Plana formula for the evaluation of the vacuum 
densities in braneworld scenarios is described in Section [20l for parallel branes on background 
of anti-de Sitter spacetime and in Section [2T1 for spherical branes in Rindler-like spacetimes. In 
Section [22l a summation formula for the series over the eigenmodes of a dielectric cylinder is 
derived and this formula is applied for the investigation of the radiation intensity from a charged 
particle moving along a helical orbit inside the cylinder. Section [23l concludes the main results 
of the paper. 

2 Generalized Abel-Plana formula 

Let f{z) and g{z) be meromorphic functions for a ^ x ^ 6 in the complex plane z = x + iy. Let 
us denote by Zf^k and Zg^k the poles of /(z) and g{z) in the strip a < x < b. We assume that 
Im Zf^k (see, however, the Remark to Lemma). 
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Lemma. If functions f (z) and g{z) satisfy condition 



h—foo 

then the following formula takes place 



b±ih 

lim / [g{z) ± f{z)] dz = 0, (2.1) 

azLih 



/ fix)dx = R[fiz),giz)] - - / [g{u)+a{z)f{u)]:zlXldz, a{z) ^ sgn(Imz), (2.2) 



where 



R[f{z),g{z)]=^i 



y~] Res g{z) + (j{zf^k) Res /( 



(2.3) 



Proof. Consider a rectangle Ch with vertices azizih, 6 it ih, described in the positive sense. In 
accordance to the residue theorem 

/ dzg{z) =2Triy^ Res g{z), (2.4) 

where the rhs contains the sum over poles within Ch- Let and denote the upper and 
lower halves of this contour. Then one has 

/ dzgiz)= V 1/ dz[g{z) + af{z)]-a f dzf{z)\. (2.5) 

By the same residue theorem 

/ dzfiz)- [ dzfiz) = 2 [' fix)dx-27ny^a{zf^k)Res fiz). (2.6) 
Jc- Jet Ja Y 

Then 

/ dz [g{z) ± f{z)] =± dz [g{u) ± f{u)Zz':Xl ^ / dz [giz] ± /(z)]. (2.7) 

Jet Jo Ja±ih 

Combining these results and allowing — > c« in (12.41) one obtains formula (12.21) . ■ 

If the functions f{z) and g{z) have poles with Rezj^k = a,b {j = f,g) the contour has to 
pass round these points on the right or left, correspondingly. 

Remark. Formula (j2.2p is also valid when the function f{z) has real poles •zj-''^, ^^■^^fn ~ 
in the region a < Re z < 6 if the main part of its Laurent expansion near these poles does not 
contain even powers of z — Zj^j^. In this case on the left of formula (|2.2|) the integral is meant in 
the sense of the principal value, which exists as a consequence of the above mentioned condition. 
For brevity let us consider the case of a single pole z = zq. One has 

dzfiz)- dzfiz) = 2 dzfiz) + 2 dzfiz) 



a J z()+p 



-2^iVa(z/,fc) Res /(z)+ /" dzfiz)+ [ dzfiz), (2.^ 
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with contours and being the upper and lower circular arcs (with center at z = zq) joining 
the points zq — p and zq + p. By taking into account that for odd negative I 



dz {z - zoY + dz{z- zoY = 0, (2.9) 
Jr- 

in the limit p ^ we obtain the required result. ■ 

In the following, on the left of ()2.2p we will write p.v. /(x), assuming that this integral 
converges in the sense of the principal value. As a direct consequence of Lemma one obtains 



Theorem 1. If in addition to the conditions of Lemma one has 

6— »oo 

then 



/•bHoo 

lim / dz [g{z) ± f{z)] = 0, (2.10) 



hm p.v. / dx fix) - R[f{z),g{z)] = - / dz [g{z) + a{z)f{z)], (2.11) 

where on the left R[f{z),g{z)] is defined as ^2.3\) . a < Re z/^fc, Re z^^fc < b, and the summation 
goes over poles Zf^^ c-i^d Zg^f^ arranged in order Rezj^^ ^ ReZj^fc+i, j = f,g. 

Proof. It is sufficient to insert in general formula (j2.2p b ^ oo and to use the condition (12.10p . 
The order of summation in R[f{z),g{z)] is determined by the choice of the integration contour 
Ch and by limiting transition 6 ^ oo. ■ 

We will refer formula ()2.1ip as generalized Abel-Plana formula (GAFF) as for b = n + a, 
< a < 1, 

g{z) = -if{z)cot7:z, (2.12) 

with an analytic function f{z), from (j2.1ip the Abel-Flana formula (AFF) |12 p i3 1 [H] is obtained 
in the form 



lim 

n— >oo 



^f{^)-l dxf{x) = Y dz f{z) {cot IT z - i) 

k=l ^ 

/ dzf{z){cotTTZ + i). (2.13) 

2« J a 

Another form of the AFF, given in |15l I16j . is obtained directly from formula ()2.2p taking the 
function g{z) in the form ()2.12p with an analytic function f{z) and a = m, b = n with integers m 
and n. In this case the points z = m and z = n are poles of the function g{z) and in the integral 
on the right of formula (j2.2p they should be avoided by small semicircles from the right and from 
the left respectively. The integrals over these semicircles give the contribution —f{m) — f{n) 
and one obtains the formula 



Y,fi.k) = / dxf{x) + -f{m) + -f[n) 

k=m Z / 



^. f{m + iz) - f{n + iz) - f{m - iz) + f{n - iz) 
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Expanding the functions in the numerator of the last integral over z and interchanging the order 
of the summation and integration, the integral is expressed in terms of the Bernoulli numbers 
B^y As a result, from (|2.14|) we obtain the formula 

" rn 1 1 °° R r 1 

Y^f(k)= \ dxf{x) + -f{m) + -f{n)-Y,7:^ f^'^^'\m) - f^'^-'\n) , (2.15) 

k=m -^"^ 11 .^^ (ZJ). 

which is a special case of the Euler-Maclaurin formula [16] (for a relationship between the APF 
and Euler-Maclaurin summation formula see also [60]). 

A useful form of (j2.13p may be obtained performing the limit a ^ 0. By taking into account 
that the point z = is a pole for the integrands and therefore has to be circled by arcs of the 
small circle Cp on the right and performing /? — > 0, one obtains 

f; fik) = r dx fix) + i/(0) + ^ Tdx ^^'''}~J}~'''\ (2.16) 
fc=o •'^ •'^ 

Note that, now condition (j2.ip is satisfied if 

lim e-2^l2^l|/(x + %)| = 0, (2.17) 

uniformly in any finite interval of x. Formula (j2.16p is the most frequently used form of the APF 
in physical applications. Another useful form (in particular for fermionic field calculations) to 
sum over the values of an analytic function at half of an odd integer points is obtained from 
(f2TT]) taking 

g{z) = if (z) tan TTZ, (2.18) 
with an analytic function f{z). This leads to the summation formula (see also [21 [3]) 

E + V2) = f dx fix) -i[dx ^^'"l~Jl~'''\ (2.19) 

By adding to the rhs of (|2.16p the term 

Res 

and to the rhs of (|2.19p the term 

aizf^k) Res 

Z=ZFl 

k 



vrV Res — -fiz), (2.20) 

^-^ z=Zf^^: smvrz 



mVa(z/,fc) Res -fiz), (2.21) 

^-^ z=zf^k cosvrz 



these formulae can be generalized to functions fiz) that have non-real poles zj^^, Rezj^^ > 0. 
The second generalization with an application to the problem of diffraction scattering of charged 
particles is given in [61]. Another generalization of the APF, given in [62], is obtained from (j2.1ip 
taking in this formula giz) = —ifiz) cot7r(z — /?), < /3 < 1, with an analytic function fiz). 

As a next consequence of ()2.1ip . a summation formula can be obtained over the points 
Zk, Rezfc > at which the analytic function s(z) takes integer values, s(zfc) is an integer, and 
s'izk) 7^ 0. Taking in (j2.1ip giz) = —ifiz) cot7rs(z), one obtains the following formula [63] 



> —r, — 7=w + / f[x)dx+ / dx 
^s'izk) io Jo 



fiix) fi-ix) 



-2-Kis{ix) ]^ g27ris(— ix) ]^ 



(2.22) 
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where 

r 0, if .(0)/0,±l,±2,... 

1 /(0)/[2s'(0)], if K0) = 0,±l,±2,... • 

For s{z) = z we return to the APF in the usual form. An example for the application of this 
formula to the Casimir effect is given in ^63j. 



3 Applications to Bessel functions 



Formula (|2.1ip contains two meromorphic functions and is too general. In order to obtain more 
special consequences we have to specify one of them. As we have seen in the previous section, one 
of possible ways leads to the APF. Here we will consider another choices of the function g{z) and 
will obtain useful formulae for the series over zeros of Bessel functions and their combinations, 
as well as formulae for integrals involving these functions. 

First of all, to simplify the formulae let us introduce the notation 



F{z) = AF{z) + BzF'{z) 



(3.1) 



for a given function F{z)^ where the prime denotes the derivative with respect to the argument 
of the function, A and B are real constants. As a function g{z) in the GAPF let us choose 



mm. 

Ju{z) 



(3.2) 



where Ju{z) and Y^{z) are Bessel functions of the first and second (Neumann function) kind. 
For the sum and difference on the right of ()2.1ip one obtains 



fiz) 



-1)'5(^) 



Hi'\z) 

Ju{z) 



fiz), k = l,2, 



(3.3) 



with hI^^ (z) and (z) being Bessel functions of the third kind or Hankel functions. For such 
a choice the integrals (|2.ip and (j2.10p can be estimated by using the asymptotic formulae for 
Bessel functions for fixed v and \z\ — > oo (see, for example, |6ll [65]). It is easy to see that 
conditions ()2.ip and (j2.10p are satisfied if the function f{z) is restricted to one of the following 
constraints 



|/(z)| <e(x)e^l^ 

where c < 2, a > 1 and £(2;) 
functions it follows that 



or 



00, 



(3.4) 



1/(^)1 < z = x + iy, 

for X — > 00. Indeed, from asymptotic expressions for Bessel 



b±ih 
bzbioo 



dz [g{z)±fiz)] 
dz [g{z)±f{z)] 





j 




J a 




Jo 



'•^ , 0\x±ih) 
dx 



dx 



Ju{x =b ih) 
(6 lb ix 



Ju{b ± ix) 



f{x±ih) 
^-f{b±ix] 



< 



< 



Nie{b) 



(3.5) 
(3.6) 



with constants Mi, M[, Ni, N[, and j = 1 {j = 2) corresponds to the upper (lower) sign. 

Let us denote by X^^k 7^ 0, A; = 1, 2, 3 ... , the zeros of Ju{z) in the right half-plane, arranged in 
ascending order of the real part. Re X^^k ^ Re X,^^k+i, (if some of these zeros lie on the imaginary 
axis we will take only zeros with positive imaginary part). All these zeros are simple. Note that 
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for real v > —1 the function Ju{z) has only real zeros, except for the case A/B + u < when 
there are two purely imaginary zeros [M]. By using the Wronskian VF[Ji,(z), yjy(z)] = 2/ttz, for 
(1231) one finds 

R[f{z),g{z)]=2Y,TAK,k)f{K,k)+riAfiz)], (3.7) 

k 

where we have introduced the notations 

^lmzk>0 Jv\z) ^lmzk<0 Jv\z) 

-vrj; Res /(z)^. (3.9) 

Here (7^ Xy^i) are the poles of the function f{z) in the region Rez > a > 0. Substituting ()3.7p 
into (j2.1ip . we obtain that for the function f[z) meromorphic in the half-plane Rez ^ a and 
satisfying condition (13. 4p . the following formula takes place 

dxf{x) ) = 



( n t 

hm <^ 2 V r,(A,,fc)/(A,,fc) + n,[f{z)] - p.v. / 



= '^^/^ 7 ^ N - 77 / dzf{z)~—, 3.10) 

2 7a -^i'(^) 2 J^(2:) 

where ReA,y,rri-i < a < ReA;y,m, ReA;^^„ < b < ReA,y,n+i, a < Rezfc < 6. We will apply this 
formula to the function f{z) meromorphic in the half-plane Rez ^ taking a ^ 0. Let us 
consider two cases separately. 

3.1 Case (a) 

Let f{z) has no poles on the imaginary axis, except possibly at z = 0, and 

/(ze^*) = -e2'^^7(z) + 0(2'^''), z^O (3.11) 

(this condition is trivially satisfied for the function f{z) = 0(2^")), with 

^r2|Rez^|-l, for integer I/, 

^ \ Heu + |Rez/| — 1, for non-integer u. 

Under this condition, for values v, for which Ju{z) has no purely imaginary zeros, the rhs of Eq. 
(j3.10p in the limit a — > can be presented in the form 



1 f°° K (t) r 

TT Jp I^[x) L 



+ dzf{z)^-^- dzfiz)^-^, (3.13) 

with 7^ and 7" being upper and lower halves of the semicircle in the right half-plane with radius 
p and with center at point z = 0, described in the positive sense with respect to this point. In 
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we have introduced modified Bessel functions luiz) and Ky{z) It follows from (j3.1ip 
that for z — > 

^TT^^^) = ''f!^"-r? /(^^"") + (3-14) 

Jv{.z) Jv{ze '^'j 

From here for p — > one finds 

Du = dzf{z) - dzf{z) = -7rRes/(z)^— ■ (3.15) 



Indeed, 



= / dzf{z)^^il^+ [ dz/(ze— 



7p Jv{z) J'Y+^ Jv{ze 

[ dzf{z)^i]^+[ dzo{z-') 



dzf{z)Sf\+[ dzo{z-'), (3.16) 

where 7^ (7ip! see below) is the upper (lower) half of the semicircle with radius p in the left 
half-plane with the center at 2; = (described in the positive sense). In the last equality we 
have used the condition that integral P-v.Jq dxf{x) converges at lower limit. In a similar way 
it can be seen that 

D, = i[ dzf{z)^M+j dzo{z-^). (3.17) 

•'ip+llp Jp\z} J^i^^ 

Combining the last two results we obtain (j3.15p in the limit p ^ 0. By using (|3.10p . (|3.13p and 
(|3T3|) we have [23: 

Theorem 2. If f(z) is a single valued analytic function in the half-plane Rez ^ (with possible 
branch point at z = 0) except the poles z^ ^u,i), Rezfc > (for the case of function f{z) 
having purely imaginary poles see Remark after Theorem 3), and satisfy conditions 1^3.4^ and 
^3.11\) . then in the case of u for which the function Ju{z) has no purely imaginary zeros, the 
following formula is valid 

lim \2S2T,{X,,k)f{K,k) + riAf{z)]-Y>.^.tf{x)dx\ = 

= - - r d-Wr U—f{xe^) + e^-^f{xe-H , (3.18) 

2z=o J^{z) vr Jo luix) L J 

where on the left ^eXy^n < b < ReAi^^„+i, < Rez^ < 6, and T^^Xp^k) md riy[f{z)\ are 
determined by relations 13. 8\) and 13. 9\) . 

Under the condition p. lip the integral on the right converges at lower limit. Recall that we 
assume the existence of the integral on the left as well (see section 2). Formula (I3.18P and 
analogous ones given below are especially useful for numerical calculations of the sums over X^^k^ 
as under the first conditions in ()3.4p the integral on the right converges exponentially fast at the 
upper limit. Below, in Section [13] we will see that summation formula (|3.18p may be generalized 
for the case when the coefficients in ()3.ip are functions of z. 
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Remark. Deriving formula (|3.18|) we have assumed that the function f{z) is meromorphic in 
the half-plane Rez ^ (except possibly at z = 0). However this formula is valid also for some 
functions having branch points on the imaginary axis, for example, 



/(^)=/i(^)n(^'+^/) 



(3.19) 



1=1 



with meromorphic function fi{z). The proof for (j3.18p in this case is similar to the one given 
above with the difference that branch points ibic; have to be circled on the right along contours 
with small radii. In view of further applications to the Casimir effect let us consider the case 
k = 1. By taking into account that 



(x^e 

from ()3.18p one obtains 



g±m/2^ - c 



if ^ X ^ c, 
'2 ii X > c, 



(3.20) 



k=l 



hm I 2 V T,(A,,fe)/(A,,fe) + n,[f{z)] - p.v. T dx f{x) I = ^Res/( 

7i(xe"^/2)_e^-Vi(^e-"*/2 



1 



dx 



dx 



K,{x) r 



(3.21) 



where f{z) = /i(2:) (2;^ + c^)^"*^^^ , c > 0. In Section [12] we apply this formula with an analytic 
function fi{z) to derive the Wightman function and the vacuum expectation values of the 
energy-momentum tensor in geometries with spherical boundaries. ■ 



For an analytic function f[z) formula (I3.18P yields 

'i^v,kf{K,k) 



A2. 



k=i y^7.,k - Ji{K,k) + XlkJu{X.,k) 

1 

vr 



r°° TT Y (z) 

/ dx fix) + -Res f{z)^ 
Jo 22=0 J^{z) 







/(xe"^/2) + e^"7(xe-"^/2> 



(3.22) 



By taking in this formula i' = 1/2, A = 1, B = (see notation (|3.ip ). as a particular case we 
immediately receive the APF in the form (j2.16p . In like manner, substituting = 1/2, A = 1, 
B = 2, we obtain APF in the form (j2.19p . Consequently, the formula (j3.18p is a generalization 
of the APF for general u (with restrictions given above) and for functions /(z) having poles in 
the right half-plane. 

For further applications to the Casimir effect in Sections [l2] and [T31 let us choose in ()3.22p 

f{z) = F{z)J^^^{zt), t>0, Rei/^0, (3.23) 
with m being an integer. Now the conditions (j3.4p formulated in terms of F(z) are of the form 



\F{z)\ < \z\ee 



{c-2t)\y\ 



or 



\Fiz)\ < , z = x + iy, \z\^ 00, (3.24) 
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with the same notations as in (|3.4p . In like manner, from condition (j3.1ip for F{z) one has 

F(ze"^) = -F{z) + o(z-2™-i)^ ^ ^ 0. (3.25) 

Now, as a consequence of ()3.22p we obtain that if the conditions (j3.24p and (j3.25p are satisfied, 
then for a function F{z) analytic in the right half-plane, the following formula takes place: 

k=i -^0 

-I Tdx ^j^ll^^ixt) [F(xe-/2) + F(xe— , (3.26) 
TT Jq I^{x) ^ L J 



for Re u ^ and Re + m ^ 0. 
3.2 Case (b) 

Let f{z) be a function satisfying the condition 

/(2;e"^/2) = _e2'^-7(2;e~"*/2)^ (3_27) 

for real x. It is clear that if f{z) has purely imaginary poles, then they are complex conjugate: 
=t*yfc) Vk > 0. By (|3.27p the rhs of p.lOh for a — > and argA^^fc = 7r/2 may be written as 



where p+ = 1, = 2, Cp{ak) denotes the right half of the circle with radius p and with the 
center at the point (Tfc, described in the positive sense, and the contours 7^ are the same as 
in (j3.13p . We have used the fact the purely imaginary zeros of Ju{z) are complex conjugate 
numbers, as Jy{ze'^'^) = e^'^'^Jy{z). We have also used the fact that on the right of (|3.10p the 
integrals (with a = 0) along straight segments of the upper and lower imaginary semiaxes are 
cancelled, as in accordance of (13.27P for arg z = 7r/2 

Let us show that from (j3.29p for zq = xqc'^'/^ it follows that this relation is valid for any z in a 
small enough region including this point. Namely, as the function f{z)H^\z)/Jy{z), p = 1,2 
is meromorphic near the point (— l)^+^xoe'^*/^, there exists a neighborhood of this point where 
this function is presented Laurent expansion 

~tW^^"^ = [.-(-l)-+ixoe-/2]- • 

From (j3.29p for z = xe'^*/^ one concludes 



00 



„(1) „-?i7ri/2 °° ( -i\nJ2) -n-Ki/2 

^ {x — Xq)"' ^ [x — Xq 



E \n , (3.31) 



n=— no n=—no 
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and hence On"* = (— l)"an^\ Our statement follows directly from here. By this it can be seen 
that 



E « / dz^^tli^m = 2vriRes4^/(-), (3-32) 



where (Jk = iyk, ^u,k, argA^^fc = tt/2. Now, by taking into account (j3.15p and letting p — > we 

get (sail]: 

Theorem 3. Let f{z) be a meromorphic function in the half-plane Rez ^ (except possibly at 
z = 0) with poles z^, Rez^ > and Hyk, yk > 0, k = 1,2, ... ^u,p)- If this function satisfies 
conditions and {3.21^ then 

hm I 2 V r,(A,,fc)/(A.,A:) + ri,[f{z)] - 

= -y E (2-SJRes/(z)4^, (3.33) 
where on the left < Re z^ < b, ReA^^„ < b < ReA,y^ri+i and ri^[f{z)] is defined by (3, 




Note that the residue terms in (|3.32p with ak = Xu,k, arg X^^k = 7r/2, are equal to 4:T^(X,^^k)fiXu,k) 
and are included in the first sum on the left of (|3.33p . 

Remark. Let ±iyk, yk > 0, and ±A,^^fc, argAj^^fc = 7r/2, be purely imaginary poles of function 
f{z) and purely imaginary zeros of Ju{z), correspondingly. Let function f{z) satisfy condition 

/(z) = -e2'^-/(ze--*) + o((z-afc)-i), z ^ a^, o^ = iy^,X,^^. (3.34) 

Now in the limit a — > the rhs of (|3.10p can be presented in the form p.28p plus integrals along 
the straight segments of the imaginary axis between the poles. Using the arguments similar those 
given above, we obtain relation (j3.32p with an additional contribution from the last term on the 
right of (j3.34p in the form J"^ (-o-j.) ^ ((-^ ~ (yk)~^^dz. In the limit p ^ the latter vanishes and 
the sum of the integrals along the straight segments of the imaginary axis gives the principal 
value of the integral on the right of (j3.18p . As a result formula (13.18P can be generalized for 
functions having purely imaginary poles and satisfying condition (j3.34p . For this on the right of 
(I3.18P instead of residue term we have to write the sum of residues from the rhs of ()3.33p and 
take the principal value of the integral on the right. The latter exists due to condition (I3.34p . ■ 

An interesting result can be obtained from (|3.33p . Let A^^^ be zeros of the function J^(z) + 
BizJ'^{z) with some real constants Ai and Bi. Let f{z) be an analytic function in the right 
half-plane satisfying condition (j3.27p and f{z) = o{z^) for z — > 0, where /3 = max(/3^, /3,y) (the 
definition see ()3.12p ). For this function we get from ()3.33p : 

oo oo 

Y.^t^(^%)fi>^^ii) = Y.^u(^-^k)f{X.,k), f^ = i^ + m. (3.35) 

k=l k=l 

For the case of Fourier-Bessel and Dini series this result is given in [M] . 

Let us consider some applications of formula (|3.33p to the special types of series. First we 
choose in this formula 

f{z) = Fi{z)J^{zt), t>0, (3.36) 
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where the function Fi{z) is meromorphic on the right half-plane and satisfies the conditions 

|Fi(z)| < ei(rE)e('=-*)lfl or |Fi(z)| < M|zr"ie(2-*)lf I , |z| ^ oo, (3.37) 
with c < 2, ai > 1/2, ei(x) = o{^/x) for x +00, and the condition 

From (j3.37p it follows that the integral p.v.J^ dx Fi{x) J fj_{xt) converges at the upper limit and, 
hence, in this case formula ()3.33p may be written in the form 

^Ty{\y^k)Fi{K,k)Jfi{K,kt) = -p-v. / dxFi{x)J^,{xt) - -n^ [Fi{z)J^{zt)] 
k=i ■'^ 

"T (2-5o,jResFi(z)J,(zt)^||M. (3.39) 

Note that the Ihs of this formula is known as Dini series of Bessel functions [63]. For example, 
it follows from here that for t < 1, Recr, Rez^ > — 1 one has 

^^J"'^^ Ja+u+l{^u,k)Ju{K,kt) = 7; [ —Ja+u+l{z)Ju{zt) 
k=l \k ^JO ^ 



= 2W + I) ' (^-^0) 

(for the value of the integral see, e.g., [M])- For B = this result is given in |66j . In a similar 
way, taking 1^1 = v + m, 



(z^ + zl) 



with Re ^ and Re + m ^ 0, from p.39p for a < 2 — t, Re cr > Re + m one finds 



^T,y{X,y^k)K,k"^ Jv + m{K,kt)— 



*^=i [Kk + 4 



^ JO (x2 + zf) ' 

= [-Zl) 7 — , ^, ,„ , a> t, (3.42) 

and the sum is zero when a < t. Here we have used the known value for the Sonine integral |64j . 
If an addition to (I3.37[) . (I3.38[) the function Fi{z) satisfies the conditions 

Fi(a;e^^/2) = -e^'^'Fiixe''''/'^) (3.43) 

and 

\Fi{z)\ < ei{x)e^'^\y\ or |Fi(z)| < Mlzr^ie*!^'!, |z| ^ 00, (3.44) 

then formula (j6.7p (see below) with B = may be applied to the integral on the right of ()3.39p . 
This gives ^27^281: 
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Corollary 1. Let F{z) be a meromorphic function in the half-plane Rez ^ (except possibly 
at z = 0) with poles z^, Rez^ > and Hyk, yk > ^u,i)- If F{z) satisfy condition 

F(xe"^/2) = (-l)'"+ie'^"^F(xe-"^/2)^ (3.45) 

with an integer m, and one of the inequalities 

\F{z)\ <ei{x)e''\y\ or |F(z)| < M|z|-"ie"ols'l, \z\ ^ 00, (3.46) 

with a < min(i, 2 — t) = uq, ei{x) = o{x^/'^), x —>■ +00, qi > 1/2, the following formula is valid 

00 

k=l rii,=0,iyk,z^. 



X 



Res I [j,+mizt)Y,{z) - Y,+Uzt)Mz)] B^X . (3.47) 



Recall that for the imaginary zeros Xu^k^ on the Ihs of (j3.47p the zeros with positive imaginary 
parts enter only. By using formula (j3.47p a number of Fourier-Bessel and Dini series can be 
summarized (see below). 

Remark. Formula ()3.47p may also be obtained by considering the integral 

- f dz W^{zt)H^^){z) - H';^Uzt)H^;^){z)\ Bp- (3.48) 

where Ch is a rectangle with vertices iti/i, bzizih, described in the positive sense (purely imaginary 
poles of F{z)/ Ju(z) and the origin are circled by semicircles in the right half-plane with small 
radii). This integral is equal to the sum of residues over the poles within (points z^, Xu,kj 
(Rezk, KeXp^k > 0)). On the other hand, it follows from (|3.45p that the integrals along the 
segments of the imaginary axis cancel each other. The sum of integrals along the conjugate 
semicircles give the sum of residues over purely imaginary poles in the upper half-plane. The 
integrals along the remaining three segments of Ch in accordance with (j3.46p approach to zero 
in the limit b, h ^ 00. Equating these expressions for ()3.48p . one immediately obtains the result 

dHZl). ■ 

From (I3.47|) . for t = 1, F{z) = Ji,{zx), m = 1, one obtains \13\ I64j 

00 ^ 

J^7;(A,,fc)J,(A,,fcx)J,+i(A,,fc) = y, 0^x<l. (3.49) 



k=l 



In a similar way, choosing m = 0, F{z) = zJy{zx)/ (z^ — a^), S = 0, we obtain the Kneser- 
Sommerfeld expansion 



V UK,kt)MX.,kx) ^ ^_Uax) _ YM)Ua)] , (3.50) 
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for ^ X ^ t ^ 1. In (j3.47p as a function F{z) one may choose, for example, the following 
functions 

n n 

1=1 1=1 

n 

Re < ^ Re /i; + n/2 + 2p + 3/2 - m - Sbao', (3-51) 
1=1 

n 

zf-^'^-^Ylil- Joibiz)], Reu <2n + 2p + 3/2- m-6bao; (3-52) 
1=1 

zP~^ n (^^ + ^"lY'^ {bi^z'^ + zi\ , W > O-half of an odd integer, 
1=1 ^ ^ 

n 

Re 1/ < - ^ /i, + n/2 + 2p + 3/2 - m - ^a^; (3.53) 

n 

^ JJzl^'l [J^^^k^{aiz)Y^^{biz) - y^,+fc,(a/z)J^,(6iz)] , A;/ - integer, 
=1 

n 

Rei/ < n + 2p + 3/2 - m - ^ |A;i| - 6a,ao, a = '^\ai - bi\ < oq; (3.54) 



z 

1=1 



1=1 



with p = + m — 2p (p is an integer), as well as any products between these functions and with 
Y\i (z^ — cf) ^\ provided condition (I3.46P is satisfied. For example, the following formulae take 
place 

k=l ■Ju+lKJv.k) 



< t 1, 



c = |a/ — 6/1 ^ t, ai,6/>0, Re//;^0, Re < n + 3/2 — Jcj; (3.55) 
1=1 

Jv{.3v,ki)Jv+i{^3v,k) Y\\^ J (h n \] - ^""^^ ~ ^'^"^ TT ^2 

n 

A + ^6i^t^l, A,6i>0; (3.56) 
1=1 

Jti{ju,kb)Ju+i{>^ju,k)Ju{ju,kt) 7rJ^+i{aX) , . j , t / w f +\] 



Ju{a) 



X + b^t^l, \,b>0, Re/i > -7/2 + 6x+b,t, (3.57) 



where jj^^^ are zeros of the function Jy{z). The examples of the series over zeros of Bessel 
functions we found in literature (see, e.g., |131 [Ml [66l I67j). when the corresponding sum was 
evaluated in finite terms, are particular cases of the formulae given in this section. 
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4 Summation formulae over zeros of combinations of cylinder 
functions 



Here we will consider series over the zeros of the function 



(4.1) 



where and in what follows for a given function f{z) the quantities with over bars are defined by 
the formula 

f j^z) = Ajf{x) + Bjzf'iz), j = a, b, (4.2) 

with constant coefficients Aj and Bj. This type of series arises in calculations of the vacuum 
expectation values in confined regions with boundaries of spherical and cylindrical form. To 
obtain a summation formula for these series, in the GAPF we substitute 



1 

2i 



' 0'\iiz) ^ 0'\iiz) 
hI"'\z) 0''\z)^ 



h{z) 



h{z) 



Hi^''\z)0''\z) 



(4.3) 



where for definiteness we shall assume that i] > 1. The sum and difference of these functions 
are given by the formula 



g{z) - (-l)V(^) 



.Hi'"'\Xz) h{z) 



k = 1,2. 



The condition for the GAPF written in terms of the function h(z) takes the form 

1/1(2)1 < ei(x)e^^'^' l^l ^ 00, z = x + iy, 



(4.4) 



(4.5) 



where ci < 2(r/ — 1) and x^^'^'-^'^^^bO ^ei(x) — > for x +00. Let ^y^^k be zeros of the function 
C'^^{ri, z) in the right half-plane. In this section we will assume values of Aj, and Bj for which 
all these zeros are real and simple, and the function hI^"'\z) {H^°'\z)) has no zeros in the right 
half of the upper (lower) half-plane. For real v and Aj, Bj the zeros ju,k Sire simple. To see 
this note that the function J^{tz)Yj"'\z) — Yi^(tz)jjf\z) is a cylinder function with respect to 
t. Using the standard result for indefinite integrals containing the product of any two cylinder 
functions (see jMlES]), it can be seen that 



dtt 



Mtz)Yj'^\z)-Y,itz)ji-\z) 



TT'^zTf{rj,z) 



where we have introduced the notation 



On the other hand 



j!>'^\vz) 



lci\,,z) 



2/^2 



K - Bt{z 



Combining the last two results we deduce that for real v, Aj, Bj the derivative 
and, hence, the zeros z = "y^^k are simple. By using this we can see that 

Res g{z) = ^Tf{r],-i,,k)h{-i,,k). 



(4.6) 



(4.7) 



(4.8) 



IS nonzero 



(4.9) 



19 



Hence, if the function h{z) is analytic in the half-plane Re z ^ a > except at the poles 
(7^ lu,i) and satisfies condition (j4.5p . the following formula takes place 



{ 2 rb 
-Y,K\Xn.,k)hhu,k)+r2Ah{z)]-p.v. 
I Ji) 



h{x)dx 



k=n 

'^+'°° 0^\\z) h{z) 
Hi^'Hz) Ci\\z) 



jr^(x) + yj'^)'(x) 



dz 



Hi^'\Xz) h{z) 
hI'^\z) Ci\\z^ 



-dz. 



(4.10) 



where we have assumed that the integral on the left exists. In formula (|4.10p . ")u,n-i < 0, < ^u,n, 
7u,m < b < 7i/,m+i, a < Ke Zk < b, Rez^ ^ Rez^+i, and the following notation is introduced 



V Res 

^-^ Im2i.=0 



+Tr > Res 



k,l=l,2 



(-l)'Im2fc<0 



t'^\z)+Y!;''^\z) 

Hi"'\Xz) hjz) 
m^\z) CfiKz) 



CfiX,z) 



(4.11) 



General formula (j4.10p is a direct consequence of the GAPF and will be used as a starting point 
for further applications in this section. In the limit a — > one has j271 EHJ [38] : 

Corollary 2. Let h{z) be an analytic function for Re z ^ except the poles Z]^ ), Re > 

(with possible branch point z = Q). If it satisfies condition |^.5D and 



and the integral 



hize""') = -h{z) +o{z-^), z^O, 
h{x)dx 



p.v. 



j('^)2^:E) + y, 



(4.12) 
(4.13) 



exists, then 



{2 rb 
y ^ h{^,^k)TfiX, 7.,fc) + r2Ahiz)] - p.v. / 
k=i ( 



h{x)dx 



—Res 
2 z=o 



h{z)0''\Xz) 
C-\X,z)Hi'''\z 



(x)+yj")'(x) 
dx il,au{x, Xx) h{xe^) + h{xe~^) 



(4.14) 



In formula ()4.14p we have introduced the notation 

^au{x, Xx) = i^^'^(Ax)/i^^"|(x) 

Ki''\x)W {Xx) - K^^\Xx)t\x) 

In the following we shall use this formula to derive the renormalized vacuum energy-momentum 
tensor for the region between two spherical and cylindrical surfaces. Note that ()4.14p may be 
generalized for functions h{z) with purely imaginary poles Hyk, Uk > 0, satisfying the condition 

h{ze-^') = -h{z) + o{{zTiykr^), z^±iyk. (4.16) 

The corresponding formula is obtained from ()4.14p by adding residue terms for z = iyt in the 
form of (j4.18p (see below) and taking the principal value of the integral on the right. The 
arguments here are similar to those for Remark after Theorem 3. 
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In a way similar to (j3.18p , one has another result [271 1 



Corollary 3. Let h{z) be a meromorphic function in the half-plane Rez ^ (with exception the 
possible branch point z = 0) with poles Zf., Hyk (7^ 7i/,jJ; Rezfc, > 0. If this function satisfies 
the condition 

/i(xe^*/2) ^ -h{xe-'''/^), (4.17) 

and the integral ^4-^3 ) exists, then 



{ 2 ^ pb 
Y,Klu,k)Tf{\lu,k) + r2u[h{z)]-^.w. / 



h{x)dx 



x) 



rik=0,iyk 

where on the left ^y^m < b < ^u,m+i ■ 



z=Vk 



Hi^'\Xz) h{z) 
Hi"'\z) Ci\\z) 



(4.18) 



Let us consider a special application of formula ()4.18p for Aj = 1, Bj = 0. Under the 
conditions given above, the generalizations of these results for general Aj,Bj are straightforward. 

Theorem 4. Let the function F[z) be meromorphic in the right half-plane Rez ^ (with the 
possible exception z = 0) with poles zp^^^iyj. (^ lu,i), Vk^^^^k > 0. If it satisfies condition 

F(xe"^/2) ^ (-l)™+iF(xe-"'/2), 



with an integer m, and one of two inequalities 

\F{z)\ < e{x)e''^\y\ or |F(z)| < Mlzr^e^^l^'l, 1^1 ^ 00, 
with ai < min(2A — a — 1, a — 1) = 02, o" > 0, e(x) for x +00, a > 1, then 

Yl 72/ ^wf2^7r^^ ^ 7 [J,y{l,y,k)Y^+m{(T-fu,k) - Y^{j^^k)J,y+m{(T7u,k)] 

fr[ Ju{iu,k)/J^[>'iu,k) - 1 



(4.19) 
(4.20) 



E 

Vk=o,iyk,Zk 



^QVk ^^^^ Yy{\z)Jy+jn{crz) - Jv{Xz)Yu+jn{(yz) ^^^^ 



vr 



z=rik Jy{z)Y^{\z) - J^{\z)Y^{z) 



(4.21) 



Proof. In (j4.18p let us choose 

h{z) = F{z) [Mz)Y^+m{(Jz) - Y^{z)J^+m((Jz)] , (4.22) 

which in virtue of (j4.20p satisfies condition (j4.5p . Condition ()4.17p is satisfied as well. Hence, 
h{z) satisfies the conditions for Corollary 3. The corresponding integral in ()4.18p with h{z) from 
(|4.22p can be evaluated by using formula (j8.7p (see below). Putting the value of this integral 
into (|4.18p after some manipulations we obtain formula (j4.2ip . ■ 

Remark. Formula (j4.2ip may also be derived by applying to the contour integral 

[ Y^{\z)Jy+n,{az) - Jy{\z)Yy+rn{(yz) 

Jc, Uz)Y,{Xz) - UXz)YAz) ^^^^^ ^^-'^^ 

the residue theorem, where Ch is a rectangle with vertices itz/i, b it ih. Here the proof is similar 
to that for Remark to Corollary 1. ■ 

By using ()4.18p a formula similar to ()4.2ip can also be obtained for series of the type 

T,kLiG{iu,k)Jf.iiu,kt). 

As a function F{z) in (I4.2ip one can choose, for example. 
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• function (j3.5ip for p = m — 2p, J2i h < a2, m < 2p + Re pi + n/2 + 1, p - integer; 

• function (|3.52p ioic p = m — 2p, J2i h < 02, m < 2p + 2n + 1; 

• function p.54p for p = m — 2p, ai > 0, Re/i/ ^ (for Re/U/ < 0, A;; > |Re///|), 
Ya=i \ai-hi\<a2,m<2p + n-Yji\h\ + '^- 

Taking F{z) = I/2, m = 0, one obtains 

£ — Ji{i.,k)/Ji{xi.,k)-i — =v^^^' (^-^^^ 

In a similar way it can be seen that 

^ 7^,fc ['/i^(7i.,fc)>"i'(<^7j/,fc) - yu{lu,k)Jv{'yiu,k)] 



ti (7^,,-c2)[J2(7.,fc)/J2(A7.,fc)-l] 

^ 1 Y^iXc)Mac) - MXc)Y^iac) 
~ TT Mc)Y^{Xc) - MXc)Y^ic) ' 



(4.25) 



J\T(Ji;TTy ^ = E^'<^-l'^«W + ^ + l> (4-26) 



where Rec ^ 0, 6; > 0, A ^ a > 1, pi 7^ —1,-2,.... Physical applications of the formulae 
derived in this section will be considered below. 



5 Summation formulae over the zeros of modified Bessel func- 
tions with an imaginary order 

5.1 Summation formula over the zeros of the function Kiz{rj) 
We denote by z = a;„ = coniil) the zeros of the function 

= AK,M + Br]d^K,M (5.1) 
in the half-plane Re z > for a given ry: 

K,^Jr])=0, n = l,2,... (5.2) 

A summation formula for the series over these zeros can be obtained by making use of the 
GAPF. For this, as functions f{z) and g{z) in formula (j2.1ip we choose 

/(.) = ^F(z) sinhvr., g{z) = ^^^^^^^U^^F(.), (5.3) 

with a function F{z) meromorphic in the right half-plane Rez > 0. For the sum and difference 
of these functions one finds 

g{z)±f{z) = 2F{z)^f^. (5.4) 
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By using the asymptotic formulae for the modified Bessel functions for large values of the index, 
the conditions for the GAPF can be written in terms of the function F{z) as follows: 

|F(2)| < e(|z|)e-™' (|z|/r/)^l^l , z = x + iy, x > 0, \z\ ^ oo, (5.5) 

where |z|e(|z|) — > when \z\ — > oo. Let us denote by ZF,k, A; = 1, 2, ... , the poles of the function 
F{z). By making use of the fact that the zeros cok are simple, one finds 



R[f{z),9{z)] = 2v^^ |E ^^1^1^, ""'^ \ , (5.6) 

with the notation 



dzKiz{v)\z=uJk 



sr^ Res 

fc,Im 2j?J;^0 



rK[f{z)] = }2 Res ^^fff(l^F(z) 



+ Y. R.S ^''!^ff"" F(.), (5.7) 



and the function a{z) is defined in (j2.2p . Substituting (j5.6p into (|2.1ip . for the function -F(2;) 
meromorphic in the half-plane Rez ^ a and satisfying the condition ()5.5p one obtains the 
following formula 

( f; ^^f^^f^ - ^p.v. rdxF(x)sinhvrx + r;,[/(z)]| 



1 

"2Tii 



""°^..F(.)^-/™..F(.)i4^ 



(5.8) 



where Rewm-i < a < Rewm, Rea;„ < & < Recj^+i, and in the expression for rK[f{z)\ the 
summation goes over the poles with a < Re ZF,k < b. Below we consider this formula in the 
limit a ^ 0, assuming that the function F{z) is meromorphic in the right half-plane Re z > 0. 

From the relation K-iz{rj) = Kiz{rj) it follows that the purely imaginary zeros of the function 
Kiz{r]) are complex conjugate. We will denote them by ^it^'f^, ^^'^ > 0, and will assume that 
possible purely imaginary poles of the function F{z) are also complex conjugate. Let us denote 
the latter by ±iyk, Uk > and assume that near these poles F{z) satisfies the condition 

F{zen = -F{z) + o{{z-iyk)-^), z ^ iyk. (5.9) 

In the limit a — > the rhs of formula (|5.8p is presented as the sum of integrals along the straight 
segments of the imaginary axis between the purely imaginary poles plus the terms 

^"^^dz, (5.10) 




K^ziv) 

where the contours 7^ and Cp{ak) are the same as in (j3.28p . Assuming the relation 

F{ze'^) = -F{z)+o{l), z^iu'k, (5.11) 

in the limit p — > the integrals in (jS.lOp are expressed in terms of the residues at the corre- 
sponding poles. Noting also that in this limit the sum of integrals over the straight segments of 
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the imaginary axis gives the principal value of the integral, one obtains 



lim 

6— »+oo 



E 



\Z=U)U,%U), 



1 



p.v. / dxF{x) smhiTX + rK[f{z)] 



1 



p.v. / ax -= — — 



'2Ko{v) 



y Res F{ 

f ^ 7=1111- 



z=Wk 



F{xe 2 ) + F{xe 



(5.12) 



where Fq = \\mz^QzF{z). The application of this formula to physical problems will be given 
below in Section [THl Taking in formula (j5.12p F{z) = Kiz{ri)G{z), we obtain a formula relating 
two types of integrals with the integration over the index of the modified Bessel functions. In 
particular, from the latter formula it follows that for an integer m one has 



dx x^^ sinh nxKix {vj) 



cosh ax 
sinh hx 



mO 



2h 



k=i 



2m 



X cos I ^ ) I^k/biv), 



(5.13) 



under the condition |Re6| — |Rea| > 7r/2. The special case of this formula with B = and 
m = is given in [66]. 



5.2 Summation formula over the zeros of the function Ziz{u,v) 

In this subsection we will derive a summation formula over the zeros z = J7„ of the function 



ZUu,v)=kI:\u)iI^\v) 



(5.14) 



where we use the notation defined by formula (14. 2p . As we will see in the applications (see 
Section [T9|) . the vacuum expectation values of physical observables in the region between two 
plates uniformly accelerated through the Pulling-Rindler vacuum are expressed in terms of series 
over these zeros. To derive a summation formula for this type of series, we choose in the GAPF 
the functions 



2i 



f{z) = — sinh7rzi^(z 
vr 

fib) 

g[z' - 



m\vyi^l{^)+it\u)i'^l{v) 



(5.15) 



Ziz{u,v) 



F{z), 



with a meromorphic function F{z) having poles z = Zk ^n), Im-2fc 7^ 0, in the right half-plane 
Re z > 0. The zeros 0„ are simple poles of the function g{z). By taking into account the relation 



9iz) ± f{z) 



24t(n)7 



(b) , 
±iz\ 



Ziz{u,v) 



-F{z), 



(5.16) 



for the function R[f(z),g{z)] in the GAPF one obtains 



2TTi 



d^Zi^{u,v)U=Q„ 



Ziz{u,v) Mzk)z^ ' -Mzk)z 



(a) 



(5.17) 
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where the zeros 0„ are arranged in ascending order, and the function cr{z) is defined in formula 
(j2.2p . Substituting relations (|5.16p and (j5.17p into the GAPF, as a special case the following 
summation formula is obtained [49\ [5U] : 

Jq 27rZx[u,v) 
_y Res-£l£l-j(^; , (t;)/^'^.^ , (u). (5.18) 

Here the condition for the function F{^z) is obtained from the corresponding condition in the 
GAPF by using the asymptotic formulae for the modified Bessel function and has the form 

|F(z)| < e(|z|)e-^^ , x>^,z = x^iy, (5.19) 

for 1^1 ^ oo, where |z|e(|z|) when \z\ oo. Formula (15.180 can be generalized for the case 
when the function F(z) has real poles, under the assumption that the first integral on the right 
of this formula converges in the sense of the principal value. In this case the first integral on 
the right of Eq. (j5.18p is understood in the sense of the principal value and residue terms from 
real poles in the form Resg{z), Imzk = 0, have to be added to the right-hand side of this 

formula, with g{z) from (j5.15p . 

6 Applications to integrals involving Bessel functions 

The applications of the GAPF to infinite integrals involving combinations of Bessel functions 
lead to interesting results [27\ [28] . First of all one can express integrals over Bessel functions 
through the integrals involving modified Bessel functions. Let us substitute in formula ()3.18p 

f{z) = F{z)Mz). (6.1) 

For the function F{z) having no poles at z = X^^k the sum over zeros of Jy{z) is zero. The 
conditions ()3.4p and ()3.1ip may be written in terms of -F(2;) as 

\F{z)\ < ei{x)e^'\y\ or \F{z)\ < M\z\-''^ e^^^ |z| ^ oo, (6.2) 

with ci < 1, x^/^~''^oei(x) ^ for x ^ cxd, qi > ao = 3/2 — 5bq, and 

F(ze^*) = -e^^'Fiz) + o{z\^^''\~^) . (6.3) 

Hence, for the function F{z) satisfying conditions ()6.2p and ()6.3p . from ()3.18p it follows that 

p.v. / dxF{x)Ju{x) = riu \F{z)J4z)] + -ResF{z)Y4z) 

Jo ^ ' 2 z=0 

+ - / dxi^^(x) e~^"^/2^(xe"*/2) + e^"/2^(xe-"'/2) . (6.4) 
Jo L J 

In expression (j3.9p for ri^ [/(z)], the points are poles of the meromorphic function F{z), Kezk > 
0. On the basis of Remark after Theorem 3, formula (j6.4p may be generalized for the functions 
F{z) with purely imaginary poles ±iyk, Uk > and satisfying the condition 

F(ze-^) = -e^-*F(z)+o((zT%fe)-i), z ^ ±iyk. (6.5) 
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The corresponding formula is obtained from (|6.4p by adding residue terms for z = iy^ in the 
form of (j6.7p (see below) and taking the principal value of the integral on the right. 

The same substitution of (|6.1|) into formula (|3.33p . with the function F{z) satisfying condi- 
tions (j6.2p and 

F{xe'''l'^) = -e''^'F{xe-^'l^), (6.6) 
for real x, yields the following result 

p.v. / dxF{x)J,{x) = [F{z)Uz)] + - V (2 - ResF(z)^(i)(z). (6.7) 

Vk=0,Wk 

In ()3.9I) . the summation goes over the poles Zk, Rezfc > of the meromorphic function F{z), 
and ±iyk, Uk > are purely imaginary poles of this function. Recall that the possible real poles 
of F{z) are such, that integral on the left of ()6.7[) exists. 

For the functions F{z) = z'^^^F(z), with F{z) being analytic in the right half-plane and 
even along the imaginary axis, F(ix) = F{—ix), one obtains 

oo 

dxx''+'^F{x)Ju{x) = 0. (6.8) 





This result for B = (see (jS.ip ) has been given previously in ^68j. Another result of [68j is 
obtained from ()6.7p choosing F{z) = z'^'^^ F{z)/{z'^ — a^). 

Formulae similar to (j6.4p and (j6.7p can be derived for the Neumann function Yu(z). Let for 
the function F{z) the integral p.v.J^ dx F{x)Y,^{x) exists. We substitute in formula (12. lip 

f{z) = Y,{z)F{z), g{z) = -iMz)F{z), (6.9) 

and consider the limit a +0. The terms containing residues may be presented in the form 

R[f(z),g(z)] = Res //(^)(z)F(z) + vr V Res Hi^\z)F(z) 

+7ry Res Mz)F{z)=r34F{z)], (6.10) 

^ lmzfe=0 
k 

where z^ (Rez^ > 0) are the poles of F{z) in the right half-plane. Now the following results can 
be proved by using ()2.1ip : 

1) If the meromorphic function F{z) has no poles on the imaginary axis and satisfies condi- 
tion ^6.2\) . then 

POO ^ POO 

p.v./ dxF{x)Y^ix) = r3^[F{z)]-- dx K^{x) 

Jo Jq 

and 

2) If the meromorphic function F{z) satisfies the conditions 

F(xe"*/2) = e^"*F(xe-"^/2) (6.12) 

and i6.2\) . then one has 

POO 

p.v./ dxF{x)Y^{x) =r3^[F{z)]+ ttJ" Res hI^\z)F{z), (6.13) 
Jo , ^=*yfe 
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where =biy^., > are purely imaginary poles of F{z). 

From (|6.13p it directly follows that for F{z) = z'^F{z), with F{z) being even along the 
imaginary axis, F{ix) = F[—ix), and analytic in the right half-plane one has |68| : 

l-OO 

/ dxx^F(x)y^(x) = 0, (6.14) 

if condition (j6.2p takes place. 

Let us consider more general case. Let the function F{z) satisfy the condition 

F{ze-^') = -e-^^^F(z), (6.15) 

for argz = 7r/2. In the GAPF as functions f{z) and g{z) we choose 

f{z) = F(z)[J^(z)cos(5 + y^(z)sin5], 

g{z) = -iF{z)[JAz)s].n5 -Yy{z)cos5], 6 = {\-v)tt/2, (6.16) 

with g{z) — {—l)^f{z) = h\!'\z)F{z) exp[(— l)*''i(5], k = 1,2. It can be seen that for such a choice 
the integral on the rhs of (I2.11|) for a — > is equal to 

y ResFW(^)F(z)e-*^ (6.17) 



TTl 

'nk=Wk 



where ±iyk, Vk > 0, as above, are purely imaginary poles of F{z). Substituting (j6.16p into (|2.1ip 
and using (f2^ we obtain [271128] : 

Corollary 4. Let F{z) he a meromorphic function for Rez ^ (except possibly at z = 0) with 
poles Zk, iiyk! y/oR-ezfc > 0. // this function satisfies conditions \6.2\) (for B = 0) and ^6.15]) 
then 



v.f dxF{x)[J^{x)cos5 + Y^{x)sin6]=T^i\Y^ Res H^^\z)F{z)e-'^ 

Jo ^lmzfe>0 

-V Res Hi^\z)F(z)e'^ -iy" Res \JJz) sin 6 - YJz) cos 6] F(z) 

Zk Zk 

+ y ResF«(z)F(z)e-^4, (6.18) 

— ^ z=rik 

Vk=iyk ) 



Vk=Wk 

where it is assumed that the integral on the left exists. 



In particular for 6 = vrn, n = 0, 1, 2..., formula (j6.7p follows from here in the case B = 0. One 
will find many particular cases of formulae (16. 7p and ()6.18p looking at the standard books and 
tables of known integrals with Bessel functions (see, e.g., [I3l [Ml [STJ [661 [Ml [ZOl [ZD [22]). Some 
special examples are given in the next section. 



7 Integrals involving Bessel functions: Illustrations of general 
formulae 

In order to illustrate the applications of general formulae from the previous section, first we 
consider integrals involving the function Ju{z). Let us introduce the functional 

POO 

A,^[Giz)]=p.v. dzz^-^"'-^G{z)Uz), (7.1) 
Jo 
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where m is an integer. Let -^1(2:) be an analytic function in the right half-plane satisfying the 
condition 

Fi(xe^^/2) = Fi(xe-^*/2)^ ^^(g) / 0, (7.2) 

(the case when -^1(2;) ~ z'^, z ^ 0, with an integer q, can be reduced to this one by redefinitions 
of Fi{z) and m). Prom (j6.7p the following results can be obtained [271 [28] : 



A, 
A. 



Fi{z) 
— 



4(2m) 



vr 



-a- -y,(a)Fi(a) + 40^ 



z^-a^ 



a 



v—2m—A 



-Y,ia)Fi{a) - (-l)'"i^,(a)Fi(m) 



+ 



I 4(0) 



A, 



F(z) 



(z2 



z^-a^ 



d 



(7.3) 
(7.4) 

(7.5) 



+^(0) ^1^^ 



Fimi'Hc) 



fz2 _ c2)P+l 



(7.6) 



Fiiz) 



{z"^ + a 



2\P+1 



(_l)m+p+l / 



2P-p! Vot^a 



P r 



,!/-2m-2 



ir,(a)Fi(ae'^*/2) 



I 4(0) 



Fi{z) 



^2 + a2)^'+^ 



(7.7) 



and etc. Note that A^^m = for m < 0. In these formulae a > 0, < argc < 7r/2, and we have 
assumed that Rez^ > 0. To secure convergence at the origin the condition Rei/ > m should be 
satisfied. In the last two formulae we have used the identity 



dz 



zF{z) 




_(z + 6)P+i_ 





d 



Fib). 



(7.8) 



Note that (17. 7p can also be obtained from ()7.6p in the limit Rec — > 0. For the case Fi = 1, m 
— 1 in (j7.7p see, for example, [M]. In (|7.3p - (|7.7p as a function ^1(2) we can choose: 



• function ([331]) for p = 1, Reu < Y^^Re ni + 2m+{n + l)/2-6bi + 6Bo, b = J^k ^ k > 0; 

• function ([332]) with p = 1, Rez^ < 2(m + n) + 1/2 - 6bi + Sbo, 6 = 6« ^ 1; 

• function ([3.53p for p = 1, Re < 2m — ^ Rep; + (n + l)/2 — 5bi + Sbo, p; > is half of 
an odd integer, 6 = ^6/^1; 

• function p.54p for Re < 2m + n — J2\h\ + ^/"^ + ^BO — ^ai, d = Y^\ai — bi\ ^ 1, ai ^0, 
ki - integer. 

Here we have written the conditions for (|7.3p . The corresponding ones for (|7.4p . (|7.5p . (|7.6p . ([7.7p 
are obtained by adding on the rhs of inequalities for Rei^, respectively 2, 4, 2(p + 1), 2(p + 1). 
In ()7.3p - (l7.7p we can also choose any combinations of functions ()3.5ip - ([3.54p with appropriate 
conditions. 



For the evaluation of A^m in special cases the following formula is useful 

2m / 1 \ m 



lim 



d 



V dz 



fi{z) = (2m - 1)!! lim 



d 



V zdz 



(7.9) 
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valid for the function fi{z) satisfying condition fi{—z) = fi{z) + 0(2;^™), z ^ 0. From here, for 
instance, it follows that for z ^ 



2m 



^ ) [z'^Y,{bz)F,{z)] 



-{2m - 1)!! 



k=0 



T{u-m + k) / d \ 



zdz I 



Fi{z), 



(7.10) 



where we have used the standard formula for the derivative {d/ zdz)"^ of cylinder functions (see 
[65]). From ()7.3p one obtains {B = 0) 



/ dzz'^-^^-\uz)ii{z' + zfy''^'j,,i 

71 



-vr 



zdz 



1=1 

n 



Y^{z)J{{z' + z^)-^^''j,,{hJz^ + z^) 



1=1 



(7.11) 



z=Q 



for m ^ and the integral is zero for m < 0. Here Rez^ > 0, 6 = 6; ^ 1, 6^ > 0, 
m < Rei^ < Y^^=i^^ + + (n + 3)/2 — 81,1. In the particular case m = we ob- 
tain the Gegenbauer integral [131 I64j . In the limit zi — > 0, from (j7.1ip the value of integral 



Juiz) nr=i ^^Jfj^ii^) is obtained. 



By using (j6.18p . formulae similar to (j7.3p - ()7.7p may be derived for the integrals of type 

poo 

B^[G{z)]=p.v. dz G{z) [J^{z) cos 6 + Y^{z) sin 6], 5 = {\-v)t^/2. (7.12) 
Jo 

It directly follows from Corollary 4 that for function F{z) analytic for Re 2; ^ and satisfying 
conditions (|6.2p and (j6.15p the following formulae take place 



Bu 



Bu[F{z)] 

z"^ — 

F{z) 
z^-a^ 







7rF(a) [Juio) sin 5 — Yy{a) cos 5\ /2, 



vr 

^-gF(a) [Ju{a.) sin 5 — Yjj{a) cos 5] 



B„ 



B,. 



Fi{z) 



(z2 - c2)f+l 
^l(^) 



d 



(z2 + a2) 



2P+1 .pi \cdc 

(-1)P+^ / d 
2P -pi [jida 



c-'F{c)Hl^\c) 
a-^F{ae''''/^)K^{a) 



-inX/2 



(7.13) 
(7.14) 

(7.15) 
(7.16) 

(7.17) 



where a > 0, < argc ^ 7r/2. To obtain the last two formulae we have used identity 
Formula (j7.13p generalizes the result of [6S\ (the cases X = v and A = u + l). Taking F{z) = 
from the last formula we obtain the result given in [64j. In (j7.13p - ()7.17p . as a function F{z) one 
can choose (the constraints on the parameters are written for formula (|7.13p : the corresponding 
constraints for (j7.14p . (j7.15p . (j7.16p . (|7.17p are obtained from the ones given by adding 2, 4, 
2{p + 1), 2(p + 1) to the rhs of inequalities for Re 1^, correspondingly): 



function p.5ip for p = A, |Rez^| < Rep < ^ Re//; + (n + 3)/2 — (J^i, 6 = 6/ ^ 1; 
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• function (jc{.52p for p = A, |Rez^| < Rep < 3/2 — ^51, 6 = ^^6/ ^ 1; 

• function 



1=1 



1=1 

ai > 0, \Reu\ + ^ \ki\ <Rep <n + 3/2 - 6ci, 
c = '^\ai- bi\ ^ 1, Rem ^ 0, 

(for Re /// < one has ki > | Re p; | ) . 



(7.18) 



Any combination of these functions with appropriate conditions on parameters can be chosen 
as well. 

Now consider integrals which can be expressed via series by using formulae ()6.7p and ()6.18p . 
In (16.711 let us choose the function 



F{z) 



sinh TTZ 



(7.19) 



where Fi{z) is the same as in formulae ()7.3p - ()7.6p . As the points ±i,±2i, . . . are simple poles 
for F(z), from (16.71) one obtains 



°° , z^-^'^JJz) 
az 



sinh (ttz) 



zFi{z) 

sinh(7r2;) J ' vr 



TT ^ ^ 



Ky{k)Fi{ik), (7.20) 



where 74Sri[/(2:)] is defined by (j7.3p and Reu > m. The corresponding constraints on Fi{z) 
follow directly from (j6.2p . The particular case of this formula when -^1(2;) = sinh(az)/2; and 
m = — 1 is given in [64J. As a function Fi{z) here one can choose any of functions (I3.51l )- (l3.54p 
with p = 1 and a, Yli^i < 1- From (|7.20p it follows that 



dz 



sinh (vrz) 



Mz)llz-^^I,,{kz) = A^^l 



1=1 



sinh(7rz^ 



1=1 



+-J2^-'i-r^'K,ik)l[k-^'^J^^ibik), bi>0,Y,bi< TT, Reu > 

k=l 1=1 1=1 

In a similar way the following formula can be derived from (j6.18p : 



m. 



(7.21) 



00 7F(7) 2? 

dz . \ [J^(z)cos(5 + y^(z)sin(5] = -^e-'^''/^^{-l)''kK^{k)F{ik), (7.22) 

k=l 



sinh(7rz) 



TT 



where A is defined by relation (j6.15p . The constraints on the function F{z) immediately follow 
from Corollary 4. Instead of this function we can choose functions (j3.5ip . (j3.52p . p.54p . 

As it has been mentioned above, adding the residue terms tti Res F[z)hI^\z) to the rhs of 

z=iyk 

(16. 4p this formula may be generalized for functions having purely imaginary poles Hyk, Uk > 0, 
provided condition (|6.5p is satisfied. As an application let us choose 



F{z) 



z^F^jz) 



, Fi{-z) = Fi{z), 6>0, 



(7.23) 
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with an analytic function Fi{z). Function (j7.23p satisfies condition (|6.5p and has poles zizikb, 
k = 0,1,2 .. .. The additional constraint directly follows from (|6.2p . Then one obtains 



/ dx , ' Fi{x) = - Y {bkYK4bk)Fi{ibk) - - / dxx''K^{x)Fi{ix), (7.24) 
Jq e^-^^l^ - 1 vr ^ TT Jo 

where the prime indicates that the m = term is to be halved. For the particular case F\(^z) = 1, 
using the relation 



oo . oo 

Y^{bkYK,{bk) = ^lyViv + 1/2) [(27rn/6)2 + l] 



-v~\l2 



(7.25) 



fc=0 



n=0 



and the known value for the integral on the right, we immediately obtain the result given in 
[M] . Relation (j7.25p can be proved by using the formulae (for the integral representation of 
McDonald function see p3] ) 



vrz" 



+ 00 



Ku{z) = / , e^'=^ = 27r ^ 5(z-27rn), (7.26) 

k=—oo n=-oo 



Jo (i^ + l) 
where (5(z) is the Dirac delta function 



8 Formulae for integrals involving Jjy{z)Y^{Xz) — Yi,{z)J^{Xz) 

In this section we will consider applications of the GAPF to integrals involving the function 
Ju{z)Yf^{Xz) — Yu{z)J^{Xz). In formula (12. lip we substitute 



1=1 



2i 



1=1 h2\z) 



For definiteness we consider the case A > 1 (for A < 1 the expression for g{z) has to be chosen 
with opposite sign). Conditions (j2.ip and ()2.10p are satisfied if the function F{z) is constrained 
by one of the following two inequalities 



\F{z)\ < eix)e 



c\y\ 



c < A — 1, e{x) ^ 0, X ^ +0O, 



or 



|F(z)| < M|z|-"e(^-^)l^l, a> 1, \z\^ oc, z = x + iy. (8. 
Then, from ()2.1ip it follows that for the function F{z) meromorphic in Rez ^ a > one has 

Mx)Y^{Xx) - J^{Xx)Y,ix 



■2) 
.3) 



nOO 

p.v. / dx ■ 

/o 



+- 



1 

2i 



Jd{x) + Y^{x 
dzF{z) 



Hi'\z) 



-F{x) = nf,^[F{z)] 

r(2) 



dzF(z] 



HPiz) , 



(8.4) 



where we have introduced the notation 



nf.u[F{z)] 



V Res 

k 

2 

+7r "S^ "S^ Res 

. (-l)'lm^fe<0 



1=1 hI'\z) 



F{z) 



hI!\xz) 

Hli\z) , 



^.5) 
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The most important case for the apphcations is the Umit a — > 0. The fohowing statements take 
place [271128]: 

Theorem 5. Let the function F{z) he meromorphic for Kez ^ (except the possible branch 
point z = 0) with poles z^, Hyk (Vk^^^Zk > 0). If this function satisfies conditions 18. or 
and 

then, for values of v for which the function Hi^\z) (Hj^\z) ) has no zeros for ^ argz ^ 7r/2 
(—tt/2 ^ arg z ^Q), the following formula is valid 



oo 



J,{x)Y^{\x) - Y,{x)J^{\x) 



P.v.y^ dx j2(,)+y^2(,) = ri,.[F(z)] + 



+ f E (8.7) 



where it is assumed that the integral on the left exists. 

Proof. From condition (j8.6p it fohows that for argz = 7r/2 one has 

H}^\z) H)f\ze-^'^) 

and that the possible purely imaginary poles of F{z) are conjugate: izy^, Ut > 0. Hence, on 
the rhs of (j8.4p . in the limit a ^ the term in the square brackets may be presented in the form 
(it can be seen in a way similar to that used for (I3.32p ) 

E ( [ +E / ]dz^^^^^Fiz), (8.9) 

with the same notations as in (j3.32p . By using (jS.Sp and the condition that the integral converges 
at the origin, we obtain 



E / dz ^\^\^'^ F{z) = (2-^o.JvrzRes4!y^^(^). (8-10) 



where (0) = 7^, Q^{iyk) = Cp{±iyk). By using this relation, from (j8.4p we receive formula 



Note that one can also write the residue at z = in the form 



Integrals of type (j8.7p which we have found in literature (see, e.g., [131 [661 EO]) are special cases 
of this formula. For example, taking F{z) = Ji,{z)Yu-\-i{X' z) — Yi,{z)Ju+i{X' z) for the integral 
on the left in (jHZD we obtain -X-^X'-"'^ for A' < A and X^X'^"-^ - X-^X'-"-^ for A' > A (see 



32 



[13]). By taking z^'^^^ l{z^ + a^), — (?) as for fi = u and integer m ^ one 

obtains 

°- J.(x)y.(Ax) - y.(x) J.(Ax) x^'^+i _ ...^ o^TT^MAo) 

/"^ J^(x)y^(Ax) - y^(x) J^(AX) X^'^+l ^ TT 2^ 

^ y.(c)j.(Ac) + y.(c)y.(Ac) 

where Re a > 0, c > 0, A > 1. Special cases of this formula for = m = are given in [13] . 
In ()8.12p taking the limit a ^ and choosing m = 0, we obtain the integral of this type given 
in [66j. In (j8.7p as a function F{z) we can choose (|3.5ip . (|3.52p . (|3.54p (the corresponding 
conditions for parameters directly follow from (18. 2p or (18. 3p ) with p = fj, — v — 2m [m is an 
integer), as well as any products between them and with H/LiC-^^^ ~ For instance, 

dx J^{X)Y^{\X) - y^(x)J^(Ax) -A- Jf^li^l^/^^ + Zf) _ TT 
~ J2(x)+y,2(x) (x2+z2)/^^/2 " 2A^ 

n 

x[|z-'^'J^^(6,^0, Rez^> 0, Rez, ^ 0, A > 1, 
1=1 

n n 

Re w + n/2 + 1 > S^^x-i, 6 = 6/ < A - 1. (8.14) 
1=1 1=1 

As another consequence of formula ()8.4p one has: 

Theorem 6. Lei -F(2;) 6e meromorphic in the right half-plane (with possible exception z = 0) 
with poles Zk, Rez^ > 0, and satisfy conditions ^8. ^|) or \8.3\) and 

F{ze'^') = -e^^'-'^^^'Fiz) + o(zlR«'^l-|RH-i)^ ^ ^ (g_-^5) 

then for values of v for which the function hI^\z) (hI^\z) ) has no zeros for ^ argz ^ 7r/2 
(—tt/2 ^ arg z ^Q) the following formula takes place 

r , JAx)Y^{>^x)-Y,{x)J^{\x) ^, TT 4'^(Az) 

P"7o J.^(x) + y,^(x) = + 2S^a)^^(^) 



{z.-^)^V2^(a;e'ri/2) ^ g(/.-i.)^i/2^(^g-^i/2)l ^ ^ > 1^ (8.16) 



2 70 K,{x) 
provided the integral on the left exists. 

Proof. This result immediately follows from ()8.4p in the limit a ^ and from (j8.10p with 
% = 0. ■ 



For example, by using (j8.16p one obtains 







n 

xJJ/^,(6,x), ^Rew + |Rei/| > |Re//| - 1, 6 = 6/ ^ A - 1, 6^ > 0, 

n 

n > 4,A-1, jUs = - ^ + y^/i;. (8.17) 
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Such relations are useful in numerical calculations of the integrals on the left as the integrand 
on the right at infinity goes to zero exponentially fast. 

We have considered formulae for integrals containing J,y{z)Yf^{Xz) — Yi^{z)J^{Xz). Simi- 
lar results can be obtained for integrals with the functions Jl{z)Yfj_{Xz) — Yl{z)J^{\z) and 

jI{z)y;^{\z)-yi{z)j'^{\z). 

9 Wightman function, the vacuum expectation values of the 
field square and the energy-momentum tensor on manifolds 
with boundaries 

In this and following sections we will consider applications of the GAPF to the problems in 
quantum field theory with boundaries, mainly to the Casimir effect. The Casimir effect is one of 
the most interesting macroscopic manifestations of the non-trivial structure of the vacuum state 
in quantum field theory [2l [3l [H O [6l . The effect is a phenomenon common to all systems 
characterized by fluctuating quantities and results from changes in the vacuum fluctuations of 
a quantum field that occur because of the imposition of boundary conditions or the choice of 
topology. It may have important implications on all scales, from cosmological to subnuclear, 
and has become in recent decades an increasingly popular topic in quantum field theory. In 
particular, the recent measurements of the Casimir forces between macroscopic bodies (see, for 
instance, [6l[73]) provide a sensitive test for constraining the parameters of long-range interac- 
tions predicted by modern unification theories of fundamental interactions [74|- In addition to 
its fundamental interest the Casimir effect also plays an important role in the fabrication and 
operation of nano- and micro-scale mechanical systems [6J. 

Before to consider special problems we first give general introduction for the procedure of 
the evaluation of the vacuum expectation values (VEVs) of the functions bilinear in the field on 
manifolds with boundaries. Consider a real scalar field 99 with curvature coupling parameter 
on a (-D -|- l)-dimensional background spacetime M with boundary dM described by the metric 
tensor gik- The corresponding field equation in the bulk has the form 

(ViV^ + + Ci?) (/J = 0, (9.1) 

where R is the scalar curvature for the background spacetime, m is the mass for the field quanta, 
Vj the covariant derivative operator associated with the metric gik- On the boundary of the 
manifold the field satisfies some prescribed boundary condition F[ip[x)\\x<^dM = 0- Here and 
below we adopt the conventions of Birrell and Davies [75] for the metric signature and the 
curvature tensor. The values of the curvature coupling parameter <^ = and C = Cd with 
C,£) = [D — 1)/4:D correspond to the most important special cases of the minimal and conformal 
couplings. It is convenient for later use to write the corresponding metric energy-momentum 
tensor (EMT) in the form 



(9.2) 



with Rik being the Ricci tensor. This form of the EMT differs from the standard one by the term 
(pgik (V/V^ -|- -|- C-R) V^/2 which vanishes on solutions of the field equation (see, for instance, 
|76j ) . It can be easily checked that the following relation for the trace of the EMT takes place: 

T^ = D{C-CD)ViV'ip^ + m^ip\ (9.3) 



Note that on manifolds with boundaries the EMT in addition to bulk part (j9.2p contains a con- 
tribution located on the boundary. For arbitrary bulk and boundary geometries the expression 
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for the surface EMT is derived in |76| by using the standard variational procedure for the action 
with boundary term. For solutions of Eq. (j9.ip . ipi and ip2, the scalar product is defined by the 
relation 

{fi,'f2) = -i dY. ^/-gT.n'- [fidnpl - (diipi) ip^] , (9.4) 



where dS is the volume element on a given spacelike hypersurface S, n* is the timelike unit 
vector normal to this hypersurface, and qy, is the determinant of the induced metric. The 
quantization of the field can be done by using the standard canonical quantization procedure 
(see, for instance, [H Eg EZl [78] ) . 

Let {ifa{x),ip^{x)} is a complete orthonormal set of classical solutions to the field equation 
satisfying the boundary condition, 

{^a{x),(Pa'{x)) = 6^a', F[ip^{x)]\j:(zdM = 0, (9.5) 

where the collective index a can contain both discrete and continuous components. On the right 
of the orthonormalization condition, Sua' is understood as the Kronecker symbol for discrete 
components and the Dirac delta function for continuous ones. In quantum field theory we 
expand the field operator in terms of the complete set: 

fix) = ^K'/'<t(2;) + a^(fl{x)], (9.6) 

with annihilation and creation operators a^, a^' satisfying the commutation relation [oo-,a^;] = 
(^o-cr'- The vacuum state |0) is defined by the relation a^lO) = for any a. Our main interest 
below will be the VEVs of the field square and the EMT. These functions give comprehensive 
insight into vacuum fluctuations and are among the most important quantities characterizing 
the properties of the quantum vacuum. Though the corresponding operators are local, due to 
the global nature of the vacuum, the VEVs describe both local and global properties of the bulk 
and carry an important information about the background geometry. In addition to describing 
the physical structure of the quantum field at a given point, the VEV of the EMT acts as the 
source of gravity in the Einstein equations, and therefore plays an important role in modelling 
a self-consistent dynamics involving the gravitational field |75j . 

The VEVs of the field square and the EMT can be obtained from a two-point function 
taking the coincidence limit in combination with the renormalization procedure. As a two-point 
function we will take the positive frequency Wightman function defined as the VEV 

W{x,x') = {0\ip{x)ip{x')\0). (9.7) 

Our choice of this function is related to that it also determines the response of the Unruh- 
DeWitt type particle detector at a given state of motion [75]. Substituting into formula (j9.7p 
the expansion of the field operator (j9.6p and using the commutation relation for the annihilation 
and creation operators, the following mode-sum formula is obtained 

W{x,x') = Y,Mx)^:ix')- (9.8) 

We can evaluate the VEVs of the field square and the energy-momentum tensor as the following 
coincidence limits: 

{0\ip'^ix)\0) = lim Wix,x'), (9.9) 

x'—*x 

{0\Tikixm = limdidiWix,x') 

+ 9ikViV - CV.Vfc - CRii\ {{)\^''{x)\Q), (9.10) 
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where di is the derivative with respect to x' and d'^ is the derivative with respect to x"^. Sim- 
ilar formulae can be written for fermionic and electromagnetic fields. The expressions on the 
right of formulae (|9.9|) . (|9.1U|) are formal as they diverge and to obtain a finite physical re- 
sult some renormalization procedure is needed. The consideration of the factors in the field 
bilinear products at different spacetime points corresponds to the point-splitting regularization 
procedure. An alternative way is to introduce in the mode-sums of the field square and the 
energy-momentum tensor a cutoff function, which makes them convergent, and to remove this 
function after the renormalization. Other regularization methods for the VEVs of the field 
square and the EMT are based on the use of the Green function (for the Green function method 
in the Casimir effect calculations see [7] and references therein) and the local zeta function [79] 
(see also [80]). The discussion for the relations between various regularization techniques can 
be found in [21 H [3 [75l |80] . 

The important thing is that, for points away from boundaries the divergences in the VEVs 
of the field square and the EMT are the same as those on background of manifolds without 
boundaries and, hence, for these points the corresponding renormalization procedure is also 
the same. However, the local VEVs diverge for the points on the boundary. These surface 
divergences are well known in quantum field theory with boundaries and are investigated for 
various types of bulk and boundary geometries (see, for example, [75l [771 EH ESI (Ml [84]). 
A powerful tool for studying one-loop divergences in the VEVs is the heat-kernel expansion 
|801 [85l I86j . In general, the physical quantities in problems with boundary conditions can be 
classified into two main groups. The first group includes quantities which do not contain surface 
divergences. For these quantities the renormalization procedure is the same as in quantum 
field theory without boundaries and they can be evaluated by boundary condition calculations. 
The contribution of the higher modes into the boundary induced effects in these quantities is 
suppressed by parameters already present in the idealized model. Examples of such quantities 
are the energy density and the vacuum stresses at the points away from the boundary and the 
interaction forces between disjoint bodies. For quantities from the second group, such as the 
energy density on the boundary and the total vacuum energy, the contribution of the arbitrary 
higher modes is dominant and they contain divergences which cannot be eliminated by the 
standard renormalization procedure of quantum field theory without boundaries. Of course, 
the model where the physical interaction is replaced by the imposition of boundary conditions 
on the field for all modes is an idealization. The appearance of divergences in the process of 
the evaluation of physical quantities of the second type indicate that more realistic physical 
model should be employed for their evaluation. In literature on the Casimir effect different 
field-theoretical approaches have been discussed to extract the finite parts from the diverging 
quantities. However, in the physical interpretation of these results it should be taken into 
account that these terms are only a part of the full expression of the physical quantity and the 
terms which are divergent in the idealized model can be physically essential and their evaluation 
needs a more realistic model. It seems plausible that such effects as surface roughness, or the 
microstructure of the boundary on small scales can introduce a physical cutoff needed to produce 
finite values for surface quantities. 

Below we are interested in the VEVs of the field square and the EMT at points away from 
boundaries. These VEVs belong to the first group of quantities and, hence, they are well-defined 
within the framework of standard renormalization procedure of quantum field theory without 
boundaries. We expect that the same results will be obtained in the model where instead of 
externally imposed boundary condition the fluctuating field is coupled to a smooth background 
potential that implements the boundary condition in a certain limit [9j. For a scalar field we 
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will assume that on the boundary the field satisfies Robin boundary condition 

(^A + Bn'Vi^(p{x) = 0, xedM, (9.11) 

where A and B are constants, and n* is the unit inward-pointing normal to the boundary. Of 
course, all results will depend only on the ratio of the coefficients in (19. lip . However, to keep 
the transition to Dirichlet and Neumann cases transparent we will write the boundary condition 
in the form (j9.1ip . Robin type conditions are an extension of Dirichlet and Neumann boundary 
conditions and appear in a variety of situations, including the considerations of vacuum effects for 
a confined charged scalar field in external fields, spinor and gauge field theories, quantum gravity, 
supergravity and braneworld scenarios [Ml EH ESI [89]. In some geometries, Robin boundary 
conditions may be useful for depicting the finite penetration of the field into the boundary with 
the 'skin-depth' parameter related to the Robin coefficient [901 191j . In examples with fermionic 
fields we will consider the bag boundary condition and in the case of the electromagnetic field 
we assume perfect conductor boundary conditions. 



10 Examples for the application of the Abel-Plana formula 
10.1 Casimir effect in x 

First we consider an example of the application of the standard Abel-Plana summation formula 
for the evaluation of the vacuum characteristics in the flat spacetime assuming that one of the 
spatial coordinates, the coordinate = y, is compactified to a circle with the length a and 
the scalar field satisfies periodic boundary condition ip{y + a) = (p{y) (for the Casimir effect in 
topologically non-trivial spaces and its role in cosmological models see [HI [U [71 [92] and references 
therein). The corresponding normalized eigenfunctions have the form 

gikny+ik-x-iuj„t 27]- 

[2a;„a(27r)-"~i]-^/^ a 

and (jj"^ = k"^ + k'^ + m? . The vector x = (x^, x^, . . . , x^~^) specifies the uncompactified spatial 
dimensions and k = |k|. From mode-sum formula (|9.8p for the Wightman function one finds 



Wpix,x',a) = - / dk-—^}2 cos(/c„Ay), (10.2) 



where Ax = x — x'. At = t — t', Ay = y — y' ■ We apply to the series over n the APF in the form 
(j2.16p . The corresponding conditions are satisfied if |Ay| + |At| < 2a. In particular, this is the 
case in the coincidence limit assuming that \y\ < a. It is easily seen that the contribution coming 
from the first integral on the right of the APF is the Wightman function for the Minkowski 
spacetime, Wm{x,x'). As a result, the Wightman function is written in the form 

Wp{x,x',a) = WM{x,x') + 2 dk—-^ dx J \^ (10.3) 

with the notation km = Vk"^ + m?. In the second term on the right of this formula we first 
integrate over the angular part of the vector k. Further, we introduce polar coordinates in the 
plane (A;, x) and integrate over the polar angle. The corresponding integrals can be found in 
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. In this way we obtain the formula 



4 



Wp{x,x' ,a) = Wm{x,x')-\ / dx ■ 

X cosh(xAy)j£,/2_i(Az\/x2 - m^), (10.4) 

where Az = [jAxp - (At)2]V2. By using the expansion (e"^ - 1)"^ = J2n=i e"""^, in the 
separate terms in the series are exphcitly integrated and we find an equivalent representation 

Wp(x,j;,o) = - — 2^ (jj_iw2 . (10-5) 

(27r) 2 

with Un = m[|Ax|2 + (Ay + an)2 — (At)^]^/^^ The n = term in this formula is the Minkowskian 
Wightman function. Formula (jl0.5p presents the Wightman function as an image sum of the 
corresponding Minkowskian functions. 

The renormalized VEV of the field square is obtained subtracting from the Wightman func- 
tion the Minkowskian part and taking the coincidence limit. As we have already separated the 
latter, this procedure here is trivial and one obtains 

/ 2v(P) 21-^ (:,2_^2^D/2-l 

i^nel = n/9T..../n^ / dx 



7r^/2r(Z)/2) Jm - 1 

(2^)(D+i)/2 (^a^){i3-i)/2 • ^ 
In the similar way, for the components of the vacuum EMT we find (no summation over i) 

(27r)(-'^+i)/2 ^ (man 

)(D+l)/2 ' ^^"-'^ 

for i = 0, 1, . . . , D — 1, and the component (T^)ron is obtained from the trace relation (j9.3p : 

(r/)ren = Tn"^ i^"^) ren ■ We have considered periodic boundary conditions. For the field with 
antiperiodic conditions (twisted scalar field) we have ip{y + a) = —ip{y). Now, in the expressions 
of the eigenfunctions kn = 27r(n + l/2)/o. The corresponding VEVs are evaluated in the way 
similar to that for the periodic case with the application of the APF in the form (j2.19p . Already 
on these simple examples we have seen two important features of the application of the APF. 
First, this enables to extract from the VEVs the part corresponding to the geometry with 
decompactified dimensions and second, to present the parts induced by the non-trivial topology 
in terms of rapidly convergent integrals. 



10.2 Casimir effect for parallel plates with Dirichlet and Neumann boundary 
conditions 

In considered examples, due to the homogeneity of the background spacetime the VEVs of the 
field square and the EMT do not depend on coordinates. This is not the case in the Minkowski 
spacetime with boundaries. The presence of boundaries breaks the homogeneity and leads to 
coordinate-dependent VEVs. Simplest examples of this kind are scalar fields with Dirichlet and 
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Neumann boundary conditions on two parallel plates with distance a. Assuming that the plates 
are located at y = and y = a, the corresponding eigenfunctions are obtained from (jlO.ip by 
the replacement e^^"-^ — > \/2sin(A;„y) with = vrn/a, n = 1,2, . . ., for Dirichlet case, and by 
the replacement e*'^"^ \/2 — 6no cos(/c„?/) with A;„ = im/a, n = 0, 1, 2, . . ., for Neumann one. 
Now, applying to the series over n in the mode-sum for the Wightman function Wj{x,x') (with 
J = D, N for Dirichlet and Neumann scalars) the APF (j2.16p . we can see that the term with the 
first integral on the right of the APF corresponds to the Wightman function wf'\x,x') for the 
geometry of a single plate at y = 0, and the term with the second integral is the part induced 
by the presence of the second plate at y = a. It can be seen that the following relations take 
place 

Wi''\x,x') = Wm{x,x')+6jWm{x,x[), (10.8) 
Wj{x,x') = Wp{x,x',2a) + 6jWp{x,x[,2a), (10.9) 

with (5d = —1, 5n = 1> and x[ is the image for the point x' with respect to the plate at y = 0: 
x[ = {t,^',-y'). 

For the geometry of a single plate at y = 0, by using the Wightman function (jlO.Sp one finds 

2\(0) _ Sjm^-^ K(g-i)/2(2m|y|) 
(2^)(^+i)/2 (2m|y|)(^-i)/2 ' 

for the VEV of the field square, and (no summation over i) 



VV' /J,ren - /n_\m4-lV9 ^o„, I , .1 N m-^ /9 ' l^iU.iUj 



, (0) _ 4<5jm^+HC - Cd) 

\ «/J,ren (27r)(^+l)/2 
-D\(0) 



K{D+i)/2{2m\y\) ^ -g(D+3)/2(2m|y|) 
(2m|y|)(^+i)/2 (2m|2/|)(«-i)/2 



+ ^{*'M°L. (10.11) 



i = 0, 1, . . . , L> — 1, {Tf^) j = 0, for the VEV of the EMT. In particular, for a conformally 
coupled massless scalar the latter vanishes. For points on the plate the VEVs diverge with the 
leading terms 

,(0) WP-l)/2) .(0) ^_JID5^{C^Cd^ (D + l\ 

W /J,rcn (4^)(D+l)/2|y|D-l' /J,i-en 2^7r(^+l)/2|y|^+l V 2 J' ^ ^ ^ 

for i = 0,1, . . . , D — 1. Note that the expressions on the right of these formulae are the VEVs 
of the field square and the EMT for a massless scalar. 

For the geometry of two plates at y = and y = a the VEVs are presented in the form 

(^^>J.en = + ^^15751^ I dx ^ ^,^^1^ [1 + 6, cosh(2xy)] , (10.13) 

for the field square, and (no summation over i) 

(22 ^ 

X I ^ ~^ + 5j [4(C - Cd)x^ - m^/D] cosh(2xy) | , (10.14) 

(TE) J ren = J, , , , / dxx^ ^ , (10.15) 

i = 0, 1, . . . , D — 1, for the EMT. Note that by using the formula 

,2a. _i cosh(2xy) = (-l)f _^2^r 

X (D-l)/2 \z=2m\an+y\^ (10.16) 

n=— oo 
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where the prime means that the n = term should be omitted, we can present the corresponding 
results in terms of series containing the McDonald function. The vacuum forces acting on the 
plates are determined by the component (r^)j].en and they are the same for Dirichlet and 
Neumann scalars. Formula (jlO.lSp for the vacuum forces acting on the plates can also be 
obtained by differentiating the corresponding vacuum energy from [93] (see also [?])• The forces 
are attractive for all values of the interplate distance. In particular, in the case of a massless 
field for the modulus of the force acting per unit surface of the plate one finds 

^.^ ,.°?°i":i. rm,j^D.N. (10.17) 



(47r)(^+i)/2a^+i 

where Cr(2^) is the Riemann zeta function. Local analysis of a quantum scalar field confined 
within rectangular cavities is presented in |94j . 

In the dimension D = 3 with perfect conductor boundary condition on two parallel plates, 
the electromagnetic field is presented as a set of two types of modes corresponding to massless 
Dirichlet and Neumann scalars. The corresponding VEV for the EMT does not depend on 
coordinates and in the region between the plates is given by [95j 

{T.hEl,ren = " 7^diag(l, 1, 1, -3). (10.18) 

Due to the cancellation between Dirichlet and Neumann modes, the electromagnetic EMT is 
uniform in the region between the plates and vanishes in the outside regions. 



11 Casimir effect for parallel plates with Robin boundary con- 
ditions 

As we have demonstrated in the previous section the application of the APF provides an efficient 
way for the investigation of the boundary-induced effects. Here we will see that already in the 
case of two parallel plate geometry with Robin boundary conditions this procedure needs a 
generalization [32]. This can be done by using the GAPF. As before, we assume that the plates 
are located at y = ai = and y = a2 = a, and the field obeys boundary conditions ()9.1ip with 
coefficients Aj, Bj for the plate at y = aj. In the region between the plates the eigenfunctions 
are presented in the form 

ip^ = l3{ky)e'''-^-'^^''y^' cos [kyy + a{ky)) , (11.1) 



where oj{ky) = ^/k^ + k"^ + m?, and a{ky) is defined by the relation 

\ [5j = {-iy-^B,/Aj. (11.2) 



,2ia{ky) ^ iPlky - 1 

iPiky + 1'' 



The corresponding eigenvalues for ky are obtained from the boundary conditions and are solu- 
tions of the following transcendental equation: 

F{z) = (l — 6162-^^) sin z — {bi + 62) z cos z = 0, z = kya, bj = f3j/a. (11-3) 

The expression for the coefficient P{ky) in (jll.ip is obtained from the normalization condition: 



1 H , ^ ' cosikyU + 2aiky)) 

kya 



(11.4) 
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The eigenvalue equation (jll.3p has an infinite set of real zeros which we will denote by ky = Xn/a, 
n = 1, 2, . . .. In addition, depending on the values of the coefficients in the boundary conditions, 
this equation has two or four complex conjugate purely imaginary zeros ±zy/, > (see |32|). 

Substituting eigenfunctions (jll.ip into mode-sum formula (|9.8p . for the positive- frequency 
Wightman function in the region between two plates one finds 



W{x,x') 



dke 



ik-Ax 



^ /?2(A;j^)e-^^(^^)^* 

ky=Xn/a,iyi/a 

X cos{kyy + a{ky)) cos{kyy' + a{ky)). (H-S) 
To obtain the summation formula over the zeros A^, in the GAPF as a function g{z) we choose 



9[z) 



[(l — 61622;^) cos z + {bi + 62) z sinz] 



F{zy 



(11.6) 



Substituting (jll.6p into (j2.1ip and taking the limit a — > 0, under the assumption that the 
poles Hyi are excluded by small semicircles on the right half-plane, one obtains the following 
summation formula 132 1 



E r 

z=X„,iyi 



vr/(z) 



-|- cos(2; -|- 2a) sin z/z 



+ 



^/(0)/2 
61 62 - 1 



dt- 



dzfiz) 

I - /(te— /2) 



Ojbl) 
2bi 



{bit-l){b2t-l) Of. _ 1 
(f)lt+l)(b2t+l) 

h{e'^ /hi) + h{cie~'^/bi) 



(11.7) 



where h{z) = {h\z'^ + l) f{z). For the case 62 = 0, 61 < 0, with an analytic function f(z) this 
formula is derived in |91j . 

For the summation over the eigenmodes A„ in (jll.5p as a function f{z) in (jll.7p we take 

^—iLu{z/a)At 

f(z) = — cosizyla + aiz/a)) cosizy' la + aizla))^ (H-S) 

aujyz/a) 

with first-order poles at z = zizi/bj. By making use of the definition for a{k) we see that 
g2«a(o) _ _-|^^ hence cos(2a(0)) = —1, which implies that /(O) = 0. The resulting Wightman 
function from (jll.5p is found to be 

4 



W{x,x') = Wo{x,x') + 



jk-Ax 



du 



{2tt)D 

cosh(yn -I- d(n)) cosh(y''u -|- a{u)) ^^^^ 



h2 



{I3xu-l)(l32u~l) 2au 
(/3it+l)(/32t+l) 



1 



11. 



where km = Vk"^ + m? and the function dij (t) is defined by the relation 

2aiu) ^ PlU - 1 



PlU+1 



(11.10) 



In formula pT9]) . 

Wo{x,x') = Wm{x,x') + 



dk 



(2vr) 



D 



00 ^—iuj{u)At 
du- 



uj{u) 



COS [u{y + y') + 2a{u)] 



+ 



e{pi)e 

(27r)^-i/3 



dk 



exp(«k • Ax - iAt^/kl - 



(11.11) 
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is the Wightman function for a single plate located at y = and WMix, x') is the Wightman 
function in the Minkowski spacetime without boundaries. The last term on the right comes 
from the bound state present in the case /5i > 0. To escape the instability of the vacuum state, 
we will assume that mf3i ^ 1. On taking the coincidence limit, for the VEV of the field square 
we obtain the formula 

(^^)R.en = (^^)gL + ^VtW / W.ll)W.t-l) - ^^'^"(ty + «(^))' (11-12) 

^ ' ^ ™ (/3ii+l){/32t+l)^ ^ 



where 



is the corresponding VEV in the region y > for a single plate at y = 0. The surface divergences 
on the plate at y = are contained in this term. The second term on the right of formula (jll.l2p 
is finite at y = and is induced by the second plate located at y = a. This term diverges at 
y = a. The corresponding divergence is the same as that for the geometry of a single plate 
located at y = a. Note that in obtaining (I11.13P from (lll.lip we have written the cos function 
in the second integral term on the right of (jll.lip as a sum of exponentials and have rotated 
the integration contour by the angle 7r/2 and by — 7r/2 for separate exponentials. For /3i > 
the corresponding integrals have poles ±i//3i on the imaginary axis and the contribution from 
these poles cancel the part coming from the last term on the right of (lll.lip . 

The VEV of the EMT is evaluated by formula (I9.10p . By taking into account formulae 
(jll.Op . (jll.l2p . for the region between the plates one finds (no summation over I, see j32j for the 
case of a massless field) 

{TbK,re. = {Tl)Z.. + 26i^j^ dt ,,U,,,J,, -m,y), (11.14) 

where 

fi{t, x) = [At^ (Cd - C) + m'^/D] cosh {2ty + 2a{t)) - [t^ - m^)/D, (11.15) 
for / = 0, 1, . . . , L> - 1, and //)(*, y) = In formula (fTITi]) . 

x|^[4(CD-C)t^ + mVl?], (11.16) 

for i = 0, 1, . . . , D — 1, and {T£)'^^^^ = 0, are the VEVs in the region y > induced by a single 
plate at y = 0, and the second term on the right is the part of the energy-momentum tensor 
induced by the presence of the second plate. In [50j , formulae (I11.12P and (I11.14P are obtained 
as limiting cases of the corresponding results for the geometry of two parallel plates uniformly 
accelerated through the Fulling-Rindler vacuum (see below). For a conformally coupled massless 
scalar field the vacuum EMT is uniform and traceless. Note that we have investigated the 
vacuum densities in the bulk. For Robin boundary conditions in addition to this part there is a 
contribution to the EMT located on the plates. 

Vacuum forces acting on the plates are determined by (r^)R^ren- This component is uniform 
and, hence, is finite on the plates. The latter property is a consequence of the high symmetry 
of the problem and is not valid for curved boundaries. In dependence of the values for the 
coefficients (dj the vacuum forces can be both attractive or repulsive. In particular, the vacuum 
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forces are repulsive in the case of Dirichlet boundary condition on one plate and Neumann 
boundary condition on the other. In the case of a massless scalar field for the modulus of the 
corresponding force from (|11.14|) with i = I = D one finds = (1 - 2-^)Fd, where is 
given by formula (jlO.lTp . Taking D = 3 and summing over the polarization states, from here 
we obtain the result of [96j for the electromagnetic Casimir force between two parallel plates 
one of which is a perfect conductor and the other is infinitely permeable. 

Series over the zeros of the function F{z) from (I11.3P appear also related to the Casimir 
effect for two parallel plates with non-local boundary conditions 

nl^^dM^^^) + I dxji /j(|x|| - ^\\\Mx'n =0,y = aj, (11.17) 

where n'^j-^ is the inward-pointing unit normal to the boundary at y = aj. In the region between 
the plates the corresponding eigenfunctions have the form (jll.ip and the eigenvalues for k 
are zeros of the function F{ka) in ()11.3p . where now the coefficients are given by the formula 
l/bj = (-lp-iaFj(A;||) with 

f C2-jr-\{D-l)/2 foo 

F,(A;||) ^ / dx|| /.(IxijDe^i'Fll = ^ (^,3^^, / dnu'^^-'y^ fj{n)J^n-3)/2{nk\^). (11.18) 
The corresponding VEVs are investigated in |33j . 

12 Scalar vacuum densities for spherical boundaries in the global 
monopole bulk 

Historically, the investigation of the Casimir effect for a spherical shell was motivated by the 
Casimir semiclassical model of an electron. In this model Casimir suggested that Poincare 
stress, to stabilize the charged particle, could arise from vacuum quantum fluctuations and the 
fine structure constant can be determined by a balance between the Casimir force (assumed 
attractive) and the Coulomb repulsion. However, as it has been shown by Boyer [97], the elec- 
tromagnetic Casimir energy for the perfectly conducting sphere is positive, implying a repulsive 
force. This result has later been reconsidered by a number of authors [98^ [81] I99j . More re- 
cently new methods have been developed for this problem including direct mode summation 
techniques jlOOj llOll I102j (see also O El [7| and references therein). However, the main part 
of the studies was focused on global quantities such as the total energy and the vacuum forces 
acting on the sphere. The investigation of the energy distribution inside a perfectly reflecting 
spherical shell was made in |103j in the case of QED and in [104] for QCD. The distribution 
of the other components for the electromagnetic EMT inside as well as outside the shell can 
be obtained from the results of \105\ I106j . In these papers consideration was carried out in 
terms of Schwinger's source theory. In |34l \35\ [36] the calculations of the the VEVs for the 
EMT components inside and outside the perfectly conducting spherical shell are based on the 
GAPE. The same method was used in [38j for the evaluation of the Wightman function and the 
VEVs for the EMT of a massive scalar fleld with general curvature coupling and obeying Robin 
boundary condition on spherically symmetric boundaries in D-dimensional Minkowski space. In 
[107J the local VEVs outside a spherical boundary and the corresponding energy conditions are 
investigated using the calculational framework developed in [108j. More complicated problems 
with spherical boundaries in the background spacetime of a global monopole are considered in 
[51] [52[ [53] [Ml [55] for both scalar and fermionic flelds. As the corresponding results in the 
Minkowski bulk are obtained as special cases, here we discuss the application of the GAPE for 
the geometry of spherical boundaries in the global monopole bulk. 
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Global monopoles are spherically symmetric topological defects created due to phase tran- 
sition when a global symmetry is spontaneously broken and they have important role in cos- 
mology and astrophysics. The global monopole was first introduced by Sokolov and Starobin- 
sky |109j . A few years later, the gravitational effects of the global monopole were considered 
in |110j . where a solution is presented which describes a global monopole at large radial dis- 
tances. Neglecting the small size of the monopole core, in the hyper spherical polar coordinates 
{r,'&,4>) = {r,9i,02, . . . On,(p), n = D — 2, the corresponding metric tensor generalized to an 
arbitrary number of spatial dimensions can be approximately given by the line element 

ds^ = - dr^ - a^r^dnl,, (12.1) 

where d^lj^ is the line element on a surface of the unit sphere in Z)-dimensional Euclidean space, 
the parameter a is related to the symmetry breaking energy scale in the theory. The solid angle 
corresponding to (|12.ip is a Sd with Sd = 2tt^/'^/T{D/2) being the total area of the surface of 
the unit sphere in D-dimensional Euclidean space. The nonzero components of the Ricci tensor 
and Ricci scalar for the metric corresponding to line element ()12.ip are given by expressions (no 
summation over i) 

\-o? 1 - a2 

K = n^r^^ i? = n(n+l)^^, (12.2) 

where the indices i = 2, 3, . . . , D correspond to the coordinates ^i, 6*2, . . . , respectively. Note 
that for a 7^ 1 the geometry is singular at the origin (point-like monopole, for vacuum polariza- 
tion effects by a global monopole with finite core see |lllj ). In this section we will consider the 
applications of the GAPF for the evaluation of the Wightman function, VEVs of the field square 
and the EMT for a scalar field with general curvature coupling parameter under the presence of 
spherical boundaries concentric with the global monopole. 

12.1 Region inside a single sphere 

For the region inside a sphere with radius a the complete set of solutions to (|9.ip is specified by 
the set of quantum numbers a = (A, m^), where m/c = (mg = m-i, . . . , m„) and mi, m2, . . . , m„ 
are integers such that ^ m„_i ^ 'm„_2 ^ ••• ^ ra\ ^ — m„_i ^ m„ ^ Tan-\- The 
corresponding eigenfunctions have the form 



^.{x) = /3,r-"/2 J,, {Xr)Y{mk; ^, 0)6-'"*, uj = V^^ + m^, / = 0, 1, 2, . . . , (12.3) 
where the order of the Bessel function is given by the formula 



1 

1^1 = - 

a 



1 + ^) +{l-a'^)n({n + l)C-^ 



nl/2 



(12.4) 



and Y{mk; (p) is the hyper spherical harmonic of degree The coefficient f3a is found from the 
normalization condition and is equal to 



^9 ATy, (Aa) 

= ... ? D-i ^ 12.5 



with Tiy{z) defined by (jS.Sp . From boundary condition (j9.1ip for eigenfunctions (jl2.3p . we see 
that the eigenvalues for Xa are zeros of the function J^^{Xa) with the barred notation from ()3.ip 
and with the coefficients 

A = A-nB/2, B = n^B/a, (12.6) 
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where = — 1 for the region inside the sphere. In the notations of Section 3, for the eigenvalues 
of A and m one finds 

A = A,,,fc/a, LJ = ^\lja^ + m?. (12.7) 

Substituting the functions (jl2.3p into (j9.8p and using the addition formula for the hyperspherical 
harmonics [TSl, for the Wightman function we obtain 



1=0 



W{x,x') = ^_^ Y.(^l + n)cf\cos9) 



E 7' "'^ J..(A.,,r/a)J.,(A,,,r7a)e"--^^^^ (12.8) 

k=i yA^,,fc + "^^«^ 

where At = t — t' . In this formula, Cp{x) is the Gegenbauer polynomial of degree p and order 
q, and 9 is the angle between the directions (■!?,(?!>) and (??',0'). To sum over we will use 
summation formula (jS.lSp . taking 



f{z) = zJ^^ {zr/a)J^^ (zr7a)e-^^V^Va^W_ (^2.9) 

For \z\/a < m this function satisfies the relation /(ze"'^*/^) = —e~'^'^'^^f{ze'^^^'^). As a result, the 
integral on the right of formula (j3.18p over the interval (0, ma) vanishes. Conditions (j3.4p are 
satisfied if r + r' + |At| < 2a. After the application of formula ()3.18p the Wightman function is 
presented in the form |51| 



(x, x') = Wrnix, x') + {^{xMx'))^, (12.10) 

where 

W^m(x,x') = f^J]/7^Cr/'(cose) / dz ' J,^{zr)J,^{zr'), (12.11) 

and 



2l + n n/2 K^iiaz) 



X . 

\J z^ — rn?' 



cosh(At\/z2 _ m2). (12.12) 



The contribution of the term (jl2.1ip to the VEV does not depend on the sphere radius, whereas 
the contribution of the term (|12.12p vanishes in the limit a ^ oo. It follows from here that 
expression (jl2.1ip is the Wightman function in the unbounded global monopole space (for 
the corresponding Euclidean Green function see [112^ 1113] ). This can also be seen by explicit 
evaluation of the mode-sum using the eigenfunctions for the global monopole spacetime without 
boundaries. As we see, the application of the GAPF allows to extract from the bilinear field 
product the contribution due to the unbounded monopole spacetime, and the term (jl2.12p can 
be interpreted as the part of the VEV induced by the spherical boundary. 

The VEV of the field square is obtained by using formula ()9.9p . For < r < a the divergences 
are contained in the first summand on the right of Eq. (112. lOp only. Hence, the renormalization 
procedure for the local characteristics of the vacuum, such as field square and EMT, is the 
same as for the global monopole geometry without boundaries. This procedure is discussed in 
a number of papers |112l 1113^ I114j and is a useful illustration for the general renormalization 
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prescription on curved backgrounds. Here we are interested in the parts of the VEVs induced 
by the presence of a spherical shell. Using the relation C"'^^(l) = r(/ + n)/r(n)/!, for the 
corresponding boundary part in the VEV of the field square we get 



a 



^ Iui{az) 



where 



Di = {2l + D- 2) 



T{l + D- 2) 



(12.13) 



;i2.14) 



is the degeneracy of each angular mode with given As it has been noted before, expression 
(|12.13p is finite for all values < r < a. For a given I and large z the integrand behaves as 
^2z{r/a-i) j ^ and, hence, the integral converges when r < a. For large values I, introducing a new 
integration variable y = z/vi in the integral of Eq. (|12.13p and using the uniform asymptotic 
expansions for the modified Bessel functions [65], it can be seen that the both integral and sum 
are convergent for r < a and diverge at r = a. For points near the sphere the leading term of 
the corresponding asymptotic expansion over the distance from the boundary has the form 

(i-2<5Bo)rr^-i^ 



(47r)(^+i)/2(a-r; 



D-l ' 



(12.15) 



This term does not depend on the mass and on the parameter a and is the same as that for a 
sphere on the Minkowski bulk. As the purely gravitational part is finite for r = a, we see that 
near the sphere surface the VEV of the field square is dominated by the boundary-induced part. 
The boundary- induced VEV (jl2.13p is zero at the sphere center for (a"^ — 1)^ > 0, non-zero 
constant for (a"^ — !)( = 0, and is infinite for (a^^ — !)( < 0. 

Substituting the Wightman function (|12.10p into Eq. ([QJOp . for the VEV of the EMT inside 
the spherical shell one finds [S^ 



^12.16) 



where {T^)^ is the corresponding VEV for the monopole geometry when the boundary is absent, 
and the part {T^)h is induced by the presence of the spherical shell. From Eqs. (jl2.1ip . ()12.12p . 
and (jl2.13p , for the boundary- induced parts of the EMT components one obtains (no summation 
over i) 



a 



^^^^,K,,iza)Fy/ [L^izr)] 



where for a given function f{y) we have introduced the notations 



r < a, 



;i2.17) 



f!^'> [f{y)] 

f!^'> [fiy)] 
Fif [f{y)] 



(1-4C) /'(y) - ^/(y)/'(2/) 



(y), 



(12.18) 
(12.19) 



{4c-i)f'Hy)--f{y)f'{y) 
y 

(n + Ijy^ 



f\y), 



;i2.20) 
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with = 4(n + 1)^ — n and in (|12.20p i = 2,3, . . . , D. It can be seen that components (|12.17p 
satisfy the continuity equation ^k{T^)h = and are finite for < r < a. They are zero at the 
sphere center for nC(a~^ — 1) > 1, are non-zero constants for n^(a~^ — 1) = 1, and are infinite 
otherwise. These singularities appear because the geometrical characteristics of global monopole 
spacetime are divergent at the origin. However, note that the corresponding contribution to the 
total energy of the vacuum inside a sphere coming from {TQ)h is finite due to the factor r^~^ 
in the volume element. 

Expectation values (jl2.17p diverge at the sphere surface, r ^ a. The corresponding asymp- 
totic behavior can be found using the uniform asymptotic expansions for the modified Bessel 
functions, and the leading terms are determined by the relations 

,^0^ /r^^ D<^{Ti\ Dr{{D + l)/2){C-CD) ... 

As for the field square, these terms do not depend on mass and parameter a, and are the same 
as for a spherical shell in the Minkowski bulk. In general, asymptotic series can be developed 
in powers of the distance from the boundary. The corresponding subleading coefficients will 
depend on the mass, Robin coefficient, and parameter a. Due to the surface divergencies near 
the sphere surface the total vacuum EMT is dominated by the boundary- induced part (T^^)b- 

12.2 Vacuum densities in the region between two spheres 

In this subsection we consider the VEVs of the field square and the EMT on background of 
the geometry described by ()12.ip . assuming that the field satisfies Robin boundary condition 
(|9.1ip on two concentric spherical boundaries with radii a and b, a < b. We will consider the 
general case when the coefficients in the boundary conditions for the inner and outer spheres are 
different and will denote them by Aj and Bj with j = a, b. For the region between two spheres 
the corresponding eigenfunctions can be obtained from formula ()12.3p with the replacement 

J,,(Ar) ^ gu,iXa,Xr) = J,^{Xr)Yj^\Xa) - J^f iXa)Y,,{Xr), (12.22) 

where the functions with overbars are defined by (j4.2p with 

Aj = Aj - Bjn/2, Bj = rijBj/j, j = a, b, Ua = 1, rib = -1- (12.23) 

Substituting eigenfunctions (jl2.3p into the normalization condition and using the value for the 
standard integral involving the square of a cylinder function [66], one finds 



^ ■K'^\T";\b/a,\a) 



where T^^{'q, z) is defined by (j4.7p . 

The functions chosen in the form (jl2.22p satisfy the boundary condition on the inner sphere. 
From the boundary condition at r = 6 one obtains that the corresponding eigenmodes are 
solutions to the equation 



~iab 



Cf(Va,Aa)=0, (12.25) 



with the function {t],z) defined by formula ()4.1|) . Hence, in the region between two spheres 
the eigenvalues for A and uj are 



A = 7.„fc/a, = ^lljci^ + m^- (12.26) 
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Substituting the eigenfunctions into mode-sum (|9.8|) and using the addition formula for the 
spherical harmonics, for the Wightman function one finds 



with the function 



2/ l\—n/2 ^ 

= J2(^l + n)C;^\cose)^hij,^^,)T-;>{b/a,^,,,k), (12.27) 

1=0 k=l 



h{z) = -g^^{z,zr/a)g^^{z,zr'/a). (12.28) 

V + m^a^ 

To sum over k we will use summation formula (j4.14p . The corresponding conditions for function 
(|12.28p are satisfied if r + r' + | Af| < 26. In particular, this is the case in the coincidence limit 
for the region under consideration. For the Wightman function one obtains 



h{z)dz 2 f°° ^ zQau{az,bz) 



xG^,"-* (az, zr )G^y^ iflz, zr') cos h{AtVz^ -m^)}, (12.29) 

where we have introduced the notation 

Gl^\z,y) = Uy)Ki^\z) - fJ\z)K,{y), j = a,b, (12.30) 

and the function ^}au{cLz,bz) is defined by relation (j4.15p . In the limit 6 — > oo the second 
integral on the right of (jl2.29p tends to zero, whereas the first one does not depend on b. It 
follows from here that the term with the first integral in the figure braces corresponds to the 
Wightman function for the region outside a single sphere with radius a on background of the 
global monopole geometry. The latter we will denote by W^°'\x, x'). So, the Wightman function 
in the region between two spheres is presented in the form [53] 

i\ jxAa)/ f\ a^~^ 2l + n n/2, f°° , z^aui{az,bz) 



X G^^^ (az, rz)G^^I^ (az, r'z) cosh(At\/z2 _ m^). (12.31) 

To extract from the function W^"'\x,x') the part induced by the presence of the sphere we use 
the identity 

iyfliv'l = My)My') -\jz^H[^'\y)H^'\y')- (12.32) 
Ji{z) + yu{z) ^f^^Hr{z) 

The contribution from the first term on the right of this relation with y = xr/a and y' = xr' /a 
gives the Wightman function VFm(x, x') for the geometry without boundaries. In the part coming 
from the second term we rotate the contour for the integration over z by the angle 7r/2 for s = 1 
and by the angle — 7r/2 for s = 2. Introducing the modified Bessel functions, one finds 

W(^\x,x') = W^ix,x')-^j:^Cr%ose)j^dzz 

^I^Ed0£d^,osHAtyG^^^). (12.33) 
Kjj^{az) -Jz^-ni^ 
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Note that in the coincidence hmit, x' = x, the second summand on the right hand-side of (|12.3ip 
will give a finite result for a ^ r < b, and is divergent on the boundary r = h. It can be seen 
that for the case of two spheres the Wightman function in the intermediate region can also be 
presented in the form 



W{x,x') = Wi''\xy)-^Y.i^.Cr%ose) (12.34) 

^ 1=0 ^ ' 

^^ zVL ^^ XazM ^ip) (^^^ (^^^ ^/^) cosh(AtVz2 - m2), 

y — w?' 

with W^^{x^ x') being the Wightman function for the vacuum inside a single sphere with radius 
6, and 

nux, Xx) = -r-, 5 Sm""^ — r^- (12.35) 

In the coincidence limit, the second summand on the right of formula (jl2.34l) is finite for a < 
r ^ b and diverges on the boundary r = a. 

The VEV of the field square is obtained computing the Wightman function in the coincidence 
limit x' — > X. In the region between the spheres, from (112. 31|) and (112. 34|] we obtain two 
equivalent forms [53] 

(OI^^IO) = - ^^f; A r dz'-^^j4^Gii^\jz,rz), (12.36) 

with j = a, b, and (99^) '-•'^ is the VEV for the case of a single sphere with radius j. 

Using the Wightman function from (jl2.3ip and the VEV for the field square from (jl2.36p . 
the components of the vacuum EMT can be presented in two equivalent forms corresponding to 
j = a,b (no summation over i) [53] : 

mtlO) = {Tt)^^^ - f^EA rdz^!%I^F(;)[Gg)(,.,rz)]. (12.37) 

Here {T^Y^^i is the vacuum EMT for the case of a single sphere with radius j and the functions 
Fu'^ [/(y)] are defined by relations ()12.18p - ()12.20p . In the region outside a single sphere the 
expressions for the boundary-induced parts are obtained from the corresponding formula for the 
interior region given in the previous subsection by the replacements I ^ K. The scalar vacuum 
densities for spherical boundaries in the Minkowski spacetime (see [38]) are obtained from the 
formulae in this section as a special case taking a = 1. Note that in this case vi = I + n/2. 

The vacuum force per unit surface of the sphere at r = j is determined by the \ - component 
of the vacuum EMT at this point. For the region between two spheres the corresponding effective 
pressures can be presented as a sum of two terms: 

P^'^=P?+P%^, j = a,b. (12.38) 

The first term on the right is the pressure for a single sphere at r = j when the second sphere 
is absent. This term is divergent due to the surface divergences in the VEVs of local physical 
observables. Here we will be concerned with the second term on the right of Eq. (jl2.38p which is 
the pressure induced by the presence of the second sphere, and can be termed as an interaction 
force. Unlike to the single shell parts, this term is free from renormalization ambiguities and is 
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determined by the last term on the right of formula (|12.37p . From the formula for the EMT, for 
the interaction parts of the vacuum forces per unit surface one finds [53] 



rijO 



l-D 



^ 1 
dj 



oo 

E 

1=0 



Di 



zdz 



az) 



Rl^Xhz) 



B]iz''f + iyf 



li^\bz)Kl':>iaz) 



;i2.39) 



where C is defined after formula p2.2U|) . The expression on the right of this formula is finite for 
all non-zero distances between the shells. For Dirichlet and Neumann scalars the second term 
in the square brackets is zero. In the case of Dirichlet scalar the interaction forces are always 
attractive. For the general Robin case the interaction force can be either attractive or repulsive 
in dependence on the coefficients in the boundary conditions. 



13 Casimir effect for a fermionic field in a global monopole back- 
ground with spherical boundaries 

13.1 Vacuum energy-momentum tensor inside a spherical shell 

Our interest in this section will be the VEV of the EMT for a fermionic field with mass M 
induced by spherical boundaries in the global monopole spacetime. The dynamics of the field 
on a curved spacetime is governed by the Dirac equation 

ij^{di+Ti)i;-Mi; = , (13.1) 

where 7^ are the Dirac matrices and F/ = 7fcV/7^/4 is the spin connection with being the 
standard covariant derivative operator. We will assume that the field satisfies bag boundary 
condition 

(1 + iVnbi) V = , r = a, (13.2) 

on the sphere with radius a, is the outward-pointing normal to the boundary. Expanding 
the field operator in terms of a complete set of single-particle states ■*} and making 

use of the standard anticommutation relations, for the VEV of the EMT one finds the following 
mode-sum formula 

{0\T,k\0) = Y.Tik{i^jf\x),xl:'f\x)}. (13.3) 

/3 

For the geometry under consideration the eigenfunctions are specified by the set of quantum 
numbers /3 = {kjma), where j = 1/2, 3/2, . . . determines the value of the total angular momen- 
tum, m = —j, . . . ,j is its projection, and o" = 0, 1 corresponds to two types of eigenfunctions 
with different parities. These functions have the form 

u} = ±E, E = Vk^ + M^ , u^ = i±ll^-'^ (13.5) 

a 2 

where n = r/r, a = {a^,a'^,a^) with the curved space Pauli matrices a^. In Eq. ()13.4p . Zi^{x) 
represents a cylinder function of the order and Oj;^m(0,(^) are the standard spinor spherical 
harmonics with la = j — ?T'(t/2. As in mode-sum formula ()13.3p the negative frequency modes 
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are employed, in the discussion below we will consider the eigenfunctions p3.4p with the lower 
sign. 

For the region inside a spherical shell one has Z^{x) = Jy{x). The imposition of the boundary 
condition on the eigenfunctions ()13.4p leads to the following equations for the eigenvalues 



Ju^{ka) = , 

where for a given function F{z) we use the notation 

F{z) = zF\z) + {f,a - + f^l - nav)F{z) , a = 0, 1, 



(13.6) 



(13.7) 



with = Ma. Let us denote by A[,^ ^ = ka, s = 1, 2, . . . , the roots to equation (jl3.6p in the 
right half-plane, arranged in ascending order. By using the standard integral for the Bessel 
functions, for the normalization coefficient one finds 



with 



Tl(z) 



+ (/X - n^u){fj.a - V^^ + mI) + 



(13.^ 



:i3. 



Substituting eigenfunctions (jl3.4p into Eq. (|13.3p , the summation over the quantum number 
m can be done by using the standard summation formula for the spherical harmonics. For the 
EMT components one finds (no summation over i) 



m\o) 



E 



2j + l 



j=l/2 <j=0,l s=l 



E J2^Liz)f!^l[z,JuAzr/a)l 



where we have introduced the notations 



f!n)[z,9Ay)] 
f!^hz,9Ay)] 



z^{2u + na) 



aliy) + gl+nAy) 



2v + n^j 

y 



9y{y)9v+nAy) 



9v{y)9v+nAy) > ^ = 2,3. 



(13.10) 

(13.11) 
(13.12) 

(13.13) 



In order to obtain a summation formula for series over the zeros in the GAPF 
function g{z) we choose 

■ Yu{z) 



9[z) 



I — 



Mz) 



f{z) 



;i3.i4) 



with a function f{z) analytic in the right half-plane Rez > 0. By making use of the asymptotic 
formulae for the Bessel functions for large values of the argument, the conditions for the GAPF 
can be written in terms of the function f{z) as follows: 



|/(z)| < e(x)e''l^l , z = x + iy, \z\ 



oo 



;i3.15) 



where c < 2 and e(x) for x oo. By using the Wronskian relation for the Bessel functions, 
one can see that Y{Xl^ g) = 2/['kJ{XI^ g)]. This allows to present the residue term coming from 
the poles of the function g{z) in the form 



vri Resg{z)=Tl{Xi^,)f{Xl^, 



(13.16) 
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Substituting p3.14p and p3.16p into the GAPF ([XTT]) and taking in this formula the hmit a ^ 
(the branch points z = ii^Ua are avoided by semicircles of small radius), we obtain that for a 
function f(z) analytic in the half-plane Rez > and satisfying condition (|13.15|) the following 
formula takes place [53] 



lim 

1 

vr 



u,sJ 



.s=l 

oo 



dx f{x) 



dx 



{+)/ 



-Res/(z)^— 
2z=o J^(z) 

(-)/ 



(13.17) 



where on the left < 6 < A[,^ , j^. In formula (|13.17|) we use the notations 



{z) = zF'{z) + (Ma - V ^2e±-'' + //2 - n,z.)F(z) (13.18) 

for a given function F{z). 

We apply formula ()13.17p to the series over s in formula (jl3.10p for the VEVs of the energy 
density and vacuum stresses. As it can be seen from expressions (113.1ip - (|13.13p . the corre- 
sponding functions f{z) satisfy the relation 

e-'^"7(xe"*/2) ^ _e^'-Y(xe-"*/2) ^ ^j. q ^ x < . (13.19) 

By taking into account that for these values x one has F(+)(x) = F(-) (x), we conclude that the 
part of the integral on the right of Eq. (113. 17p over the interval {0, Ha) vanishes. 

As a result, after the application of summation formula (jl3.17p . the components of the 
vacuum EMT can be presented in the form 

{0\T,'\0) = {Thm + {Tth, (13.20) 

where {T^)-[a does not depend on the radius of the sphere a and is the contribution due to 
unbounded global monopole spacetime. The corresponding quantities for the massless case are 
investigated in ^115j . The second term on the right of formula ()13.20p is induced by the presence 
of the spherical shell and can be presented in the form (no summation over i) 



with 



Fl"[:c, /„(!,)] = (13.23) 

Here and below I = j + 1/2, v = vi = l/a + 1/2, and for given functions /(x) and g{x) we use 
the notation 

W[f{x),g{x)] = [xf{x) + {iia + y)f{x)\ 

X [x5'(x) + {na + y)g{x)] + (x2 - ^ll)f{x)g{x). (13.24) 



52 



It can be easily checked that for a massless spinor field the boundary-induced part of the vacuum 
EMT is traceless and the trace anomalies are contained only in the purely global monopole part 
without boundaries. These expressions diverge in a non-integrable manner as the boundary is 
approached. The energy density and azimuthal pressure vary, to leading order, as the inverse 
cube of the distance from the sphere. At the sphere center the boundary parts vanish for the 
global monopole spacetime {a < 1) and are finite for the Minkowski spacetime (a = 1). In the 
limit of strong gravitational field, corresponding to small values of the parameter a, describing 
the solid angle deficit, the boundary- induced part of the vacuum EMT is suppressed by the 
factor exp[— (2/a)| ln(r/a)|] and the corresponding vacuum stresses are strongly anisotropic: 
{T^^^ ~ a^Tl)^. Having the components of the energy-momentum tensor we can find the 
corresponding fermionic condensate (0|'i/'V'|0) making use of formula Tji = M'ijj^ for the trace 
of the EMT. Note that the formula for the VEV of the EMT in the exterior region (r > o) is 
obtained from (jl3.2ip by the replacement Ky in the denominator and in the argument of 
the function F^^ [M]. 

Fermionic vacuum densities induced by a spherical shell in the Minkowski bulk are obtained 
from the results of this section as a special case with a = 1. The previous investigations on the 
spinor Casimir effect for a spherical boundary (see, for instance, [31 HI [U [71 1116] and references 
therein) consider mainly global quantities, such as total vacuum energy. For the case of a 
massless spinor the density of the fermionic condensate {^\'4)tjj\Q) is investigated in [117j (see also 

|Z!)- 

13.2 Vacuum expectation values of the energy-momentum tensor between 
two spherical shells 

In this subsection we consider the region between two spherical shells concentric with the global 
monopole on which the fermionic field obeys bag boundary conditions: 

(l + iVn[")) ^l^^ = 0,w = a,b, (13.25) 

where a and b are the radii for the spheres, a < b, nf"'^ = is the outward-pointing normal 

to the boundaries. Here and below we use the notations n^^^ = — 1, n^^^ = 1. The corresponding 
eigenfunctions have the form ()13.4p . Note that in terms of the function Zy^[kr) the boundary 
conditions (|13.25p take the form 

Zy^{kw) = -Zy^+nAkw)^^^ , (13.26) 

with w = a,b. In the region between two spherical shells, the function Zy^{kr) is a linear 
combination of the Bessel and Neuamnn functions. The coefficient in this linear combination is 
determined from the boundary condition (|13.26p on the sphere r = a and one obtains 

ZyAkr)=gi''J{ka,kr) ^ Jy^ikr)Y,^^\ka) - Y,^{kr)jil\ka), (13.27) 

where for a given function F{z) we use the notation 

F(-)(z) = zF'iz) + [n(-)(^^ - ^/z^ + fil) - n^v.]F{z) , (13.28) 

with w = a,b, and = Mw. Now, from the boundary condition on the outer sphere one finds 
that the eigenvalues for k are solutions to the equation 

Cl{b/a,ka) ^ ji:Hka)Y(^){kb) - fjf {ka)j(,''J (kb) = 0. (13.29) 
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Below we denote by 7^ 



ka, s = 1,2,..., the positive roots to equation (|13.29|) . arranged 



in ascending order, Tu^ s < Tua s+i- Substituting the eigenfunctions into the normahzation 
integral and using the standard integrals for cylinder functions (see, for instance, |66|). for the 
normalization coefficient of the negative frequency modes one finds 



A 



Tl-\7j,ka), ka = -fl^,, 



where we have introduced the notation 



j{b)2 



{r]z) 



with 



w 



k^ + {M- E){M - n^m^'"^/w) + 



n 



2wE 



(13.30) 



(13.31) 



(13.32) 



and ka = z. Substituting eigenfunctions (jl3.4p with (jl3.27p into Eq. (jl3.3p . the summation over 
the quantum number m can be done by using the standard summation formula for the spherical 
harmonics. For the EMT components one finds (no summation over i) 



2j + l 



, 32a2a3r 

3=1/2 o-=0,l s=l 



^ Y,Tl:\v:z)f^l[z,g^^}{z,zr/a)] 



(13.33) 



where the functions fa} [z,gu{y)\ are defined by formulae (jl3.1ip - (jl3.13p with 

9u+nAy) = giinAz,y) ^ Ju+nAy)Y.^^\z) - y,+„,(y)4")(z). 



;i3.34) 



The VEV given by formula ()13.33p is divergent and needs some regularization procedure. To 
make it finite we can introduce a cutoff function ^x{z), z = Ju^^s, with the cutoff parameter A, 
which decreases sufficiently fast with increasing z and satisfies the condition $a — ^ 1; A — > 0. 

To evaluate the VEV of the EMT we need to sum series over the zeros of the function 
Cl{r],z). A summation formula for this type of series can be obtained by making use of the 
GAPF. In the GAFF as functions g{z) and f{z) we choose 



r(2f>). 



Hi'''\z) 
with the sum and difference 

9(z) - ( 



Hz) 



Hi''''\r]z) h{z) 
Hi'"'\z) Ct{v,z)' 



Hz) 



Hi"'\z)Hi"'\z) 



;i3.35) 



n 



1,2. 



The conditions for the GAPF written in terms of the function h(z) are as follows 

\h{z)\ < ei{x)e'^^^^ , |z| ^ cxd, z = x + iy, 



(13.36) 



;i3.37) 



where c < 2{rj — 1) and ei{x)/x ^ for j; ^ +00. To find the residues of the function g{z) at 
the poles z = jl^^ we need the derivative 



d_ 

dz 



CUrj,z) 



Ji'\^z) 



Z = lu 



;i3.38) 
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By using this relation it can be seen that 



Res5(^) = ^Tr(r?,7L)/i(7L 



;i3.39) 



Let h{z) be an analytic function for Rez ^ except possible branch points on the imaginary 
axis. By using the GAPF, it can be seen that if this function satisfies condition (jl3.37p . 



and the integral 



exists, then [55] 



h{ze''') = -h{z) + o{z-^), z^O, 
h{x)dx 



I 

Jo 



j('^)2(x) + yr^(x) 



r>(a)2. 



;i3.40) 

;i3.4i) 



lim 

L— >+oo 



•2 



s=l 



h{x)dx 



vr 

—Res 
2 2=0 



h{z)Hi^'\r,z) 



j('^)2(x)+yj'^)'(x) 

POO 

-Y^J^ dxnii\x,vx)h{xe^-^ 



/3==f 



Cl{v,z)Hi"'\z 

where 7^^ < L < ^l^+i- Here the function i^i^J {x,t]x) is defined as 



;i3.42) 



Ki'"'\x)li"'\vx) - &^\x)Ki'^\vx) ' 



;i3.43) 



and for a given function -F(2;) we use the notation 

F^^^\z) = zF'iz) + [n("')(/i^ - ^z^e^-' + fil) - n^v]F{z) 



;i3.44) 



Now we apply to the sum over s in (I13.33P the summation formula ()13.42|) . As a func- 
tion h{z) in this formula we take h{z) = flji^^lz , g"u} {z , zr / a)]<^ \{z) . The function f{z) = 
fai^a[z,gu} {z, zr /a)] satisfies the relation 



f{xe 



,vri/2\ 



-f{xe 



-7ri/2\ 



, for ^ X ^ . 



;i3.45) 



By taking into account that for these values x one hasF("'+)(tt>2;/a) = F^''^-\wx/a), we conclude 
that in this case the part of the integral on the right of Eq. (|13.42p over the interval (0, ^q) 
vanishes after removing the cutoff. As a result the components of the vacuum EMT can be 
presented in the form [55] 

{0|I^'=|0) = (7;'=)i„ + (7^^),,, 



with separate parts 



j=l/2 a=0,l-^0 



^J^l[x,g)^J{x,xr/a)] 



{a)2. 



and 



(^i )ab 



1=1 /3=± 



Y^lT^ / dxnilHax,bx)Fif[ax,Gi-fHax,rx)], 



;i3.46) 



;i3.47) 



;i3.48) 



55 



where z/ = z/i = //a + l/2, and the notation p3.44|) is specified to 



= zF'{z) + {u + n^^'^fiy, T m^'^V^' - )^(^) • 



(13.49) 



In formula (|13.48p we have introduced the notations 



- M^ + PiM)Gi^'^^^{ax,y) 
-{Vx'^-M^ - piM)G^J-fl^ {ax, y) 



X 



+ - 



2 _M2 

2v - 1 



G^f^\ax,y)-G^:^\\ax,y) 



y 

(2v - l)x^ 



G^f\ax,y)G^:!l{ax,y) 



2yyjli 



M2 



G'^:^\ax,y)G^:!^l{ax,y) 



where 



= U{y)Ki^^\x)-K,{y)I^:"^\x),w = aA 



(13.50) 

(13.51) 
(13.52) 



(13.53) 
(13.54) 



In these formulae, for a given function F{z) we use the notation F^'^^\z) defined by formula 

As it has been shown in Ref. [5l], the term (jl3.47p presents the vacuum EMT in the case 
of a single spherical shell with radius a in the region r > a. After the subtraction of the part 
coming from the global monopole geometry without boundaries, this term is presented in the 
form 

{T^la = {Thm + {Tt)a, (13.55) 

where the part 



{Tha 



1=1 /3=± 



M 



;i3.56) 



is induced by a single sphere with radius a in the region r > a. This quantity diverges on the 
boundary r = a with the leading divergence (r — a)~^ for the energy density and the azimuthal 
stress, and (r — a)~^ for the radial stress. 
Note that by using the identities 

Il^J^pm ^ (^^^, (bx) (,^) 



+ n^'"^9.^l^l{ax,hx)F[f\x,G^^^\wx,rx)] 



(13.57) 



with the notation 



dfJ{axM) 



li'"'\ax)/ir'{bx) 



;i3.58) 



Kl''^\ax)li'^\bx) - ll^f'\ax)Kl'^\bx) ' 
the vacuum EMT in the region between the spheres can also be presented in the form 

{Tt) = {Tt)m + {Thb + {Thba, (13.59) 
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where 



{Thba = :^^^lYl / dxnif^\ax,bx)Fif[x,G^^^\bx,rx)], (13.60) 
" 1=1 p=±Jm 

and the quantities {Tf)^ are the vacuum densities induced by a single shell with radius b in the 
region r < b. The formula for the latter is obtained from ()13.56p by the replacements a ^ b, 
I K and coincides with the result derived in the previous subsection. The surface divergences 
in vacuum expectation values of the EMT are the same as those for a single sphere when the 
second sphere is absent. In particular, the term (T^^^^^ i{'Fi)ba) ^iiite on the outer (inner) 
sphere. In the formulae above taking a = 1 we obtain the corresponding quantities for a spinor 
field in the Minkowski bulk. In this case v = 1 + 1/2 and the Bessel modified functions are 
expressed in terms of elementary functions. 

As in the case of a scalar field, the vacuum forces acting on the spheres can be presented in 
the form of the sum of self-action and interaction parts. The latter is determined by formulae 
(113. 48p and (113. 60p with i = k = 1 substituting r = a and r = b respectively. By making use 
of the properties of the modified Bessel functions, the interaction force on the sphere r = w is 
presented in the form 

- —^l / dx 

^(int) ^2«2^2 9y;A^ J^^ - M'^ 

where w = a,b, and defined after formula (jl3.25p . As it has been shown in [55], these 

forces can also be obtained from the Casimir energy by standard differentiation with respect to 
the sphere radii. 



ir\ax)Kr\bx) 
Kjf^\ax)li'''\bx) 



(13.61) 



14 Electromagnetic Casimir densities for conducting spherical 
boundaries 

14.1 Energy-momentum tensor inside a spherical shell 

In this section we consider the application of the GAPF to the investigation of the VEV for the 
EMT of the electromagnetic field induced by perfectly conducting spherical boundaries in the 
Minkowski spacetime. For the region inside a perfectly conducting sphere with radius a, in the 
Coulomb gauge the corresponding eigenfunctions for the vector potential are presented in the 
form 

{ V X [ji{ujr)Xim[0,ip)], A = l 

where A = and 1 correspond to the spherical waves of magnetic and electric type (TE and TM 
modes respectively). They describe photon with definite values of the total angular momentum 
i, its projection m, energy io and parity (—1)'+'^'*"^. Here the spherical vector harmonics have 
the form 

Xim{e,^) = -i^^=^l=Yi^{0,^), //O, (14.2) 
V H' + Ij 

with Yim{9, (p) being the spherical functions, and ji{x) = \J j^xJiJ^xjiix) is the spherical Bessel 
function. The coefficients fia are determined by the normalization condition 



/ 



27r 

dV K,-Kl, = —b,,,, (14.3) 
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where the integration goes over the region inside the sphere. Using the standard formulae for 
the spherical vector harmonics and spherical Bessel functions (see, for example, |118j ). for the 
normalization coefficient one finds 

(31 = 8r^(wa)/wa, u = 1 + 1/2, (14.4a) 

where Ty[z) is defined in ()3.8p . 

Inside the perfectly conducting sphere the photon energy levels are quantized by the standard 
boundary conditions: 

nxE = 0, n-B = 0, r = a, (14.5) 

where E and B are the electric and magnetic fields and n is the normal to the boundary. They 
lead to the following eigenvalue equations with respect to uj: 

d^[rji{ur%^^ = Q, A = 0,1. (14.6) 

It is well known that these equations have an infinite number of real simple roots. 
By substituting the eigenfunctions into the mode-sum formula 

(0|r,fc|0) = ^r,fc{A,,A:}, (14.7) 

with the standard expression for the electromagnetic EMT and after the summation over m by 
using the standard formulae for the vector spherical harmonics (see, for example, [118j ). one 
obtains 

(0|7;*^|0) = diag (e, -p, -p^, -p±) (14.8) 

(index values 1,2,3 correspond to the spherical coordinates r,9,(p with the origin at the sphere 
center). Here the energy density, e, effective pressures in transverse, p±, and radial, p, directions 
are determined by the relations 

'?("'^) = -^^''T,{Loa)Di'\ji{u;r)], q = e, p, p±, (14.9) 

ujlX 

where the following notations are introduced 

and p = e — 2p± . In the sum (jl4.9p , uj takes a discrete set of values determined by equations 
(flAHD . 

The VEVs ()14.9p are infinite. The renormalization of (0|Tjfc|0) in flat spacetime is affected 
by subtracting from this quantity its singular part {OMlTiklOu), which is precisely the value it 
would have if the boundary is absent. Here |0m) is the amplitude for the Minkowski vacuum 
state. To evaluate the finite difference between these two infinities we will introduce a cutoff 
function ipfj_{uj), which makes the sums finite and satisfies the condition ~^ 1' ^• 

After the subtraction we will allow /x — > 0: 

{Tik)ron = lim [(0|r,fc|0) - (OM|Tifc|OM)] . (14.11) 

Hence, we consider the following finite quantities 

^(^'«'^) = E^il E Ei'f ^^(i'2)^M(iiVa)i?{^^bKji>)], (14.12) 
1=1 A=0,lfc=l 
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where x = r/a, and lo = jl^l/a are solutions to the eigenvalue equations (jl4.6p for A = 0, 1, 
respectively. The summation over k in (jl4.12l) can be done by using formula (I3.22P and taking 
A = 1, B = iov TE modes (A = 0) and ^ = 1, 5 = 2 for TM modes (A = 1) (recall that 
in (j3.18p Xu^k are zeros of Ju{z) with barred quantities defined as (I3.1jl ). Let us substitute in 
formula (fA]l 

fiz)=z^^^{z/a)Dl'^\jlizx)], (14.13) 

with Dl^^fiy)] defined from (nTTOll . We will assume the class of cutoff functions for which 
the function ()14.13p satisfies the conditions for Theorem 2, uniformly with respect to ^ (the 
corresponding restrictions for ip^ can be easily found from these conditions using the asymptotic 
formulae for the Bessel functions). Below for simplicity we will consider the functions without 
poles. In this case, function (jl4.13p is analytic on the right-half plane of the complex variable 
z. The discussion on the conditions to cutoff functions, under which the difference between 
divergent sum and integral exists and has a finite value independent any further details of cutoff 
function, can be found in |TT]. For TE and TM modes by choosing the constants A and B as 
mentioned above one obtains 

oo ^ — , 

dwLU^iP^{Lu)Dl'\ji{u;r)] 



where the functions -F/^^[/(y)] are defined as 



Ft\f{y)] = f'\y) + m + i)/y'-i]f{y), (i4.i5) 
Fh\f(y)] = + x,{y) = [M^y) + M-w)]/2, (u.ie) 

and F/^^[/(y)] = F^'^fiy)] - 2f/^^^ [/(y)]. In (fTmp we have introduced the Ricatti-Bessel 
functions of the imaginary argument, 

si{z) = ^I,{z), ei{z) = ^K,{z),u = l + 1/2. (14.17) 

As {^M\Tik\^M) = li™a-*oo(0|Tjfc|0), the first integral in (|14.14p represents the vacuum EMT for 
empty Minkowski spacetime: 

POO 

^ 1=1 •'^ 

This expression can be further simplified. For example, in the case of the energy density one 
has 

^ oo 

= T2Y. duu^ij^iu;) [/jf+i(a;r) + (/ + l)jti(u;r) + {21 + l)jfiujr)] 
^ 1=1 

2^ /-oo POO 

= ^ duJUJ^i>^[uj)Y,{2l + l)jf{i^r) = duu'^^ioj). (14.19) 



1=0 



As we see, the use of the GAPF allows to extract from infinite quantities the divergent part 
without specifying the form of the cutoff function. Now the renormalization of the EMT is 
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equivalent to omitting the first summand in (|14.14p . which corresponds to the contribution of 
the spacetime without boundaries. For r < a the second term on the right (jl4.14p is finite in 
the hmit ^ ^ and for the renormalized components one obtains 



oo 



21 + 1 

1=1 



dz z 



ei{az) ^ e[{az) 



si{az) s[{az) 



Fi''\si{zr)l r<a. 



[U.20) 



From here it is obvious the independence of the renormahzed quantities on the specific form 
of the cutoff, on the class of functions for which (jl4.13p satisfies conditions for (j3.22p . The 
derivation of the vacuum densities (jl4.20p given above uses the GAPF to summarize mode- 
sums and is based on [3H[35]. One can see that these formulae for the case of exponential cutoff 
function may also be obtained from the results of [1051 HQS] , where the Green function method 
is used. 

We obtained the renormalized values (jl4.20p by introducing a cutoff function and subse- 
quently subtracting the contribution due to the unbounded space. The GAPF in the form 
()3.22p allows to obtain immediately this finite difference. It should be noted that by using the 
GAPF in the form p.lSp we can derive the expressions for the renormalized azimuthal pressure 
without introducing any special cutoff function. To see this note that for x < 1 the function 
(114. ISp with q = p± and ip^ = 1 satisfies conditions of Theorem 2. It follows from here that we 
can apply formula (j3.18p directly to the corresponding sum over uj in (jl4.9p or over k in (jl4.12p 
without introducing a cutoff function. This immediately yields to formula (114. 20p for q = p± 
with ■0„ = 1. 



The VEV given by formula (jl4.20p diverges on the sphere surface, r = a, and this divergence 
is due to the contribution of large I (note that the integral over z converges at r = a). The 
corresponding asymptotic behavior can be found by using the uniform asymptotic expansions 
for the Bessel functions and the leading terms have the form 



2p± 



-1 



-1 



307r^a(a — r)^ ' 



607r^a^(a — r) 



2 ■ 



(14.21) 



These surface divergences originate in the unphysical nature of perfect conductor boundary con- 
ditions and are well-known in quantum field theory with boundaries. As we have mentioned 
before, they are investigated in detail for various types of fields and general shape of the bound- 
ary. Eqs. (jl4.2ip are particular cases of the asymptotic expansions for the VEV near the smooth 
boundary given in [82]. In reality the expectation values for the EMT components will attain 
a limiting value on the conductor surface, which will depend on the molecular details of the 
conductor. From the asymptotic expansions given above it follows that the main contributions 
to 

qren{i") sre duc to the frequencies w < (a — r) ^. Hence we expect that formulae (114. 20p 
are valid for real conductors up to distances r for which (a — r)~^ <^ ujq, with ujq being the 
characteristic frequency, such that for uj > ujq the conditions for perfect conductivity fail. 

At the sphere center, in (|14.20p I = 1 multipole contributes only and we obtain \35\ I103| 



e(0) = 3p(0) =3px(0) 



-0.0381a- 



(14.22) 



At the center the equation of state for the electromagnetic vacuum is the same as that for black- 
body radiation. Note that the corresponding results obtained by using the uniform asymptotic 
expansions for the Bessel functions [105^ I106| are in good agreement with (114. 22p . 

The components of the renormalized EMT satisfy continuity equation V^T^^ = 0, which for 
the spherical geometry takes the form 



p'{r) + 2{p — p±)/r = 0. 



(14.23) 



60 



From here by using the zero trace condition the following integral relations may be obtained 

1 r 2 r 

p{r)= dtt^e{t) = ^ dttp^it), (14.24) 

r'^ Jo Jo 

where the integration constant is determined from relations ()14.22p at the sphere center. It 
follows from the first relation that the total energy within a sphere with radius r is equal to 

E{r)=4TT[ dtt'^e{t) = 3V{r)p{r), (14.25) 
Jo 

where V{r) is the corresponding volume. 

The distribution for the vacuum energy density and pressures inside the perfectly conducting 
sphere can be obtained from the results of the numerical calculations given in [10511106] . In their 
calculations Brevik and Kolbenstvedt used the uniform asymptotic expansions of the Ricatti- 
Bessel functions for large values of the order. In [33| [35] (see also j28j ) the corresponding 
quantities are calculated on the base of the exact relations for these functions and the accuracy 
of the numerical results in |105[ 1106] is estimated (~ 5%). The simple approximation formulae 
are presented with the same accuracy as asymptotic expressions. Note that inside the sphere all 
quantities e, p, p± are negative. 

By the method similar to that used in this subsection, the VEVs for gluon fields can be 
evaluated in the bag model for hadrons [119] . In this model the vacuum outside the spherical 
bag is a perfect color magnetic conductor. The color electric and magnetic fields are confined 
inside the bag and satisfy the boundary conditions n • E = 0, n x B = on its boundary (for 
the role of the Casimir effect in the bag model see O HI [7] ) . 

14.2 Electromagnetic vacuum in the region between two spherical shells 

Now let us consider the electromagnetic vacuum in the region between two concentric conducting 
spherical shells with radii a and b, a < b. In the Coulomb gauge the complete set of solutions 
to the Maxwell equations can be written in the form similar to (114. ip : 

A. = /3.^| ^ ^ = ° (14.26) 

\ V x[gii{uja,ujr)Xi^], A = l ' 



where as above the values A = and A = 1 correspond to the waves of magnetic (TE modes) 
and electric (TM modes) type, 

/ N f ilMn-iix) — ji(x)ni(y), A = , 

{ Jl{y)[xni{x)] - [xjiix)] ni{y), A = l 

with ni{x) being the spherical Neumann function. From the boundary conditions at surfaces 
r = a and r = b one finds that the possible energy levels of the photon are solutions to the 
following equations 

[rgxi{uja,Lor)l^^ = 0, A = 0, 1. (14.28) 

All roots of these equations are real and simple [65]. 

The coefficient Po- in (|14.26p is determined from the normalization condition (114. 3p . where 
now the integration goes over the region between spherical shells, a ^ r ^ b. By using the 
standard relations for the spherical Bessel functions it can be presented in the form 



Pi 



bjfiub) 

1 j b[ujaji{uja)] ^ 
1 a [ujbji{ujb)]"^ 



A = 0, (14.29) 
/(/ + 1) 



CO 



252 



N -1 

1 + ^^^i , A = l. (14.30) 
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From the mode-sum formula (|14.7|) with the functions (|14.26p as a complete set of solutions one 
obtains the VEV in the form (jl4.8p with 



« = ^E(2^ + l)^^^"A^'^[5AK^a,u;r)], q = e,p,p^, (14.31) 

U)l\ 

where the frequencies uj are solutions to equations (|14.28p and the functions D'f^ [/(y)] are defined 
by formulae (114. lOp with /(y) = gxi{uja, y). It is easy to see that the eigenvalue equations (114. 28p 
can be written in terms of the function C"*, defined by ()4.ip . as 

Cf{r],uja) = 0, u = 1 + 1/2, rj = b/a, A = 1/{1 + X), B = X, X = 0,1. (14.32) 

By this choice of constants Aj and Bj the normalization coefficients (114. 29p and (jl4.30p are 
related to the function T^^ from (14. 7p by the formula 

Pl = Tf{r,,u;a). (14.33) 

This allows us to use the formulae from Section H] for the summation over the eigenmodes. 

As above, to regularize infinite quantities p4.3ip we introduce a cutoff function Tp^{ui) and 
consider the finite quantities 

« = E ^EE^'^-^''(^'^)^M(^/«)A^'^[5Az(^,^rr)]|^^^(.), (14.34) 

1=1 X=Ok=l "'^ 

where x = r/a, and uja = 7^'*^ are solutions to equations (|14.28p or ()14.32p . In order to sum 
over k we will use formula (|4.14p with 

h{z) = z^^^{z/a)D\''\gxi{z,zx)], (14.35) 

assuming a class of cutoff functions for which p4.35p satisfies conditions (|4.5p and (|4.12p uni- 
formly with respect to fi. The corresponding restrictions on -0^ can be obtained using the 
asymptotic formulae for Bessel functions. From (jl4.34p . by applying to the sum over k formula 
(|4.14p . for the EMT components one obtains 



1=1 A=0,1 



1 r e\»(vz)zx,{zla)F','\a>.,{z.zx)] \ 

where we use the notations 

4'\y) ^ d^eiiy), s^'^y) ^ d^si{y) (14.37) 

Gxiix,y) = 4^\x)siiy)-eiiy)s[^\x), (14.38) 

for the Riccati-Bessel functions, 

^h\l{Z) -Sr.. Ni'2 , r / m'2 \ 1 ' (14.39) 

and the functions i^/''^[/(y)] with f{y) = G\i{z,y) are defined by relations ()14.15p . (114. 16p . 
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In p4.36|) . taking the limit b ^ oo for fixed a and r, we obtain the VEV for the EMT 
components outside a single conducting spherical sheh with radius a. In this limit the second 
integral on the right of formula p4.36|) tends to zero, whereas the first one does not depend on 
b. Hence, one obtains 



1=1 A=0,1 



dz z'^ipf^{z/a) 



Dl''\gxi{z,zx)] 



(14.40) 



For the renormahzation of expressions (|14.40p we subtract the Minkowskian part, namely ex- 
pression ()14.18p . It can be seen that 



Dl''\gxiiz,zx)] (,) 



E 

m=l,2 



2 ^^d^[zh'i-\z)] 



D^ihrizx)], 



(14.41) 



with hj"^\z), m = 1,2 being the spherical Hankel functions. Now, for the corresponding z- 
integral we rotate the integration contour in the complex z-plane by the angle 7r/2 for the term 
with m = 1 and by the angle —tt/2 for the term with m = 2. By introducing the Ricatti-Bessel 
functions (114. ITh and for points with r > a removing the cutoff, for the renormalized components 
of the vacuum EMT outside the sphere we find 



(a,r) 



1=1 



87r2r2 



sijaz) ^ s'ljaz) 
ei{az) e[[az) 



Fi''\ei{zr)], r 



> a, 



;i4.42) 



where the functions F'^'^\f{y)\ are defined by formulae (jl4.15p . (jl4.16p . The exterior mode-sum 
consideration given in this section follows [MIES]- For the case of exponential cutoff function 
formula (jl4.42p can also be obtained from the results |105^ 1106) . where the Green function 
formalism was used. Note that the expressions for the exterior components are obtained from 
the interior ones replacing si <^ e;. 

Expressions (114. 42p diverge on the sphere. The leading terms of these divergences may be 
found using the uniform asymptotic expansions for the modified Bessel functions for large values 
of the order and these terms are given by the same formulae as those for the interior region (see 
(114. 2ip ). In particular, the cancellation of interior and exterior leading divergent terms occurs 
in calculating the total energy and force acting on the sphere. The same cancellations take place 
for the next subleading divergent terms as well. For distances far from the sphere one finds 



dz z^e\{z) 



IGvr^r''' ' 



-4p 



27r2r7' 



r a. 



;i4.43) 



The results of the numerical calculations for the vacuum EMT components outside the sphere are 
given in |35[ 1105^ I106j . In [105^ I106j calculations are carried out by using the uniform asymptotic 
expansions for the Riccati-Bessel functions. The accuracy of this approximation is estimated in 
[35], where exact relations are used in numerical calculations. Simple approximating formulae 
with the same accuracy as those for the asymptotic calculations are presented as well. In the 
exterior region the energy density and azimuthal pressure are positive, and radial pressure is 
negative. 

33|) may now be written in the following integral form 



Note that the continuity equation 



pir) 



1 



dtfe{t) 



dttp±{t). 



[U.U) 
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Prom (|14.24|) and (|14.44p it follows that 

E{a) = j dVe{r) = A'Ka^[p{a-) - p{a+)l 



;i4.45) 



where E{a) is the total vacuum energy for a spherical shell with radius a, p(a±) = lim^^o p{a^r). 
By using the expressions for p{r) given above, one can obtain the following formula for the total 
energy (the same result can also be obtained by integrating the energy density) 



E{a) 



E 

1 = 1 

OO 

E 

1=1 



2vra 7o SAz)e\{z) 



s\{z)e\{z) ^ + ^ ^ 
si{z)ei{z) z^ 



2/ + 1 
27ra 



dzXfiiz/a)z-^ln^l - [si{z)ei{z)] ^| , 



;i4.46) 



where we have restored the cutoff function. By taking the cutting function ip^{oo) = e~'^ one 
obtains the expression for the Casimir energy of the sphere derived in [99j by the Green function 
method. Note that in this method the factor ■0^(zz/a) = e~^^'^ appears automatically as a result 
of the point splitting procedure. The evaluation of ()14.46p leads to the result E = 0.092353/2a 
for the Casimir energy of a spherical conducting shell [HTl [971 [98t [99 l[T00l 1101 ]. This corresponds 
to the repulsive vacuum force on the sphere. 

As we have seen, the VEV of the EMT in the region between two concentric perfectly 
conducting surfaces is given by ()14.36p . Using this formula the components of the vacuum EMT 
can be presented in the form 



/(a,6,r) = g(a,r) +g('''')(r). 



e,p±,p, a<r<b, 



;i4.47) 



where q{a, r) are the corresponding quantities outside a single sphere of radius a given by (114. 40p . 
and [371 



1=1 



A=0.1 



- ^ Fl''\Gxi{az,Tz)]eY>{hz)/e\-'\az) 
dz z- 



{az)s'l^\bz) 



{bz)s\^\az 



(14.48) 



In (jl4.47p the dependence on b is contained in the summand q^""^^ (r) only. This quantity is finite 
for a ^ r < b and the renormalization of q{a, b, r) is equivalent to the renormalization of the 
first summand. 

It can be seen that the quantities ()14.47|) may also be written in the form 



{a,b,r) = q{b,r) + q^''^\r), q = e,p±,p, a < r < b, 



where 



g(H(r) 



^ 21 + 1 



A=0,1 



- Fi''\G,i{bz,rz)]sr{az)/sr{bz) 

(A). X (A)., X (A),, X (A). x' 

el [az)Si [bz) — el {bz)sl [az) 



(14.49) 



(14.50) 



In ()14.49p . q^^°'\r) — > when a — > and q{b,r) coincides with the corresponding quantities 
inside a single conducting shell with radius b. Note that in (jl4.50p the sum and integral are 
convergent for a < r ^b. 

As we have already mentioned the vacuum energy density and stresses inside and outside a 
single shell may be obtained from the expressions q{a, b, r) in limiting cases a ^ or 6 — > cxd. It 
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can be seen that in the hmit a,b ^ oo with fixed h = b — a, we obtain the standard result for 
the Casimir parallel plate configuration with e = — 7r^/720/i^. 
Let us present the quantities q = £,p,p± in the form 

q = q{a,r) + q{b,r) + Aq{a,b,r), a < r < b, (14.51) 

where the interference term may be written in two ways 

Aq{a,h,r) = q^"^'^ {a,b,r) - q{b,r) = q^^°'\a,b,r) - q{a,r). (14.52) 

This term is finite for all a ^ r ^ 6. Near the surface r = a it is convenient to use the first 
presentation in (jl4.52l) . as for r — > a both summands in this formula are finite. For the same 
reason the second presentation is convenient for calculations near the surface r = b. 

So far we have considered the electromagnetic vacuum in the region between two perfectly 
conducting spherical surfaces. Consider now a system consisting of two concentric thin spherical 
shells with radii a and b, a < b. In this case the VEV for the EMT components may be written 
in the form 

q{a, b, r) = q{a, r)9{a - r) + q{b, r)e{r -b) + [q{a, r) + (r)]6'(r - a)e{b - r), (14.53) 

where 9{x) is the unit step function. By using the continuity equation for the EMT, it is easy 
to see that the total Casimir energy for the system under consideration can be presented in the 
form 

E{a, b) = E{a) + E{b) + Att [ft^^M (5) _ a^pi^h) (^)l ^ ^^^ r^^^ 

where E{j) is the Casimir energy for a single sphere with radius j = a,b. As follows from 
(|14.48p and ()14.50p , the additional vacuum pressures on the spheres are equal to |28l [37] 



' ~ h'-'b'Jo '^p)(a./6).f)(.)-e;^)(.).!^)(a./5)•^''•''^ 

In [120j . the vacuum forces acting on boundaries in the geometry of concentric conducting 
spherical shells are investigated by making use of local Green function method. The results 
of the numerical evaluations of quantities Aq{a,b,r), q = £,p,p±, as well as those for p^"'^\a), 
p(*»)(5) are presented in \28\ I37j . Note that, as follows from the results of these calculations, 
the quantities (114.551) and (114. 56p are always negative, and therefore the interaction forces 
between two spheres are always attractive (as in the parallel plate configuration). Note that the 
interaction forces can also be obtained from the corresponding part of the total Casimir energy 
in the region between the spheres by standard differentiation with respect to the sphere radii 



15 Vacuum polarization induced by a cylindrical boundary in 
the cosmic string spacetime 

Cosmic strings generically arise within the framework of grand unified theories and could be 
produced in the early universe as a result of symmetry braking phase transitions |12H I122j . 
Although the recent observational data on the cosmic microwave background radiation have 
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ruled out cosmic strings as the primary source for primordial density perturbations, they are 
still candidates for the generation of a number interesting physical effects such as the genera- 
tion of gravitational waves, high energy cosmic strings, and gamma ray bursts. In the simplest 
theoretical model describing the infinite straight cosmic string the spacetime is locally flat ex- 
cept on the string where it has a delta shaped curvature tensor. In quantum field theory the 
corresponding non-trivial topology leads to non-zero VEVs for physical observables. In this sec- 
tion we will consider the vacuum polarization induced by a cylindrical boundary coaxial with a 
cosmic string assuming that on the bounding surface the field obeys Robin boundary condition 
|43j . The cylindrical problem for the electromagnetic field in the Minkowski bulk with perfectly 
conducting conditions was first considered in |123j (see also [TJ 11241 1125) ). While the earliest 
studies have focused on global quantities, such as the total energy and stress on a shell, the 
local characteristics of the corresponding electromagnetic vacuum are considered in [39] for the 
interior and exterior regions of a conducting cylindrical shell, and in [40j for the region between 
two coaxial shells (see also |28j). The vacuum forces acting on the boundaries in the geometry of 
two cylinders are also considered in |126j . The vacuum densities for a Robin cylindrical bound- 
ary in the Minkowski background are investigated in [41 j . In [127[ a cylindrical boundary with 
Dirichlet boundary condition is introduced in the bulk of the cosmic string as an intermediate 
stage for the calculation of the ground state energy of a massive scalar field in (2-|-l)-dimensions. 

15.1 Wightman function 

We consider a scalar field (p with curvature coupling parameter Q propagating on the {D + 1)- 
dimensional background spacetime with a conical-type singularity described by the line element 

ds^ = Qikdx'dx^ = dt^ - dr^ - r^dcp"^ - ^ dzi^, (15.1) 

i=l 

with the cylindrical coordinates (x^, x^, . . . , x^) = {r, (j), zi, . . . , Z]\j), where N = D — 2, O^cj)^ 
(po, and the spatial points (r, 0, zi , . . . , zj\f) and (r, (j) + zi, . . . , zjy) are to be identified. In 
the standard D = 3 cosmic string case the planar angle deficit is related to the mass per unit 
length of the string by 27r — 0o = SnGfi, where G is the Newton gravitational constant. We 
assume that the field obeys Robin boundary condition (|9.1ip on the cylindrical surface of radius 
a, coaxial with the string. 

In the region inside the cylindrical surface the eigenfunctions satisfying the periodicity con- 
dition are specified by the set of quantum numbers a = (n,7,k), n = 0, ±1, ±2, . . ., k = 
(ki, . . . , kisf), — oo < kj < oo, and have the form 

(faix) = (3aJg\n\{7r)exp{iqn(j) + i'kr^i-iujt) , (15.2) 
a; = VT^ + fe^^, k'^ = k"^ + w? , q = 2tt/(1>o, (15.3) 

where ry = {zi, . . . , z^). The eigenvalues for the quantum number 7 are quantized by boundary 
condition (I9.11|) on the cylindrical surface r = a. From this condition it follows that for a given 
n the possible values of 7 are determined by the relation 

7 = A„,,/a, / = 1,2,..., (15.4) 

where A„^/ are the positive zeros of the function Jg|„,|(-z), Jq\n\{^n,i) = 0, arranged in ascending 
order, Xn^i < Xn,i+i, n = 0, 1, 2, . . .. Here, for a given function f{z), we use the barred notation 
p.ip with the coefficients 

A = A, B = -B/a. (15.5) 
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In the following we will assume the values oi A/B for which all zeros are real. 

The coefficient Pa- in (|15.2p is determined from the normalization condition with the inte- 
gration over the region inside the cylindrical surface and is equal to 

^'^ " (27r)^^,/.oa2' ^^^-^^ 

with the notation Tyi^z) from (j3.8p . Substituting the eigenfunctions (jl5.2|) into the mode-sum 
formula (j9.8p with a set of quantum numbers a = (n/k), for the positive frequency Wightman 
function one finds 

W{x,x') = 2 / d^ke*''^''!! ^'cos(gnA(^) 

n=0 

oo 

X Yl l^lJmiir)JUir')e-''^''%=x^,/a, (15.7) 
1=1 

where Ary = ry — rj|, A(/> = (p — cp' , At = t — t' , and the prime in the sum means that the 
summand with n = should be taken with the weight 1/2. As we do not know the explicit 
expressions for the eigenvalues X^^i as functions on n and /, and the summands in the series over 
I are strongly oscillating functions for large values of /, this formula is not convenient for the 
further evaluation of the VEVs of the field square and the EMT. In addition, the expression 
on the right of (jl5.7p is divergent in the coincidence limit and some renormalization procedure 
is needed to extract finite result for the VEVs of the field square and the EMT. To obtain an 
alternative form for the Wightman function we will apply to the sum over / summation formula 
p.lSp . In this formula as a function f{z) we choose 

zJqn{zr/a)Jgn{zr' /a) 



m = ^ "^^z ^ ^ eM-^^Wkl + zya^), (15.8) 

where /c^ = k'^ + m?. The condition to formula (jS.lSp to be satisfied is r + r' + |At| < 2a. In 
particular, this is the case in the coincidence limit t = t' for the region under consideration, 
r, r' < a. Formula (jS.lSp allows to present the Wightman function in the form [13] 



W{x, x') = Ws{x, x') + {ip{x)^{x'))a, (15. 

where 



Ws{x,x') = ^ I T^^z^e^^^^ii / dz 



oo 

cos{qnA(j))Jgnizr)Jgnizr'), (15.10) 



X 
n=0 



and 



oo 



TTcPoJ (27r)^ yi^^ 

COii[qnAS(p)lqn{Zr)lgn{Zr 

n=0 



m 



X V'cos(gnA0)V(zr)V(zrO ^7^7 . (15.11) 

^ Iqn{za) 



In the limit a — oo for fixed r, r', the term {(p{x)ip{x'))a vanishes and, hence, the term Ws{x, x') 
is the Wightman function for the geometry of a cosmic string without the cylindrical boundary 
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(the integral representation of the corresponding Green function for a massive scalar field is 
considered in Refs. |128j ). Consequently, the term {ip{x)ip{x'))a is induced by the presence of 
the cylindrical boundary. For points away from the cylindrical surface this term is finite in the 
coincidence limit and the renormalization is needed only for the part coming from (jl5.10p . As 
it has been shown in Ref. [33], the boundary induced part of the Wightman function in the 
exterior region, r > a, is obtained from the corresponding part in the interior region by the 
replacements I K. 

15.2 Vacuum expectation values of the field square and the energy-momentum 
tensor inside a cylindrical shell 

Taking the coincidence limit rc' — > x in formula (|15.9p for the Wightman function and integrating 
over k, the VEV of the field square is presented as a sum of two terms: 

where jO^) is the amplitude for the vacuum state in the geometry when the cylindrical shell is 
absent. The second term on the right of this formula is induced by the cylindrical boundary and 
is given by the formula [33] 

(^')^ = -4^f;' rdzz{z^-m^)'^^S^lUzr), (15.13) 
where we have introduced the notation 

= YpII)- ^^^-^^^ 

The formula for the VEV of the field square in the region outside the cylindrical shell is obtained 
from ()15.13p by the replacements I ^ K. For points away from the cylindrical surface, r < a, 
the integral in (jl5.13p is exponentially convergent in the upper limit and the boundary-induced 
part in the VEV of the field square is finite. In particular, this part is negative for Dirichlet 
scalar and is positive for Neumann scalar. Near the string, r <C a, the main contribution to 
{^'^)a comes from the summand with n = and one has 

{^\^-y4^ r dzz{z^-M)^B^. (15.15) 

As the boundary-free renormalized VEV diverges on the string, we conclude from here that near 
the string the main contribution to the VEV of the field square comes from this part. 

The part (v'^)a diverges on the cylindrical surface r = a. Near this surface the main contribu- 
tion into (jl5.13p comes from large values of n. Introducing a new integration variable z — > nqz, 
replacing the modified Bessel functions by their uniform asymptotic expansions for large values 
of the order (see, for instance, [65]), and expanding over a — r, to the leading order one finds 



« ^oir/' \ (15.16) 



(1 - 25bo)T 



(47r) 2 (a — r 



This leading behavior is the same as that for a cylindrical surface of radius a in the Minkowski 
spacetime. As the boundary-free part is finite at r = a, near the boundary the total renormalized 
VEV of the field square is dominated by the boundary-induced part and is negative for Dirichlet 
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scalar. Combining this with the estimation for the region near the string, we come to the 
conclusion that in this case the VEV of the field square vanishes for some intermediate value of 
r. 

Now we turn to the investigation of the boundary- induced VEV given by (jl5.13p . in the 
limiting cases of the parameter q. Firstly consider the limit when the parameter q is large which 
corresponds to small values of 0o and, hence, to a large planar angle deficit. In this limit the 
order of the modified Bessel functions for the terms with n 7^ in (jl5.13p is large and we can 
replace these functions by their uniform asymptotic expansions. On the base of these expansions 
it can be seen that to the leading order the contribution of the terms with n 7^ is suppressed 
by the factor q^^~^^/'^{r /a)'^'^ and the main contribution to the VEV of the field square comes 
from the n = term: 



with the linear dependence on q. In the same limit the boundary-free part in the VEV of the field 
square behaves as q^~^ and, hence, its contribution dominates in comparison with the boundary- 
induced part. In the opposite limit when (7^0, the series over n in Eq. ()15.13p diverges and, 
hence, for small values of q the main contribution comes from large values n. In this case, to the 
leading order, we can replace the summation by the integration: ^„ f{qn) (1/^) dxf{x). 
As a consequence, we obtain that in the limit g ^ the boundary-induced VEV in the field 
square tends to a finite limiting value: 



Now we consider the VEV of the EMT for the situation when the cylindrical boundary is 
present. Similar to the case of the field square, this VEV is written in the form 



where the part {{)s\Tik\^s) corresponds to the geometry of a cosmic string without boundaries and 
{Tik)a is induced by the cylindrical boundary. The first term for a conformally coupled D = 3 
massless scalar field was evaluated in Ref. pHS]. The case of an arbitrary curvature coupling 
is considered in Refs. [1301 143j . The term induced by the cylindrical shell is obtained from the 
corresponding part in the Wightman function, acting by the appropriate differential operator and 
taking the coincidence limit (see formula (|9.10p ). For points away from the cylindrical surface 
this limit gives a finite result. For the corresponding components of the EMT one obtains (no 
summation over i) 




(15.17) 




(15.18) 



{0\Tik\0) = {Os\Tik\Os) + {Tik) 



(15.19) 




00 



(15.20) 



with the notations 




(15.21) 



(15.22) 



(15.23) 
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and [f{y)] = i^in [f{y)] for i = 3, . . . ,L>. The formula for the VEV of the EMT in the 
region outside the cyhndrical sheh is obtained from (I15.20p by the replacements I K. It can 
be checked that the expectation values (|15.20p satisfy the continuity equation for the EMT. The 
boundary-induced part in the VEV of the EMT given by Eq. (115. 20p is finite everywhere except 
at points on the boundary and at points on the string in the case q < 1. Unlike to the surface 
divergences, the divergences on the string are integrable. 

In the case > 1, near the string, r ^ 0, the main contribution to the boundary part (115. 20p 
comes from the summand with n = and one has 



2^00-"+^ J ma loiz) 

with the notations 

F^'\z) =2C-l+ '-f_'[^'\ F(^\z) =C-l,^ = l,2. (15.25) 

For q < 1 the main contribution to the boundary-induced part for points near the string comes 
from n = 1 term and in the leading order one has 



„2/l 2(jr~2p(«) foo 0-3 K 



za 



^ ''"^ 2'?7rr2(g + l) y„ ^ > ^ > 



za) 



where 



= g(2C - 1/2), f[^^ = C, ) = {2q - 1)C. (15.27) 



As we see, in this case the VEVs for the EMT diverge on the string. This divergence is integrable. 
In particular, the corresponding contribution to the energy in the region near the string is finite. 

As in the case of the field square, in the limit g ^ 1 the contribution of the terms with n 7^ 
to the VEV of the EMT is suppressed by the factor q^^~^^^'^(r /a)'^'^ and the main contribution 
comes from the n = term with the linear dependence on q. In the same limit, the boundary- 
free part in the VEV of the EMT behaves as q^~^^ and, hence, the total energy-momentum 
tensor is dominated by this part. In the opposite limit, when g <C 1, by the way similar to that 
used before for the VEV of the field square, it can be seen that the boundary-induced part in 
the vacuum EMT tends to a finite limiting value which is obtained from (jl5.20p replacing the 
summation over n by the integration. 

The boundary part (T^'^ ^ diverges on the cylindrical surface r = a. Introducing a new 
integration variable z nqz and taking into account that near the surface r = a the main 
contribution comes from large values of n, we can replace the modified Bessel functions by their 
uniform asymptotic expansions for large values of the order. To the leading order this gives 

This leading divergence does not depend on the parameter q and coincides with the corresponding 
one for a cylindrical surface of radius a in the Minkowski bulk. For the radial component to 
the leading order one has (Ti )a ~ (a — r)~^. In particular, for a minimally coupled scalar field 
the corresponding energy density is negative for Dirichlet boundary condition and is positive for 
non-Dirichlet boundary conditions. For a conformally coupled scalar the leading term vanishes 
and it is necessary to keep the next term in the corresponding asymptotic expansion. As the 
boundary-free part in the VEV of the EMT is finite on the cylindrical surface, for points near 
the boundary the vacuum EMT is dominated by the boundary-induced part. Taking g = 1, 
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from the formulae given in this section we obtain the corresponding results for the geometry of 
a cylindrical boundary in the Minkowski bulk [41j . 

The problem considered in this section is closely related to the problem of the investigation 
of the vacuum densities in the geometry of a wedge with the opening angle 4>o and with the 
cylindrical boundary of radius a [451I46|. For a scalar field with Dirichlet boundary condition the 
corresponding eigenfunctions in the region inside the cylindrical shell are determined by formula 
(|15.2p with the replacement exp{iqn(j)) — > sm{qn(l)), n = 1, 2, . . ., and q = ir/cpo. The eigenvalues 
for 7 are given by relation (jl5.4p . where now A„^/ are the zeros of the function Jqn{x). The 
corresponding Wightman function and the VEVs of the field square and the EMT are evaluated 
by making use of summation formula (j3.18p in the way similar to that used in this section for 
the geometry of a cosmic string. The corresponding formulae can be found in [46j. Note that in 
the case of a wedge the VEVs of the field square and the EMT depend on the angle (j). 



16 Electromagnetic Casimir densities induced by a conducting 
cylindrical shell in the cosmic string spacetime 

In this section we consider the application of the GAPE for the investigation of the polarization 
of the electromagnetic vacuum by a perfectly conducting cylindrical shell coaxial with the cosmic 
string [44]. This geometry can be viewed as a simplified model for the superconducting string, 
in which the string core in what concerns its superconducting effects is taken to be an ideal 
conductor. The background spacetime is described by the = 1 version of the line element 

dm}. 

16.1 Vacuum expectation values of the field square inside a cylindrical shell 

In the region inside the cylindrical shell we have two different types of the eigenfunctions cor- 
responding to the waves of the electric and magnetic types. In the Coulomb gauge, the vector 
potentials for these waves are given by the formulae 

A =R / (lA^^) (T^ea + iA;Vt) Jg|„|(7r)exp [z(5n(/> + A;z - wi)] , A = , , 

" ^""X -e^>iVt{J,\n\{ir)e^v[i{qn<t^ + kz-ut)]], A=l ' ^'^"'^ 

where es is the unit vector along the cosmic string, is the part of the nabla operator transverse 
to the string, and 

= ^2 _^ q = 27r/0o, n = 0, ±1, ±2, . . . . (16.2) 

Here and in what follows A = and A = 1 correspond to the cylindrical waves of the electric 
(transverse magnetic (TM)) and magnetic (transverse electric (TE)) types, respectively. The 
normalization coefficient in (116. ip is found from the orthonormalization condition (|14.3p . where 
the integration goes over the region inside the shell. From this condition, by using the standard 
integral involving the square of the Bessel function, one finds 

^ qT^^ 

TTLOa'J 

with T,y{x) defined by (fSlHD . 

The eigenvalues for the quantum number 7 are determined by standard boundary conditions 
(|14.5p for the electric and magnetic fields on the cylindrical shell. From boundary conditions we 
see that these eigenvalues are solutions of the equation 

J^^J|(7a)=0, A = 0,1, (16.4) 
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where a is the radius of the cyhndrical shell, jI)'\x) = Jy{x) and Jl'''{x) = Jl{x). We will 
denote the corresponding eigenmodes by 7a = j^} , / = 1, 2, . . .. As a result the eigenfunctions 
are specified by the set of quantum numbers a = {k, n, A, /). 

Substituting the eigenfunctions into the corresponding mode-sum formula, for the VEVs of 
the squares of the electric and magnetic fields inside the shell we find 



n=0 



+00 



,-(A)3 
,1 



A=0,1 1=1 



i7rHk,J,nUn}r/a)], 



(16.5) 



where F = E,B with rj^x = A, rj^x = 1 — A, and the prime in the summation means that the 
term n = should be halved. In (|16.5|) . for a given function f{x), we have used the notations 



(fcVVx^) [f^x) + u^f{x)/x^] + fix), j = 
(1 + fcVVx^) [/'2(x) + iy^f{x)/x^] , j = l 



[16.6) 



The expressions (|16.5|) corresponding to the electric and magnetic fields are divergent. They may 
be regularized introducing a cutoff function ipnioj) with the cutting parameter ^ which makes 
the divergent expressions finite and satisfies the condition ip^{uj) 1 for /i — > 0. After the 
renormalization the cutoff function is removed by taking the limit /i ^ 0. An alternative way 
is to consider the product of the fields at different spacetime points and to take the coincidence 
limit after the subtraction of the corresponding Minkowskian part. Here we will follow the first 
approach. 

In order to evaluate the mode-sum in (jl6.5p , we apply to the series over j summation formula 
(13.18p . As it can be seen, for points away from the shell the contribution to the VEVs coming 
from the second integral term on the right-hand side of (j3.18p is finite in the limit fi ^ and, 
hence, the cutoff function in this term can be safely removed. As a result the VEVs can be 
written as 

(0|F2|0) = {Os\F^\Os) + {F^)a, (16.7) 

where 



{Os\F^\Os 



E' 



n=0 



+00 



dk 



1 + 2- 







(161 



and 



4.q 



E7 ^^-E 



n=0 



A=0,1 



dx X" 



Kqn {xa) Gqn 



[k,Ign{xr)] 



4n\xa) 



Note that in Eq. (116.90 we used the notations 



k^ 



(16.9) 



(16.10) 



The second term on the right-hand side of Eq. ()16.7p vanishes in the limit a ^ 00. Thus, we 
can conclude that the term (O^jF^IO^) corresponds to the part in VEVs when the cylindrical 
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shell is absent with the corresponding vacuum state |0s). Note that for the geometry without 
boundaries one has (Os|-E^|Os) = (Os|-B^|Os). Hence, the application of the GAPF enables us 
to extract from the VEVs the boundary-free parts and to write the boundary-induced parts 
in terms of the exponentially convergent integrals. The boundary-free parts can be further 
simplified with the final result 



(g2-i)(g2 + n) 



ISOvrr^ 



:i6.111 



Changing the integration variable to y = ^/x"^ — and introducing polar coordinates in the 
(fc, y) plane, after the explicit integration over the angular part, the part in the VEV induced 
by the cylindrical shell can be written in the form 



vr 



E' 



A=0,1 



dx X" 



Iqn {xa) 



' qn 



Iqn{xr)] 



:i6.i2) 



where we have used the notation 



x) + v^f\ 

-nx) - 



x)/x^ + 2f{ 
v^f{x)/x\ 




1 



;i6.i3) 



The boundary-induced parts for the electric and magnetic fields are different and, hence, the 
presence of the shell breaks the electric-magnetic symmetry in the VEVs. Of course, this is a 
consequence of different boundary conditions for the electric and magnetic fields. The formulae 
for the VEVs of the field squares in the region outside the cylindrical shell are obtained from 
(116. 12p by the replacements I K. 

The expression in the right-hand side of (jl6.12p is finite for < r < a and diverges on the 



shell with the leading term {E 



(3/47r)(a — r) . Near the string, r/a <C 1, the 



asymptotic behavior of the boundary induced part in the VEVs of the field squares depends 
on the parameter q. For q ^ 1, the dominant contribution comes from the lowest mode n = 
and to the leading order one has {E'^)a ~ O.Slq/a"^, {B'^)a ~ — 0.742g/a^. For q < 1 the main 
contribution comes form the mode with n = 1 and the boundary-induced parts diverge on the 
string. The leading terms are given by 



(r/a)2(9-i) 
22(g-i)7rr2(g)a4 



dx X 



2q+l 



Kqjx) 

Iq{x) I'qix) 



K^x) 



(16.14) 



As for points near the shell, here the leading divergence is cancelled in the evaluation of the 
vacuum energy density. In accordance with (jl6.1ip . near the string the total VEV is dominated 
by the boundary-free part. Here we have considered the VEV for the field square. The VEVs 
for the bilinear products of the fields at different spacetime points may be evaluated in a similar 
way. 

Now, we turn to the investigation of the behavior of the boundary-induced VEVs in the 
asymptotic regions of the parameter q. For small values of this parameter, g <C 1, the main 
contribution into (|16.12p comes from large values of n. In this case, we can replace the summation 
over n by an integration in accordance with the correspondence 



^ h{qn) 



n=0 



1 

9 Jo 



dx h{x). 



(16.15) 



By making this replacement, we can see from (|16.12p that, in this situation, the boundary 
induced VEVs tend to a finite value. Note that the same is the case for the boundary-free part 
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(|16.1ip . In the limit q ^ 1, the order of the modified Bessel functions is large for n ^ 0. By 
using the corresponding asymptotic formulae it can be seen that the contribution of these term 
is suppressed by the factor exp[— 2gnln(a/r)]. As a result, the main contribution comes from 
the lowest mode n = and the boundary induced VEVs behave like q. 

16.2 Vacuum expectation value for the energy-momentum tensor 

In this subsection we consider the vacuum EMT in the region inside the cylindrical shell. Sub- 
stituting the eigenfunctions (jl6.ip into the corresponding mode-sum formula, we obtain (no 
summation over i) 

= i^fcs e7 E " (16-16) 

where we have introduced the notations 

= (-l)^(2W + l)[r(x)+i^V'(x)/x2]+/2(a;), (16.17) 
f!)'Kk,fix)] = (-iyf'^ix)-[l + {-iyi.yx^]fix), (16.18) 

with i = 0, 3 and j = 1,2. As in the case of the field square, we apply summation formula 
(j3.18p to rewrite the sum over I. This enables us to present the VEV as the sum of boundary- free 
and boundary-induced parts as follows 

(0|7^*^|0) = (0,|7;'=|0,) + (7;'=)„. (16.19) 

The part induced by the cylindrical shell may be written in the form (no summation over i) [33] 

iTt)a = 1^ E' E r dxx3^^i^«[V(xr)], (16.20) 

^ n=0 A=0,1 [Xa) 

with the notations 

= 4'\fiy)] = fiy), (16.21) 
F!)'Hfiy)] = -(-i)7"(y) - [1 - (-i)VVy'] f{y), ^ = 1,2. (I6.22) 

As it can be easily checked, this tensor is traceless and satisfies the covariant continuity equation 
Vfc(Tf )a = 0. By using the inequalities [I„{y)K^{y)]' < and /^(y) < y^l + i^'^ /y^I,y{y), and the 
recurrence relations for the modified Bessel functions, it can be seen that the boundary-induced 
parts in the vacuum energy density and axial stress are negative, whereas the corresponding 
radial and azimuthal stresses are positive. The VEV of the EMT in the region outside the 
cylindrical shell is obtained from (116. 20p by the replacements I ^ K. 

The renormalized VEV of the EMT for the geometry without the cylindrical shell is obtained 
by using the corresponding formulae for the field square (116. lip . For the corresponding energy 
density one finds |130j 

Other components are found from the tracelessness condition and the continuity equation. 

Now, let us discuss the behavior of the boundary-induced part in the VEV of the EMT in 
the asymptotic region of the parameters. Near the cylindrical shell the main contribution comes 
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from large values of n. Thus, using the uniform asymptotic expansions for the modified Bessel 
functions for large values of the order, up to the leading order, we find 

These leading terms do not depend on the planar angle deficit in the cosmic string geometry. 
Near the cosmic string the main contribution comes from the mode n = and we have 

(7^0°)b « -(r/)b « -(T|)b ^^^^l^dxx' 

Therefore, differently from the VEV for the field square, the boundary-induced part in the 
vacuum EMT is finite on the string for all values of q. 

The behavior of the boundary-induced part in the VEV of the EMT in the asymptotic regions 
of the parameter q is investigated in a way analogous to that for the field square. For g <C 1, we 
replace the summation over n by the integration in accordance with (|16.15p and the VEV tends 
to a finite limiting value which does not depend on q. Note that as the spatial volume element 
is proportional to 1/q, in this limit the global quantities such as the integrated vacuum energy 
behave as 1/q. In the limit q ^ 1, the contribution of the modes with n ^ 1 is suppressed by 
the factor exp[— 2gnln(a/r)] and the main contribution comes from the n = mode with the 
behavior oc q. Though in this limit the vacuum densities are large, due to the factor 1/q in the 
spatial volume the corresponding global quantities tend to finite value. The VEVs for the EMT 
of the electromagnetic field induced by perfectly conducting cylindrical shell in the Minkowski 
spacetime considered in [39], are obtained from the formulae given in this section taking q = 1. 

The VEVs of the field square and the EMT for electromagnetic field inside a perfectly 
conducting wedge with the opening angle 0o and with a conducting cylindrical boundary of 
radius a can be evaluated by the way similar to that used in this section for the geometry of 
a cosmic string. The corresponding results are given in [37]. In the geometry of a wedge the 
eigenfunctions for the vector potential in the region inside the cylindrical shell are given by 
formula (jl6.ip with the replacements exp{iqn(j)) — > sm{qn(j)), n = 1, 2, . . ., for TM modes (A = 0) 
and exp{iqn4>) — > cos{qn(j)), n = 1,2,..., for TE modes (A = 1), where now q = 1^/(1)^. In 
the case of a wedge the VEVs of the field square and the EMT depend on the angle ^. The 
vacuum energy density induced by the cylindrical shell is negative for the interior region and 
the corresponding vacuum forces acting on the wedge sides are always attractive. 

17 Vacuum densities in the region between two coaxial cylin- 
drical surfaces 

17.1 Scalar field 

In this subsection, we consider the positive frequency Wightman function, the VEVs of the field 
square and the EMT for a massive scalar field with general curvature coupling parameter in the 
region between two coaxial cylindrical surfaces with radii a and b, a < b, on background of the 
(D + l)-dimensional Minkowski spacetime [42]. In an appropriately chosen cylindrical system 
of coordinates the corresponding line element has the form (jlS.ip with ^ (p ^ 2tt. We will 
assume that on the bounding surfaces the field obeys the boundary conditions 

(i,- + 4n;,)Vi) V (x) l^. = 0, j = a, b, (17.1) 

with Aj and Bj being constants, n^^.^ is the inward-pointing normal to the bounding surface 
r = j. For the region between the surfaces, a ^ r ^ b, one has •^ = njS\ with the notations 



Kojx) Ki{x) 
Io{x) Ii{x) 



-0.0168- 



;i6.25) 
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Ua = 1 and rih = —1. The eigenfunctions are specified by the set of quantum numbers a = 
(7, n, k), n = 0, ±1, ±2, . . ., and have the form 

(faix) = Pag\n\iia,jr) exp{in(j) + ikr^i -iut), (17.2) 

with 

<7„(7a,7r) = yi'^)(7a)J„(7r) - JiC\ia)Y^{jr), (17.3) 

and the other notations are the same as in Section [T5l The barred notation is defined by formula 
(|4.2p with the coefficients 

Aj = Aj, Bj = rijBj/j, j = a, b. (17.4) 

From the boundary condition on the surface r = 6 it follows that the possible values of 7 are 
solutions to the equation (77 = b/a) 

C^\v,ia) ^ jt\7a)Y^'\^b) - Yt\^a)J^\^b) = 0. (17.5) 

The corresponding positive roots we will denote by 7a = 7^^;, / = 1,2, . . ., assuming that they 
are arranged in the ascending order, ^ri,i < ln,i+i- 

From the orthonormality condition for the eigenfunctions, for the coefficient f3a- one finds 



4u;a(27r)^-i 



/3. = . \:.nZ ^ (17.6) 



with the notation from (j4.7p . Substituting eigenfunctions into the mode-sum formula (j9.8p . for 
the positive frequency Wightman function one finds 

°° " °° zgn{z,zr/a)gniz,zr'/a) 



■J ^ — n I 1 



2 



z) cos(nA(/)) exp(ikAr|| — itjAt) , (17.7) 



where, as before, k"^ = /c^ + and the prime on the summation sign means that the summand 
with n = should be halved. For the further evaluation of this VEV we apply to the sum over 
I summation formula (j4.14p with 

x) = exp(-iAf Va;7a + ^m)- (17.8) 

a/ X Kr^Cl 

The corresponding conditions are satisfied if r + r' + |At| < 2b. In particular, this is the case 
in the coincidence limit t = t' for the region under consideration. Now we can see that the 
application of formula (j4.14p allows to present the Wightman function in the form 



W{x,x') = ^ cos{nA(t>) / Ae 



;kArii 



1 r ^ h{z) 2 xnau{az,bz) 
dz , ,„ TTT^ / dz- 



X G^f ) {az, zr)G(f ) {az, zr') cosh(AVz2 _ ^2^) 



(17.' 



with the notations from (I12.30p . 
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In the limit 6 — > oo the second term in figure braces on the right of p7.9p vanishes, whereas 
the first term does not depend on b. It follows from here that the part with the first term 
presents the Wightman function in the region outside of a single cylindrical shell with radius a. 
As a result the Wightman function is presented in the form 



W{x,x') = WM{x,x') + {ip{xMx'))^ 

n=0 ''^"^ V ~ ^rn 



X 



Gl^\az, zr)G^^^ {az, zr') cosh(AtV^^ - ^^), (17.10) 



where Wm{x,x') is the Wightman function for a scalar field in the unbounded Minkowskian 
spacetime, and 



{^ix)ip{x'))^ = cos{nA^) d^ke'^^'\\ dz z 

^ n=0 

ll^\az) Kn{zr)Kn{zr') 



h{Aty^^^), (17.11) 



is the part of the Wightman function induced by a single cylindrical shell with radius a in the 
region r > a. Hence, the last term on the right of (jl7.10p is induced by the presence of the 
second shell with radius b. It can be seen that the Wightman function can also be presented 
in an equivalent form which is obtained from (jl7.10p by the replacement a — > 6 except in the 
argument of the function QanictZjbz). In this representation 

2^D °° , r roo 



{(fixMx'))^ = cos(nA0) / d^ke^'^^'li / dz z 



n=0 

Ki^\bz) In{zr)In{zr') i-^ /1^7 1o^ 

X -(fcN , / 9 ,9 cosh AtV^2 _ fc2^^ 1712 

is the part induced by a single cylindrical shell with radius b in the region r < b. This formula 
is also directly obtained from (jlS.lip taking (po = 2Tr. 

By making use of the formulae for the Wightman function and taking the coincidence limit 
of the arguments, for the VEV of the field square one finds 



n=0 



u 



x{u'^-m'^) 2 nj^{au,bu)Gl^'^^{ju,ru), (17.13) 
where j = a and j = b provide two equivalent representations and 



22-D 



Bd = D+._,^_^. - (17.14) 



■K—T 



For points away from the boundaries the last two terms on the right of formula ()17.13p are finite 
and, hence, the subtraction of the Minkowskian part without boundaries is sufficient to obtain 
the renormalized value for the VEV: (v?^)rcn = (Olv^^lO) — ((/7^)m. In formula (jl7.13p the part 
is induced by a single cylindrical surface with radius j when the second surface is absent. 



77 



The formulae for these terms are obtained from the Wightman function in the coincidence limit. 
For j = a one has 



oo 
n=0 



duuiu — m 



D-3 
2 



(17.15) 



and the formula for {(p^)b is obtained from here by the replacements a ^ b, I ^ K. The last 
term on the right of formula (jl7.13p is induced by the presence of the second cylindrical surface. 

The VEV for the EMT is obtained by using the formulae for the Wightman function and 
the VEV of the field square: 



, poo 

n=0 -^"^ 



U 



X [u 



m2) 2 njn{au,bu)Fi'^G^J\ju,ru)], 
where the notations Fi^\f(z)] are defined by formulae (jl5.2ip - (jl5.23p with f{z) 



G 



(17.16) 

(^■)(jn,z) 



and Fr[f{z)] = Fr[f{z)] for i = 3, . . . , Z) - 1. In formula (fT7T6D . the term {T^)^^^ is induced 
by a single cylindrical surface with radius j. These parts for both interior and exterior regions 
are investigated in [H] and are obtained from formulae of subsection 115.21 taking = 2tx. The 

,277(0 



formula for the case j = a is obtained from (jl7.15p by the replacement (ru) 



FA^^[i^„(™)]. 



The vacuum force per unit surface of the cylinder at r = j is determined by the \ - component 
of the vacuum EMT at this point. Similar to the case of spherical geometry, for the region 
between two surfaces the corresponding effective pressures can be presented as the sum of self- 
action and interaction parts, formula (jl2.38p . The interaction forces are determined from the 
last term on the right of formula (|17.16p with i = k = \ taking r = j: 



^(int) 



Ad 
2r 



n=0 



/ du u (yV? — rn?^ ^ Qjn{au,bu) 



[(nVj' + B] + A^njA,B,/j - A]] . 



(17.17) 



The expression on the right of this formula is finite for all non-zero distances between the shells. 
It can be seen that the vacuum effective pressures are negative for both Dirichlet and Neumann 
scalars and, hence, the corresponding interaction forces are attractive. For the general Robin case 
the interaction force can be either attractive or repulsive in dependence on the coefficients in the 
boundary conditions. By using the properties of the modified Bessel functions, the interaction 
forces per unit surface can also be presented in another equivalent form j42j : 



„(i) 



. 00 

Aorij ^/ 



2j 



E 

n=0 



duu (u — m 



^^ 2 



1+ 4e-lMT^ 



-In 

dj 



liC\au)Kr{bu 



(b), 



(17.18) 



The relation between the interaction forces and the corresponding bulk and surface Casimir 
energies in the geometry under consideration is discussed in [42]. 



17.2 Electromagnetic field 

In this subsection we consider the VEV for the EMT of the electromagnetic field in the region 
between two coaxial cylindrical surfaces with radii a and b, a < b [lO]. The corresponding 



78 



eigenfunctions satisfying the boundary condition on the surface r = a have the form 

A =R / (lA^) (7^63 + iA;Vt) PAn(7«,7'^)exp [z(n</) + /cz - wt)] , A = ^^^ ^^\ 

" ^"l -eax Vt{PAn(7«,7Oexp[i(n0 + A:z-a;t)]}, A=l ' ^"■'''> 

where 

Pxn{x,y) = J„(y)i;W(x) - (17.20) 

and the other notations are the same as in (jl6.ip . Prom the boundary conditions on r = 6 one 
obtains that the eigenvalues for 7 have to be solutions to the following equations 

d^Pxn{ia,-ir)\r=b = f). (17.21) 

These equations have an infinite number of simple real solutions. The eigenvalue equation 
(117. 2ip can be written in terms of function (14. ip as 

<^f(^,7&) =0, Aa = Ab = l-\,Ba = Bi, = \ A = 0,1 (17.22) 

(see notation (13. ip ). Now the normalization coefficient (3a is given by the formulae 



0= — 



Jl{z)/Jl{zii)-l] \ A = 



(1 - nyz^r]^) J'^\z)/J'„^{zr]) - 1 + n^/z^ 



-1 , (17.23) 

A = 1 



where z = 'ja, r] = b/a. Note that this coefficients can be expressed in terms of function (14. 7p as 

Pl = ^^Tf{r,,z). (17.24) 

Using these relations and introducing the cutoff function ipfj^, the vacuum EMT can be written 
in the form of the following finite integrosums (no summation over i) 



n=o ■'-'^ A=o,i 1=1 ^-fl^)^ + k^a^ 
xr^^^,7S)7S^'"'''/«[fe,PAn(7ff,72^/a)], (17.25) 

where 7a = J^^i , I = 1,2,..., are the solutions to eigenvalue equations (|17.2ip and the ex- 
pressions for the functions fn \k, Px-nijajx)] are obtained from (jl6.17p and (|16.18p taking 
f{x) = PxnilCLyX). By choosing in formula (j4.14p 

h{z) = "^^^^^^mhPxn{z,zx)], (17.26) 

one obtains that the VEV of the electromagnetic EMT in the region between two coaxial con- 
ducting cylindrical surfaces is presented in the form [40] 

(0|i;'=|0) = (0|i;'=|0)('') + {T^)ab, a<r<b. (17.27) 

In this formula the first term on the right is given by (no summation over i) 
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and for the second term we have 



{-^i dz 2^ z (17.29) 



with the notation Fi,'\f{y)] from (116.211) . (116. 22|) and 

Qxn{z,y) = Ki^\z)In{y) - li^Hz)Kn{y). (17.30) 

Note that the part (|17.29p is finite for a ^ r < b and in this part we have removed the cutoff 
function. 

In the hmit b ^ oo the second term on the right of (jl7.27p tends to zero, whereas the first 
one does not depend on b. From here we conclude that the term (0|Tj''|0)^'^^ is the VEV in the 
region outside a single cyhndrical shell with radius a. For the renormalization of this term we 
subtract the contribution of unbounded Minkowski spacetime which can be presented in the 
form: 

{0\Tt\0)M = j^Y! dk dz 'l^P_ fi^^[k,Urz)]. (17.31) 

„_Q OO Jo V K + Z 

By using the identity 

Ju [az) + Yi) [az) ^ ^=1^2 [az) 

and rotating the integration contour for z by angle 7r/2 for s = 1 and by angle —it/2 for s = 2, 
for the renomalized components we obtain 



n=0 



In(az) ^ r„[az) 



Kn{az) K^{az) 



F«[i^„(rz)]. (17.33) 



Note that the corresponding expressions in the region inside a single cylindrical shell are obtained 
from (jl6.20p taking q = 1 and differ from (|17.33p by the replacements I ^ K. As in the case of 
the interior components, the renormalized VEV (|17.33p is divergent when r ^ a with the leading 
terms given by formulae (jl6.24p . The corresponding asymptotic behavior at large distances from 
the cylinder, r 3> a, can be found from (117. 33p introducing a new integration variable y = rz 
and expanding the integrands over a/r. In this limit the main contribution comes from the 
lowest order mode with n = and one obtains (no summation over i) 

Here, compared to the spherical case, the corresponding quantities tend to zero more slowly. 
From the continuity equation for the vacuum EMT one has the following integral relation 

{Tl)a = -J drr{T^U = ^^, (17.35) 

where -E'°y*(r) is the total energy (per unit length) outside the cylindrical surface of radius r. 
Combining this relation with the similar relation in the interior region, for the total vacuum 
energy of the cylindrical shell per unit length we obtain 

E,yi = 4^i(a) + E^^^[a) = vra^ [{Tl)a{a+) - (ri),(a-)] . (17.36) 
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By taking into account the corresponding expressions for the radial component this yields 



In{z)Kn{z) 



_1 I /■CO T 

V / dzx^{z/a)z''-\n I - z\ln{z)Kn{z)f . (17.37) 
^ -'0 dz I J 



47ra2 

n=0 

In the last expression, integrating by part and omitting the boundary term we obtain the Casimir 
energy in the form used in numerical calculations. The corresponding results are presented in 
|1231 1124[ I125j . Note that in the evaluation of the Casimir energy for a perfectly conducting 
cylindrical shell by the Green function method to perform the complex frequency rotation pro- 
cedure an additional cutoff function has to be introduced (see |123j ) . This is related to the above 
mentioned divergence of the integrals over z for r = a. The results of the numerical evaluations 
for the energy density and pressures distributions (formula (|17.33p ) are presented in [281 139j . 
The energy density, (rg)^, and azimuthal pressure, — (T|)(i, in the exterior region are positive, 
and the radial pressure, — (r/)^, is negative. The ratio of the energy density to the azimuthal 
pressure is a decreasing function on r, and 1/3 ^ ~ {"^o) a / {T2) a ^0.5. Note that this ratio is a 
continuous function for all r and monotonically decreases from 1 at the cylinder axis to 1/3 at 
infinity. 

The quantities (117.271) with (I17.29P and (117.331) represent the renormalized VEV of the EMT 
in the region between two coaxial conducting cylindrical surfaces. Let us consider the limiting 
cases of the term (117. 29|) . First let a/r, a/h <C 1. After replacing z ^ bz and expanding the 
integrand over a/r and a/6 it can be seen that 

{T^)ab « {T^)b, a/r, a/h « 1, r < 6, (17.38) 

where (Zf )& is the vacuum EMT inside a single cylindrical shell with radius h. When a ^ b the 
sum over n in (jl7.29p diverges. Consequently for b — a <^h the main contribution comes from 
large n. By using the uniform asymptotic expansions for the modified Bessel functions, in this 
limit one obtains the corresponding quantity for the Casimir parallel plate configuration. 

From (jl7.27p and ()17.29p it can be seen that the VEV of the EMT can also be written in 
the form 

{Q\T^\Q) = {Q\Tf\^)^^) + {T^),a, a<r<b, (17.39) 



where 



/^k. r V ..3 F!,'\Q,^{bz,rz)]e\az)/li'\bz) 

n=o Jo ;^^o^ K;^ ' {az)Iii ' {bz) - Ki^ ' {bz)I^i ' (az) 



The quantities (|17.40p are finite for all a < r ^ 6 and diverge on the surface r = a. 

Now we turn to the interaction forces between the cylindrical surfaces due to the vacuum 
fluctuations. The vacuum force acting per unit surface at r = j is determined by the }- 
component of the EMT evaluated at this point. The corresponding effective pressures are 
presented as a sum of the self-action and interaction parts, similar to ()17.18p . Substituting into 
formulae (|17.29p and (|17.40p with i = k = 1 the values r = a and r = b, respectively, and 
using the Wronskian for the modified Bessel functions, for the interaction parts of the vacuum 
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pressures on the cylindrical surfaces one has [40]: 



^ At^i -^0 iz)In (bz/a) - K'rC> {bz/a)l)C' {z) 



4.2^4 h hi ^° Kl^\az/b)ll^\z) - K^\z)I^^\az/b) 

In p7.4ip and p7.42p the summands with A = (A = 1) come from the electric (magnetic) waves 
contribution. By using the properties of the modified Bessel functions it can be seen that the 
quantities p)- \^ are negative and, hence the corresponding interaction forces are attractive. Note 
that interaction forces p7.4ip . (jl7.42p can also be obtained by differentiating the corresponding 
part in the vacuum energy with respect to the radii |126j . 



18 Polarization of the Fulling- Rindler vacuum by a uniformly 
accelerated mirror 

It is well known that the uniqueness of the vacuum state is lost when we work within the 
framework of quantum field theory in a general curved spacetime or in non-inertial frames. 
In particular, the use of general coordinate transformations in quantum field theory in flat 
spacetime leads to an infinite number of unitary inequivalent representations of the commutation 
relations. Different inequivalent representations will in general give rise to different vacuum 
states. For instance, the vacuum state for a uniformly accelerated observer, the Fulling-Rindler 
vacuum, turns out to be inequivalent to that for an inertial observer, the familiar Minkowski 
vacuum. An interesting topic in the investigations of the Casimir effect is the dependence of the 
vacuum characteristics on the type of the vacuum. In this section, by using the GAPF, we will 
study the scalar vacuum polarization brought about by the presence of infinite plane boundary 
moving by uniform acceleration through the Fulling-Rindler vacuum. The corresponding VEVs 
of the EMT were studied by Candelas and Deutsch |131j for conformally coupled 4D Dirichlet 
and Neumann massless scalar and electromagnetic fields. In this paper only the region of the 
right Rindler wedge to the right of the barrier is considered. In (48j we have investigated the 
Wightman function, the VEVs of the field square and the EMT for a massive scalar field with 
general curvature coupling parameter, satisfying Robin boundary condition on infinite plate 
in an arbitrary number of spacetime dimensions and for the electromagnetic field. The both 
regions, including the one between the barrier and Rindler horizon are considered. 



18.1 Wightman function 

In the accelerated reference frame it is convenient to introduce Rindler coordinates (r, ^,x) 
related to the Minkowski ones, (t, x^,x) by the formulae t = ^sinhr, = ^coshr, and 
X = (x^, . . . ,x^) denotes the set of coordinates parallel to the plate. In these coordinates the 
Minkowski line element takes the form 

ds^ = fdr^ - df - dx^, (18.1) 

and a world-line defined by ^, x = const describes an observer with constant proper acceleration 
The Rindler time coordinate r is proportional to the proper time along a family of uniformly 
accelerated trajectories which fill the Rindler wedge, with the proportionality constant equal to 
the acceleration. 
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We will assume that the plate is located at = a, with being the proper acceleration, 
and the field satisfies boundary condition (j9.1ip . We will consider the region on the right from 
the boundary, ^ ^ a. In Rindler coordinates the boundary condition takes the form 

{A + Bd^)(p = 0, i = a. (18.2) 

In the region ^ > a, a, complete set of solutions that are of positive frequency with respect to 
d/dr and bounded as — > oo is 



ip^{x) = [3„KUKy''''-'^\ A = v/fc' + m2, a = (cu,k). (18.3) 

From boundary condition p8.2|) we find that the possible values for w have to be zeros of the 
function Ki^{\a)^ where the barred notation is defined by (jS.ip with the coefficients 

A = A, B = B/a. (18.4) 

We will denote these zeros by a; = w„ = ^n{k), n = 1,2,..., arranged in ascending order: 
LOn < t^ra+i- The coefficient /3cr in (jl8.3p is determined by the normalization condition with 
respect to the standard Klein-Gordon inner product, with the ^-integration over the region 
(a,oo). From this condition one finds 

^'^ (2^)^-1 5^i^,^(Aa)U=^„- ^ ■ ' 

Substituting the eigenfunctions (jl8.3p into (j9.8p . for the Wightman function we obtain 

/ ikAx °° T- (\n) 
''kp^E3j^/^,UAf)A-„(A5>— 'U„., (18.6) 

with Ax = X — x'. At = t — t' . For the further evaluation of VEV (jl8.6p we apply to the sum 
over n summation formula (|5.12p . As a function F{z) in this formula we choose 

Fiz) = Ki,[Xe,)Ki,{\Oe-"''\ (18.7) 

Using the asymptotic formulae for the modified Bessel functions it can be seen that condition 
(15. 5p is satisfied if a'^e'^'^l < In particular, this is the case in the coincidence limit r = r' for 
the points in the region under consideration, > a. With F{z) from (jl8.7p . the contribution 
corresponding to the integral term on the left of formula (j5.12p is the Wightman function for 
the Fulling-Rindler vacuum without boundaries: 

Wk{x,x') = -^ j dkj^-yj^ duJsm\l{^TUJ)e-'^^^KUK)K^^{K')■ (18.8) 

Taking into account this and applying summation formula ()5.12p to Eq. ()18.6p . we receive |48) 

W{x, x') = VFr(x, x') + {ip{x)^{x'))a, (18.9) 
where the second term on the right is induced by the barrier: 

1 f ikAx poo J 

{v{xMx'))a = — / dk— _ / du;^^KUX()KUXacosHu;AT), (18.10) 

and is finite for ^ > a. The divergences in the coincidence limit are contained in the first term 
corresponding to the Fulling-Rindler vacuum without boundaries. The boundary-induced part 
in the Wightman function in the region < < a is obtained from (jl8.10p by the replacements 
I^K I48j (see below). 
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18.2 Vacuum expectation values of the field square and the energy-momentum 
tensor 

In the coincidence limit, from formula (llS.lOp for the boundary part of the field square we have 



dk k 



duj 



18.11 



This quantity is monotone increasing negative function on ^ for Dirichlet scalar and monotone 
decreasing positive function for Neumann scalar. Substituting the function (|18.9p into Eq. (|9.10p 
and taking into account Eqs. ()18.10p and (|18.1ip . for the VEV of the EMT in the region ^> a 
one finds 



:i8.i2) 



where the first term on the right is the VEV for the Fulling-Rindler vacuum without boundaries, 



{mm. 



7r(D+3)/2p 



;i8.13) 



and the second term is the contribution brought by the presence of the barrier: 



l2-Drfc 



dkk'^-^X^ 







vr(^+i)/2r (^) 70 
In formula (jl8.13p we have introduced the notations 



K^{\a) 



;i8.14) 



2 ^ 



dg{z) 



dz 



z dz 



1 



- - 2C + - + 2C 



21 



\9i^)\'' 



dg{z) 



dz 
2C 



dg{z) 



dz 



+ 1 



i5(^)r 

A2 



(18.15) 



m 



{D - 1)A2 



15(^)1% 



with i = 2,3, ■■■,D, and the expressions for the functions F^'^\g{z)\ are obtained by the replace- 
ment oj ^ ioj: 

F'^'^[g{z)]=f'^%{z),u:^iuj]. (18.16) 

The boundary-induced parts in the VEVs of the field square and the EMT for the region ^ < a 
are obtained from (jlS.lip and (118. 14p with the replacements I K. It can be easily checked 
that both summands on the right of Eq. (|18.12p satisfy the continuity equation VkT^ = 0, 
which for the geometry under consideration takes the form d^{^T^) — Tq = 0. 

The purely Fulling-Rindler part (jl8.13p of the EMT is investigated in a large number of 
papers (see, for instance, references given in [49]). The most general case of a massive scalar 
field in an arbitrary number of spacetime dimensions has been considered in Ref. |132j for 
conformally and minimally coupled cases and in Ref. [4^ for general values of the curvature 
coupling parameter. For a massless scalar the VEV for the Rindler part without boundaries can 
be presented in the form 



(T^Vsub = (0|7;^|0)r-(0|7;^|0)m 

2^fr^"' 
(47r)^/2r(D/2) 



duj 



^277 a; 



+ (-1) 



D 



(18.17) 
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where the expressions for the functions g^^\uj) are presented in Ref. [35], and (0|r.'^|0)M is 
the EMT for Minkowski vacuum without boundaries. Expression (jl8.17p corresponds to the 
absence from the vacuum of thermal distribution with standard temperature T = (27r^)~^. In 
general, the corresponding spectrum has non-Planckian form: the density of states factor is 
not proportional to uj^^^duj. The spectrum takes the Planckian form for conformally coupled 
scalars in = 1,2,3 with g'^^\ijj) = —Dg'^'^\uj) = 1, z = 1, 2, . . . D. It is of interest to note that 
for even values of spatial dimension the distribution is Fermi-Dirac type (see also j l33j ). For a 
massive scalar the energy spectrum is not strictly thermal and the corresponding quantities do 
not coincide with ones for the Minkowski thermal bath. 

Boundary part (|18.14p is finite for all values ^ > a and diverges at the plate surface ^ = a. 
To extract the leading part of the boundary divergence, note that near the boundary the main 
contribution into the w-integral comes from large values of lo and we can use the uniform asymp- 
totic expansions for the modified Bessel functions [65]. Introducing a new integration variable 
k Lok and replacing the modified Bessel functions by their uniform asymptotic expansions, it 
can be seen that the leading terms do not depend on the mass and the Robin coefficients and 
are the same as for the plate in Minkowski spacetime. 

Now let us consider the asymptotic behavior of the boundary part (jl8.14p for large ^, ^ » a. 
Introducing in Eq. (jl8.14l) a new integration variable 2/ = and using the asymptotic formulae 
for the modified Bessel functions for small values of the argument, we see that the integrand is 
proportional to . It follows from here that the main contribution into the cj-integral 

comes from small values of u). Expanding with respect to oj, in the leading order we obtain 

^f(-i)^.og-^-i4-)(^^) 

^ ^'"^ 2^37r(^-i)/2(l + 5^o)r(^)ln2(2e/a)' ^ ' 

where 

... roo 



If, in addition, one has ^ 1, the integral in this formula can be evaluated replacing the 
McDonald function by its asymptotic for large values of the argument. In the leading order this 
yields 

Af ^ Af ^ -2me4'^ ^ ^(1/4 - C)r (^^) (mO(''+'^/'e-2-«, (18.20) 

with i = 2,3,-- - ,Z). For massless case the integral in Eq. (jl8.19p may be evaluated and one 
obtains 

^(0) ^ ^ _^^(2) ^ 2^^(Cc - C) r4 (D + i\ 

For a conformally coupled scalar the leading terms vanish and the VEVs are proportional to 
^-■^-Mn~3(2^/a) for ^/a > 1. From Eq. (fTSlSj) we see that for a given ^ the boundary 
part tends to zero as a ^ (the barrier coincides with the right Rindler horizon) and the corre- 
sponding VEVs of the EMT are the same as for the Fulling- Rindler vacuum without boundaries. 
Hence, the barrier located at the Rindler horizon does not alter the vacuum EMT. 

And finally, we turn to the asymptotic a, ^— >oo, £, — a = const. In this limit ^/a ^ 1 and 
w-integrals in Eq. ()18.14p are dominated by large uj. Replacing the modified Bessel functions 
by their uniform asymptotic expansions, keeping the leading terms only and introducing a new 
integration variable v = cu/a, we can see that the VEVs coincide with the corresponding VEVs 
induced by a single plate in the Minkowski spacetime with Robin boundary condition on it. 
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18.3 Electromagnetic field 

We now turn to the case of the electromagnetic field in the region > a for the case D = 3. 
We will assume that the mirror is a perfect conductor with the standard boundary conditions 
of vanishing of the normal component of the magnetic field and the tangential components of 
the electric field, evaluated at the local inertial frame in which the conductor is instantaneously 
at rest. As it has been shown in |131j . the corresponding eigenfunctions for the vector potential 
Ai{x) may be resolved into the transverse electric (TE) and transverse magnetic (TM) (with 
respect to ^-direction) modes Aaa-h = ("^jk): 

Aoai{x) = (0,0,-iA:3,iA:2)v?o<x(x), a = 0, (18.22) 
Ai,i{x) = (-C6»5,icu/e,0,0)99i,(x), a = 1, (18.23) 

where the eigenfunctions </?o-o(x) are given by formula (jl8.3p with m = 0, k = {k2,k3), and 
a = 0, 1 correspond to the TE and TM waves, respectively. From the perfect conductor boundary 
conditions on the vector potential we obtain the corresponding boundary conditions for the scalar 
modes ipaaix)- 

ipo^ = 0, di:ipi^ = 0, e = a. (18.24) 

As a result the TE/TM modes correspond to Dirichlet/Neumann scalars. In the corresponding 
expressions for the eigenfunctions Aa^^i^x) the normalization coefficient is determined from the 
orthonormality condition 

dx / -^Al^Al,^,i = 5,„.<5../. (18.25) 

Ja i ^ 

The eigenvalues for w are the zeros of the function Ki^(ka), k = |k|, for the TE modes and the 
zeros of the function K-^(ka) for the TM modes. 

Substituting the eigenfunctions for Aaaiix) into mode-sum formula (114. Th with the standard 
bilinear form for the electromagnetic field EMT and applying to the sums over u;„ formula 
we find 

i:k poo roc 



where 



(0|I;'=|0)r = (0|T;'=|0)m-^ du;-2^;^diag(l, --,--,--) (18.27) 

is the VEV for the Fulling-Rindler vacuum without boundaries |131] . In formula ()18.26p the 
notations 

F^il ioiz)] = {-iyg'\z) + (1 - {-iyu;'/z')g\z), i = 0, 1, (18.28) 
F^^igiz)] = -g\z), i = 2,3, (18.29) 

are introduced. Formula (|18.26p is derived in |131] by using the Green function method. Note 
that the w-integral in (118. 271) is equal to 11/240. The VEV for the EMT of the electromagnetic 
field in the region ^ < a is obtained from (jl8.26p by the replacements ^ K^^ [38]. By using 
the properties of the modified Bessel functions it can be seen that 

lujika) ^ lL{ka) ^ ^ K^jka) ^ K'^jka) ^ ^ 



K^{ka) Ki{ka) I^ka) I^ka) 
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and F}^ [K^z)] > 0, F}^ [lUz)] > 0, FeW [K^z]] < 0, FeW [/^(z)] > 0. In particular, from 
these inequalities it follows that the energy density induced by the plate is positive in the region 
^ < a and negative in the region ^ > a. The corresponding graphs are given in [48\ I131j . The 
behavior of the plate-induced VEVs in various asymptotic regions of the parameters can be found 
in [481 1131j . In Ref . |134j the formulae derived in this section were used to generate the vacuum 
densities for a conformally coupled massless scalar field in de Sitter spacetime in presence of a 
curved brane on which the field obeys the Robin boundary conditions with coordinate dependent 
coefficients. 



19 Vacuum densities in the region between two plates uniformly 
accelerated through the Fulling-Rindler vacuum 

19.1 Wightman function 

Consider two parallel plates moving with uniform proper accelerations assuming that the quan- 
tum scalar field is prepared in the Fulling-Rindler vacuum [49 ^ I50j . We will let the surfaces 
^ = a and ^ = 6, 5 > a > 0, represent the trajectories of the plates, which therefore have proper 
accelerations a^^ and b~^. We consider the case of a scalar field satisfying Robin boundary 
conditions on the surfaces of the plates: 

= 0, j = a,b, (19.1) 

with constant coefficients Aj and Bj. The plates divide the right Rindler wedge into three 
regions: < ^ < a, ^ > 6, and a < ^ < 6. The VEVs in two first regions are the same as those 
induced by single plates located at ^ = o and ^ = 6, respectively. In the region between the 
plates the eigenfunctions satisfying boundary condition on the plate ^ = 6 have the form 

^Ax) = /3.gS2(A6, Ae)e^'^™^ (19.2) 

with the function Gu\x,y) defined by formula ()12.30p . Note that the function {Xb, X^) is 
real, G^^-^{Xb, X^) = Gf^{Xb, X^). From the boundary condition on the plate ^ = a we find that 
the possible values for oo are roots to the equation 

Zi^(Aa,A6) = 0, (19.3) 

where the function Zi^{u,v) is defined by formula (j5.14p . and the barred notations are defined 
by formula (|4.2p with Aj = Aj, Bj = Bj/j. For a fixed A, the equation (119. 3p has an infinite 
set of real solutions with respect to u. We will denote them by 0„ = Qn{Xa, Xb), 0„ > 0, 
n = 1,2,..., and will assume that they are arranged in ascending order Q„ < In the 

consideration below we will assume the values of the coefficients in Eq. ()19.ip for which the 
imaginary solutions are absent and the vacuum is stable. 

The coefficient jSa in formula ()19.2p is determined from the orthonormality condition. Taking 
into account the boundary conditions, for this coefficient one finds 

(2^)^-^lg(Aa) 
" 4^(A6)c'^Zi^(Aa,A6) 



(19.4) 
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Now, substituting the eigenfunctions into mode-sum formula (|9.8p . for the positive frequency 
Wightman function one finds 

w{x,x') = / ^%^y J'-^^^^ 

xGfJ{Xb,XC)GfJ{Xb,Xa ^ . (19.5) 

For the further evaluation of this VEV we apply to the sum over n summation formula (jS.lSp . 
As a function F{z) in this formula we choose 

p, , _ GS(Ai».AaGS(Ai,Ar) -i.A. „„ f.> 

4"(A^/«(A^) • * ' 

Condition (|5.19p for this function is satisfied if a^el^"^' < By using formula (|5.18p . for the 
Wightman function one obtains the expression |50j 



W{x,x') = W^^\x,x')- -^j—jj— dujnbU^a,Xb) 



xG^^\Xb,X0G^^H^b,XOcosh{ujAT), (19.7) 



where we have used the notation (jl2.35p . In Eq. (jl9.7p . 

.ArGfJiXb,XOGfJ{Xb,Xe) 



W^''\x,x') = I / da;sinh(vr^) 



{191 



is the Wightman function in the region ^ < b for a single plate at ^ = 5. This function is 
investigated in Ref. [48] and can be presented in the form 

W^^\x,x') = Wr{x,x') + {ipixMx'))^, (19.9) 

where Wji{x,x') is the Wightman function for the right Rindler wedge without boundaries and 
the part 

f dke^^^^ f°° K^\Xb) 
{^{x)^{x'))^ = - / / du 



^(27r)«-i 7o ll^\xb) 
x/..(Ae)/^(Ae')cosh(u;AT) (19.10) 

is induced in the region ^ < 6 by the presence of the plate at ^ = 6. Note that the representation 
(fTOI) with (119.101) is vahd under the assumption < ^^glArl^ 
By using the identity 

£\Xb) K}^\Xa) 

+ nm,U)^a,Xb)G^j\Xj,XOz!j\Xj,Xa, (19.11) 

j=a,b 
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with n^") = 1, n^'') = — 1, the Wightman function can also be presented in the form 



xGi^\Xa, XOGi^\Xa, A^') cosh(c^AT). (19.12) 



W{x,x') = W^^\x,x')- I „^ / cLunaU>^a,Xb) 



In this formula W^"'\x, x') is the Wightman function in the region ^ > a for a single plate at 
^ = a, and is investigated in the previous section. In the coincidence limit the second term on 
the right of formula (|19.7p is finite on the plate £, = b and diverges on the plate at ^ = a, whereas 
the second term on the right of Eq. (jl9.12p is finite on the plate ^ = a and is divergent for 
C = b. Consequently, the form (119. Th ( (jl9.12p ) is convenient for the investigations of the VEVs 
near the plate ^ = b = a). 

19.2 Scalar Casimir densities 

Here we will consider the VEVs of the field square and the EMT in the region between the 
plates |50] . In the coincidence limit from the formulae for the Wightman function one obtains 
two equivalent forms for the VEV of the field square, 



-Bd / dkk^-^ / dionj^{Xa, Xb)G^J^'^{Xj, X^, (19.13) 
Jo Jo 

corresponding to j = a and j = b, and Bd is defined by (jl7.14p . In Eq. (jl9.13p the part (y^)^ 
is induced by a single plate at ^ = j when the second plate is absent. The last term on the 
right of formula (jl9.13p is finite on the plate at ^ = j and diverges for points on the other plate. 
Extracting the contribution from this plate, we can write the expression p9.13p for the VEV in 
the symmetric form 



j=a,b 



with the interference part 



poc poo 

A ((^2^ = -Bd / dkk^-"^ / dujfj\Xa) 
Jo Jo 

G^^\xb,XO KliXj) 
I^^\Xb)G%\Xa,Xb) K!f\Xa) 



(19.15) 



An equivalent form for this part is obtained with the replacements a b and I ^ K in the 
integrand. Interference term (jl9.15p is finite for all values of ^ in the range a ^ ^ ^ 6, including 
the points on the boundaries. The surface divergences are contained in the single plate parts 
only. In the limit a ^ b with fixed values of the coefficients in the boundary conditions, the 
interference part (jl9.15p is divergent and for small values of 6/a — 1 the main contribution comes 
from large values of to. Introducing an integration variable x = k/u) and replacing the modified 
Bessel functions by their uniform asymptotic expansions, we can see that to the leading order 
A((/?2^ coincides with the VEV of the field square in the region between two parallel plates in 
the Minkowski bulk. 

Large values of the proper accelerations for the plates correspond to the limit a, 6 — > 0. In 
this limit the plates are close to the Rindler horizon. From formulae (jlS.lip . (|19.15p we see that 
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for fixed values of the ratios a/h, ^/b, both single plate and interference parts behave as b^~^ in 
the hmit 6 — > 0. In the hmit a ^ for fixed values ^ and b, the left plate tends to the Rindler 
horizon for a fixed world-hne of the right plate. The main contribution into the w-integral in 
Eq. (jl9.15p comes from small values to, to < l/ln(2/Aa). Using the formulae for the modified 
Bessel functions for small arguments, it can be seen that interference part (jl9.15p vanishes as 
ln-2(26/a). 

In the limit of small accelerations of the plates, a,b ^ oo, with fixed values b — a, Bj/Aj, and 
m, the main contribution comes from large values lo. Using the uniform asymptotic formulae 
for the modified Bessel functions, it can be seen that a-nd A(9?^) coincide with the cor- 

responding expressions for the geometry of two parallel plates on the Minkowski bulk. In this 
limit, corresponds to the Cartesian coordinate perpendicular to the plates which are located 
at ^ = a and ^ = b. 

For large values of the mass, ma ^ 1, we introduce a new integration variable y = X/m. The 
main contribution into the w-integral comes from the values uj ~ ^/ma. By using the uniform 
asymptotic expansions for the modified Bessel functions for large values of the order and further 
expanding over uj/ma, for the single plate parts to the leading order one finds 



« - , n , „ ^ [n , j = «, b, (19.16) 
2(47r)-c,(m)|e-j|- 



where we have introduced notations 



/ , (19.17) 



Aj + nij^Bjy 

In the similar way, for the interference part we obtain the formula 



A{ip^) « ^ ^. (19.18) 

(Air) 2 ^Ca{m)cb{m){b - a) 2 

As we could expect, the both single plate and interference parts are exponentially suppressed 
for large values of the mass. 

Making use of the formulae for the Wightman function and the field square, one obtains two 
equivalent forms for the VEV of the EMT, corresponding to j = a and j = b (no summation 
over i) [50] (see also |49] for the case of Dirichlet and Neumann scalars): 

(OII^'^IO) = {0\Tt\0)K + {Tt),-BnS'l Jdkk^-' 

POO |- _ 

xA^ / (ia;%^(Aa,A6)F« G^J\Xj,XO ■ (19.19) 
Jo L -I 

In this formula, (0|7;'=|0)r is the corresponding VEV for the Fulling-Rindler vacuum without 
boundaries, and the term {T^)j is induced by the presence of a single plane boundary located 
at ^ = j in the region ^ > a for the case j = a and in the region ^ < 6 for j = b. In formulae 
()19.19p . for a given function g(z) we use the notations F^^^g{z)] from (jl8.16p . For the last term 
on the right of Eq. (jl9.19p we have to substitute g{z) = G[i\xj,z). It can be easily seen that 
for a conformally coupled massless scalar the EMT is traceless. 
Now let us present the VEV (119. 19p in the form 

(OII^'^IO) = (0|I;'=|0)r + {Thj + A(7;^), (19.20) 

i=a,b 
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where (no summation over i) 



dk k^-'^X'^ 



dujl^^\\a) 



F«[G(t)(A5,Ae)] FW[K^(Ae)] 



(a) I 



(19.21) 



is the interference term. The surface divergences are contained in the single boundary parts and 
this term is finite for all values a ^ ^ ^ 6. An equivalent formula for A(Tj^) is obtained from 
Eq. p9.2ip by replacements a b, I ^ K . Both single plate and interference parts separately 
satisfy the continuity equation. For a conformally coupled massless scalar field they are traceless 
and we have an additional relation (T?^ = 0. 

In the limit a — > 6 expression ()19.2ip is divergent and for small values of b/a — 1 the main 
contribution comes from large values of uj. Introducing a new integration variable x = k/tu and 
replacing the modified Bessel functions by their uniform asymptotic expansions for large values 
of the order, at the leading order one receives 



(4vr) 



" 2 



r(# + i) 



dy 



y 



D 



kahe^y^^-'') 



1 + 2D (1 - 5\) (C - Cd) 



:i9.22) 



where kj = 1 — 25oBj - In the limit of large proper accelerations for the plates, a,b ^ 0, for fixed 
values a/b and (,/b, the world-lines of both plates are close to the Rindler horizon. In this case 
the single plate and interference parts grow as b"^^^. The situation is essentially different when 
the world-line of the left plate tends to the Rindler horizon, a — > 0, whereas b and ^ are fixed. 
In the way similar to that for the case of the field square, it can be seen that in this limit the 
interference part (jl9.2ip vanishes as ln~^(26/a). 

In the limit of small proper accelerations, a, 6 — > oo with fixed values b — a, Bj/Aj, and 
m, the main contribution comes from large values of uj. Using the asymptotic formulae for the 
modified Bessel functions, to the leading order one obtains the corresponding expressions for the 
geometry of two parallel plates on the Minkowski background investigated in |32j . In particular, 
in this limit the single boundary terms vanish for a conformally coupled massless scalar. 

For large values of the mass, ma 3> 1, by the method similar to that used in the previous 
subsection for the field square, it can be seen that the both single plate and interference parts 



are exponentially suppressed (no summation over i): {T- 



m 



D/2+1 



exp[-2m|^ j = a,b, 



for single plate parts and A(T?) ~ m^^'^^^ exp[2m(a — b)] for the interference part. 

Now we turn to the interaction forces between the plates due to the vacuum fluctuations. The 
vacuum force acting per unit surface of the plate at ^ = j is determined by the {-component 
of the vacuum EMT evaluated at this point. The corresponding effective pressures can be 
presented as a sum of self-action and interactions terms, p^-'^ = Pi'^ +?'(int)' j ~ ^' ^' '^^^ ^^^^ 
term on the right is the pressure for a single plate at ^ = j when the second plate is absent. 
The first term on the right is the pressure for a single plate at ^ = j when the second plate 
is absent. This term is divergent due to the surface divergences in the subtracted VEVs and 
needs additional renormalization. This can be done, for example, by applying the generalized 
zeta function technique to the corresponding mode-sum. This procedure is similar to that used 
in |135j for the evaluation of the total Casimir energy in the cases of Dirichlet and Neumann 
boundary conditions and in |136j for the evaluation of the surface energy for a single Robin 
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plate. The term P(j4) is the pressure induced by the presence of the second plate. This term 
is finite for ah nonzero distances between the plates and is not affected by the renormalization 
procedure. For the plate at ^ = j the interaction term is due to the third summand on the right 
of Eq. ()19.19p . Substituting into this term i = j and using the Wronskian for the modified 
Bessel functions one finds [501 



(i) 

^(int) 



ApA] 
2p 



dk k 



D-2 



1 



Qj^{Xa, A5), 



;i9.23) 



with Pj = Bj/{jAj). In dependence of the values for the coefficients in the boundary conditions, 
the effective pressures p9.23p can be either positive or negative, leading to repulsive or attractive 
forces. It can be seen that for Dirichlet boundary condition on one plate and Neumann boundary 
condition on the other one has P/j^) > the interaction forces are repulsive for all distances 
between the plates. Note that for Dirichlet or Neumann boundary conditions on both plates 
the interaction forces are always attractive [49j. The interaction forces can also be written in 
another equivalent form 



n 



A 



— I dkk 
2j 



D-2 



1 + 



(4C - 1) (3j 



(AV+^2)/3? + /3,-l 



I^^\\a)K!^\\b) 



[19.24) 



For Dirichlet and Neumann scalars the second term in the square brackets is zero. 

The results obtained above can be applied to the geometry of two parallel plates near the 
D = 3 'Rindler wall.' This wall is described by the static plane-symmetric distribution of the 
matter with the diagonal EMT = diag(em5 —Pm, —Pm, —Pm) (see Ref. |137j ). Below we will 
denote by x the coordinate perpendicular to the wall and will assume that the plane a; = is 
at the center of the wall. If the plane the boundary of the wall, when the external 

{x > Xg) line element with the time coordinate t can be transformed into form (jlS.ip with 



e(x) 



Xs + 



1 



27r(7, 



2TTasy/goo{xs)t. 



;i9.25) 



In this formula the parameter ag is the mass per unit surface of the wall and is determined by 
the distribution of the matter: 



+ 3pm) [g{x)/g{xs)]^^'^ dx. 



;i9.26) 



For the 'Rindler wall' one has 522(^)U=o < |137] (the external solution for the case 522(^)U=o > 
is described by the standard Taub metric). Hence, the Wightman function, the VEVs for the 
field square and the EMT in the region between two plates located at x = and x = Xb, 



> Xs near the 'Rindler wall' are obtained from the results given above substituting j = ^{x 



j = a,b and ^ = (,{x). For ag > 0, x ^ Xg one has ^(x) ^ ^(xs 
regular everywhere in the external region. 



> and the Rindler metric is 



19.3 Electromagnetic field 

As in the case of a single plate considered in subsection 118. 3|, in the region between two uni- 
formly accelerated perfectly conducting plates the eigenfunctions for the electromagnetic vector 
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potential are resolved int TE and TM modes in accordance with (|18.22p . (jl8.23|) . with Dirichlet 
and Neumann boundary conditions: 

^Oa\^=a = fOcr\^=b = 0, d^(fir^\^=a = d^(fia\^=b = 0. (19.27) 

The corresponding eigenvalues for uj are the zeros of the function 

ZaMka, kb) = lt\kh)Kt\ka) - Kt\kh)lt\ka), (19.28) 

where a = 0, 1 for the TE and TM modes respectively, and fi^\x) = fiuj{x), fi^\x) = fl^{x) 
for f = I,K. The eigenfunctions for the separate scalar modes (paaix) are given by formula 
()19.2p with m = 0, D = 3 (for the case of an arbitrary D see |135j ) and with the replacement 

gS(A6, AC) ^ Gc,,Ukt kb), a = 0,1, where 

GaA^,y) = li''\y)KUx) - Ki''\y)Ux). (19.29) 

The corresponding normalization coefficients are determined from (jl8.25p where now the ^- 
integration goes over the interval {a,b). 

Substituting the normalized eigenfunctions into the mode-sum formula for the VEV of the 
EMT and applying formula (j5.18p for the summation over the eigenvalues of uj, similar to the 
case of a scalar field one finds ^49j 

(OIT^'^IO) = {Ttf^ - ^ Hdkk^ Hdco (19.30) 
4vr2yo Jo ^tH^i C\kb)Z^Aka,kb) 

where {T^)^^") is the VEV in the region ^ < b corresponding to the geometry of a single plate 
with ^ = b and the functions Fe^[g{z)] are defined by formulae (jl8.28p . (|18.29p with g{z) = 
GoL,uj{z, kb). An alternative form for the vacuum EMT in the region between two plates is 



47r^/n ./n K!f>{ka)Z^Aka,kb) 



where {T^)^"^ is the vacuum EMT in the region ^ > a corresponding to the geometry of a single 
boundary at = a. 

For the interaction force p)r'i;„.\, j = a,b, per unit area of the plate at f = j, from Eqs. 

t> 111 \ 1 liXi J 

p^30]) and (fira]) one obtains [49] 

"^^-(-^^ " ^-'a^h ''Vo K^\ka) Z^Aka.kb) ' ^^^'^^^ 

. = -Ai r^^^ r^- E (-1)"^^^^^^^^- (19-33) 



p. 



Recalling that {—!)" Zo,^^{ka, kb) > we see that the electromagnetic interaction forces are 

(i) (?) (?) (?) (7) 

attractive. Note that ^^^^(int) = PD{mt) +?'iv(int)' "^^^^^ PD{int) ^^'^ PN{mt) interaction 

forces for Dirichlet and Neumann scalars. In the limit a ^ 6, to the leading order of 1/(6 — a) 

from these expressions the electromagnetic Casimir interaction force between two parallel plates 

in the Minkowski spacetime is obtained. Note that the interaction forces (|19.32p . (jl9.33p can 

also be obtained by the differentiation the corresponding Casimir energy ^135j. 
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20 Wightman function and Casimir densities for branes on AdS 
bulk 



Anti-de Sitter (AdS) spacetime is one of the simplest and most interesting spacetimes allowed 
by general relativity. Quantum field theory in this background has been discussed by several 
authors. Much of early interest to AdS spacetime was motivated by the questions of principle 
nature related to the quantization of fields propagating on curved backgrounds. The importance 
of this theoretical work increased when it was realized that AdS spacetime emerges as a stable 
ground state solution in extended supergravity and Kaluza-Klein models and in string theories. 
The appearance of the AdS/CFT correspondence and braneworld models of Randall-Sundrum 
type [138] has revived interest in this subject considerably. Motivated by the problems of the 
radion stabilization and the generation of cosmological constant, the role of quantum effects in 
braneworlds has attracted great deal of attention (see, for instance, [57j for relevant references). 
In this section we apply the GAPF for the investigation of the positive frequency Wightman 
function and VEV of the EMT for a massive scalar field with general curvature coupling parame- 
ter subject to Robin boundary conditions on two parallel branes located on {D + l)-dimensional 
AdS background [56j . The general case of different Robin coefficients on separate branes is 
considered. 



20.1 Wightman function 

Consider a scalar field i^ix) on background of a (D + l)-dimensional plane-symmetric spacetime 
with the line element 

ds^ = Qikdx'dx^ = e-2'^(^)??^^dx^dx^ - dy"^ , (20.1) 

and with rj^jj = diag(l, —1, . . . , —1) being the metric for the L)-dimensional Minkowski space- 
time. Here and below in this section i. A; = 0, 1, . . . , D, and /i, z/ = 0, 1, . . . , D — 1. By making a 
coordinate transformation 

s'"(f)(i?/, (20.2) 



metric ()20.ip is written in a conformally-flat form ds^ = e~'^'^^y^ iji^dx^ dx^ . 

We will study quantum vacuum effects brought about by the presence of parallel infinite 
plane boundaries (branes), located at y = a and y = b, a < b, with mixed boundary conditions 

Ay + Bydy'^ if (x) = 0, y = a,b, (20.3) 

and constant coefficients Ay, By. As a first stage we will consider the positive frequency Wight- 
man function. To apply the mode-sum formula we need the eigenfunctions for the problem 
under consideration. On the base of the plane symmetry of the problem these functions can be 
presented in the form 

- — ' (20.4) 



faix'') = =)fn{y), = (w,k). 



where co = y^k"^ + m^, k = |k|, and the separation constants are determined by the boundary 
conditions. Substituting eigenfunctions (|2U.4|) into the field equation, for the function fn{y) one 
obtains the following equation 

- e'^'^dy {e-''"dyfn) + (m2 + CR) /„ = m^e^-/,. (20.5) 
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For the AdS geometry one has a{y) = k£,y, z = e^^^^^ /k^, and R = —D{D + l)k'jy, where the 
AdS curvature radius is given hy l/kij. In this case the solution to equation (|20.5p for the region 
a < y < b satisfying the boundary condition at y = 6 is given by 

fn{y) = Cne^"/"^ gy{mnZa,rnnZ), (20.6) 
where the function g^(u,v) is defined by formula (|17.3|) . 



{D/2y - D{D + 1)C + m?/kl, Zj = e'^^^^ko, (20.7) 

and we use the barred notation ()4.2p with 

Aj = Aj + BjkDD/2, Bj = Bjkn, j = a, b. (20.8) 

We will assume values of the curvature coupling parameter for which u is real. For imaginary 
v the ground state becomes unstable [I39| . Note that for a conformally coupled massless scalar 
one has u = 1/2 and the cylinder functions in Eq. (j20.6p are expressed via the elementary 
functions. From the boundary condition on the brane y = b we receive that the eigenvalues nin 
have to be solutions to the equation 

Cf{Zb/Za,mnZa) = ji''\m„Za)Yj'\mnZb) - Yl)''\mnZa)J^J'\mnZ^) = 0. (20.9) 

The eigenvalues for m„ are related to the zeros of the function C^'^{r], z) as 

m„ = kD-f^,ne-'''^^^ = 7u,n/Za. (20.10) 

The coefficient Cn in Eq. (j20.6p is determined from the orthonormality condition 

rb 



[ (i2/e(2-^)-/„(y)/„,(y) = ^„„,, (20.11) 

J a 



and is equal to 



cl = :^!^Tf{7],u), U = J,,n, ll = Zb/Za. (20.12) 



2k£)Z 



a 



Note that, as we consider the quantization in the region between the branes, Za ^ z ^ zi,, the 
modes defined by (j20.6p are normalizable for all real values of u from Eq. (j20.7p . 

Substituting the eigenfunctions ()20.4p into mode-sum (j9.8p , for the expectation value of the 
field product one finds 

W{x,x') = ^Z+i\n-3,2 J dke^'^^-^^/i,(7.,„)T;^ (r/,7.,n) , (20.13) 

where x = (x^,x^, . . . ,x^~^) represents the coordinates in {D — l)-hyperplane parallel to the 
branes, Ax = x — x', and 



K{u) = ^====-g^{u,uz/za)gv{u,uz' /Za), (20.14) 

^Ju^/zi + k^ 

with At = t- t'. 

To sum over n we will use summation formula ()4.14p . Using the asymptotic formulae for 
the Bessel functions for large arguments when v is fixed (see, e.g., [65]), we can see that for the 
function hu{u) from Eq. ()20.14p the condition ()4.5p is satisfied if 2; + + |Af| < 2zf,- As for 



95 



\u\ < k one has hy{ue^'^) = —hy{u), the condition (|4.12|) is also satisfied for the function hy{u). 
Note that hy{u) ~ n^~'''=o for u ^ and the residue term on the right of formula ()4.14p vanishes. 
Applying to the sum over n in Eq. (|20.13|) formula (j4.14|) . one obtains 

t.D~l( /\D/2 p 

1 f°° K(u)du 2 /""^ nauiuZa,UZb) 

duu- 



xG^^\uZa,uz)G^^\uZa,uz')cosh{At^v? - A:2)|, (20.15) 

where we use notations (|4.15|) . p2.3U|) . Note that we have assumed values of the coefficients Aj 
and Bj for which all zeros for Eq. (j20.9p are real and have omitted the residue terms. In the 
following we will consider this case only. 

In the way similar to that used before, it can be seen that the first integral in the figure 
braces in Eq. (j20.15p is presented in the form 

hy{u)du f°° ^-iAWu'^+k'^ 



ji^^\u)+Yy^\u) "Jo v^^n^ 



2zl Ji''\uza) K,{uz)K,{uz') 



X 



cosh(At\/ii2 _ p), (20.16) 
Substituting this into formula (|20.15p . the Wightman function can be written in the form [56] 

W{x,x') = iyW(x,xO- ^^''yr' [dke^^^- Tdn^-^-^^^"'^^^^ 



xG(f ) {uza, uz)G(f ) (uza, uz') cosh{Aty^u^ - k^). (20.17) 

Here the term 

iyW(x,2;') = Ws{x,x') + {ip{xMx'))a (20.18) 

does not depend on the parameters of the boundary at y = b and is the Wightman function for 
a single brane located at y = a. In formula (j20.18p . the term 



h.D-1 f /\D/2 I- POO -iAiVw^+fe^ 

Ws{x,x') = ^\ / dke'^^'^ / duu^==^Uuz)Uuz'), (20.19) 

does not depend on the boundary conditions and is the Wightman function for the AdS space 
without boundaries, and the second term on the right. 



Kyiuz)Ky{uz') 



Kl^'\uZa) 



- cosh(Af\/n2 - A;2)^ (20.20) 

is induced in the region z > Za by a single brane at 2: = Za when the brane z = z^ is absent. 
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The Wightman function in the region Za ^ z ^ zi, can also be presented in the equivalent 
form 

^ ^ ^ ^ 2D~HD J - 

xG^^^uzb, uz)G^J'^ {uzb, uz') cosh(Aty^I^"^)- (20.21) 

In this formula 

W^'\x, x') = Ws{x, x') + {ip{xMx'))f„ (20.22) 

where 

{^{xMx))b = J dke ^^^4^)^ 

^h^^^^fikkilcoshi^t^f^^^), (20.23) 

is the Wightman function for a single brane at y = b. The term (|20.23p is the boundary part 
induced in the region z < Zb hy the brane at y = 6. Combining two forms, formulae ()20.17p 
and (|20.2ip . we see that the expressions for the Wightman function in the region Za ^ z ^ Zb 
is symmetric under the interchange a h and 1^ <^ Kj^. Note that the expression for the 
Wightman function is not symmetric with respect to the interchange of the brane indices. The 
reason for this is that, though the background AdS spacetime is homogeneous, the boundaries 
have nonzero extrinsic curvature tensors and two sides of the boundaries are not equivalent. In 
particular, for the geometry of a single brane the VEVs differ for the regions on the left and on 
the right of the brane. Here the situation is similar to that for the case of a spherical shell on 
background of the Minkowski spacetime. 



20.2 Casimir densities for a single brane 

In this subsection we will consider the VEV of the EMT for a scalar field in the case of a single 
brane located aX z = Za- As it has been shown in the previous subsection the Wightman function 
for this geometry is presented in the form (j20.18p . The brane induced part {ip{x)ip{x'))a is given 
by formula (120. 20p in the region z > Za and by formula (120. 23p with replacement Zb — > Za in the 
region z < Za- For points away from the brane this part is finite in the coincidence limit and in 
the corresponding formulae for the Wightman function we can directly put x = x' . Introducing 
a new integration variable v = Vu'^ — k'^, transforming to polar coordinates in the plane {v, k) 
and integrating over angular part, the following formula can be derived 



By using this formula and Eq. (j2U.20p . the boundary induced VEV for the field square in the 
region z > is presented in the form [56j 

^ / ^ D-l '-^''''-^ Kliuz), Z>Za. (20.25) 

1^ 2«"i7r^/2r(D/2) io Kl''\uza) 

The corresponding formula in the region z < Za'is obtained from Eq. (j20.23p by a similar way 
and differs from Eq. (j20.25p by replacements ly ^ Ky. 
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The VEV of the EMT can be evaluated by substituting expressions for the positive frequency 
Wightman function and VEV of the field square into Eq. (j9.10p . First of all we will consider 
the region z > Za- The vacuum EMT is diagonal and can be presented in the form 

(0|i;'=|0) = (0|I^'=|0)s + (7;')a, (20.26) 

where 

uD+ltrk / r)\ poo 

= (^1 - -J duu^-^ f^%Uu)], (20.27) 



is the VEV for the EMT in the AdS background without boundaries, and the term [56) 

is induced in the region z > Za by a single boundary sA z = Za- For a given function g{v) the 
functions F^'^[g{v)] in formula (j20.28D are defined as 

+ (^' + ^' + 2)(4y_i) )g'(^) ' i = 0,l,...,i^-l, (20.29) 

F^^\g{v)] = -'i^g'\v) + ^{AC-l)vg{v)g\v) + 

+\ [^^ + ^'^ + KD{D + 1) - dV2] 9^{v), (20.30) 

and the expressions for the functions J^*-* {9{v)\ are obtained from those for F^*) [giv)] by the 
replacement v ^ iv. Note that the boundary- induced part (120. 28p is finite for points away the 
brane and, hence, the renormalization procedure is needed for the boundary-free part only. The 
latter is well investigated in literature. 

In a similar way, for the VEVs induced by a single brane in the region z < z^, by making 
use of expression (|20.23p (with replacement z^, Za)-, one obtains 

Note that VEVs (j20.28p . (|20.3ip depend only on the ratio z/za which is related to the proper 
distance from the brane by the equation 

z/za = e''^^y-''\ (20.32) 

As we see, for the part of the EMT corresponding to the coordinates in the hyperplane parallel 
to the branes one has (T^jy)^") ~ Of course, we could expect this result from the problem 
symmetry. It can be seen that the VEVs obtained above obey the continuity equation V^T-^ = 0, 
which for the AdS metric takes the form 

z^+^d, (z-^Tg) + DTl^ = 0. (20.33) 

For a conformally coupled massless scalar v = 1/2, and by making use of the expressions for 
the modified Bessel functions, it can be seen that {Tj')^"'^ = in the region z > Za and 

/T^r'Ua) _ _r,(rrO^{a) _ {kDz/Za)^+^ j"^ dt 
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in the region z < Za- Note that the corresponding energy-momentum tensor for a single Robin 
plate in the Minkowski bulk vanishes [32] and the result for the region z > Za is obtained by a 
simple conformal transformation from that for the Minkowski case. In the region z < Za this 
procedure does not work as in the AdS problem one has < z < Za instead oi —oo < z < Za 
in the Minkowski problem and, hence, the part of AdS under consideration is not a conformal 
image of the corresponding manifold in the Minkowski spacetime. 

The brane-induced VEVs given by equations (|20.28p and (120. 3ip . in general, can not be 
further simplified and need numerical calculations. Relatively simple analytic formulae can be 
obtained in limiting cases. First of all, as a partial check, in the limit /c^ — > the corresponding 
formulae for a single plate on the Minkowski bulk are obtained (see Ref. [32]). This can be 
seen noting that in this limit z/ ~ m/A;£i is large and by introducing the new integration variable 
in accordance with u = z/y, we can replace the modified Bessel functions by their uniform 
asymptotic expansions for large values of the order. The Minkowski result is obtained in the 
leading order. 

In the limit z ^ Za for a fixed ku expressions ()20.28p and (|20.3ip diverge. In accordance 
with (120. 32p . this corresponds to small proper distances from the brane. Near the brane the 
main contribution into the integral over u in Eqs. (j20.28p . (|20.3ip comes from large values of 
u and we can replace the modified Bessel functions by their asymptotic expressions for large 
values of the argument when the order is fixed (see, for instance, [65]). To the leading order this 
yields 

where ka is defined after formula (jl9.22p . Note that the leading terms for the components with 
i = 0,1, . . . , D — 1 are symmetric with respect to the brane, and the ^ - component has opposite 
signs for different sides of the brane. Near the brane the vacuum energy densities have opposite 
signs for the cases of Dirichlet {Ba = 0) and non-Dirichlet {Ba 7^ 0) boundary conditions. Recall 
that for a conformally coupled massless scalar the vacuum EMT vanishes in the region z > Za 
and is given by expression (j20.34p in the region z < Za- The latter is finite everywhere including 
the points on the brane. For large proper distances from the brane compared with the AdS 
curvature radius, k£)\y — a\ ^ 1, the boundary induced EMT vanishes as exp[2i/A;£)(a — y)] in 
the region y > a and as exp[/i:£)(2z^ + D){y — a)] in the region y < a. The same behavior takes 
place for a fixed y — a and large values of the parameter k^. In the large mass limit, ko, 
the boundary parts are exponentially suppressed. 

20.3 Two-brane geometry 

In this subsection we will investigate the VEVs for the field square and the EMT in the region 
between two branes. Taking the coincidence limit in the formulae for the Wightman function 
and using formula (j20.24p , for the VEV of the field square one finds [56j 

(0|v.^|0) = iO\^'\0)s + {^'h- ^^jX4%/2) 

du u^^^Q.ju{uZa, uzhjG^j^"^ {uzj,uz), (20.36) 



00 







where j = a,b provide two equivalent representations. The last term on the right of this formula 
is finite for points on the brane at y = j and diverges on the other brane. In the similar way, 
substituting the corresponding Wightman function from Eq. ()20.17p into the mode-sum formula 
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dnU]), we obtain 



t.D+1 D';k 

{0\Tt\0) = {0\Tt\0)s + {Thj ^ 



2^-i7r^/2r(D/2) 
X / duu'^-'^nj^{uZa,uzb)F^'\Gi^\uzj,uz)], (20.37) 





with the functions F(*)[c/(t;)] from Eqs. (120.291) . (j20.30p . The formulae for the cases j = a,b are 
obtained from each other by the replacements a b, ^ Ky. The VEVs in the region between 
the branes are not symmetric under the interchange of indices of the branes. As it has been 
mentioned above, the reason for this is that, though the background spacetime is homogeneous, 
due to the non-zero extrinsic curvature tensors for the branes, the regions on the left and on the 
right of the brane are not equivalent. By the same way as for the case of a single brane, it can 
be seen that in the limit /cd — > from formulae (|20.36p and (|20.37p the corresponding results 
for two plates geometry in the Minkowski bulk. For a conformally coupled massless scalar field 
V = 1/2 and the formulae for the VEVs in the region between the branes can also be obtained 
from the corresponding formulae for parallel Robin plates in the Minkowski bulk by conformal 
transformation |140j . 

Now we turn to the interaction forces between the branes. The vacuum force acting per unit 
surface of the brane ai z = Zj is determined by the ^ - component of the vacuum EMT at this 
point. The corresponding effective pressures can be presented as the sum p^^^ = p'p +?'(int)' 
j = a, h, where the first term on the right is the pressure for a single brane at z = Zj when the 
second brane is absent. This term is divergent due to the surface divergences in the VEVs. The 
second term is determined by the last term on the rhs in ()20.37p evaluated at z = Zj . This term is 
the pressure induced by the presence of the second brane. Using the relations Gu\u, u) = —Bj, 
Gu^'{u,u) = Aj/u, one finds 



0) _ 



(int) 2^7r^/2r (f ; 







X [(x^ + 2m^/kl,) B] - D{4C - l)AjBj - Aj] . (20.38) 

In dependence of the values for the coefficients in the boundary conditions, these effective pres- 
sures can be either positive or negative, leading to repulsive or attractive forces. Note that due 
to the asymmetry in the VEV of the EMT, the interaction forces acting on the branes are not 
symmetric under the interchange of the brane indices. It can be seen that the vacuum effective 
pressures are negative for Dirichlet scalar and for a scalar with Aa = Ai, = and, hence, the 
corresponding interaction forces are attractive for all values of the interbrane distance. 

Let us consider the limiting cases for the interaction forces described by Eq. (j20.38p . For 
small distances compared with the AdS curvature radius, koib — a) <C 1, the leading terms 
are the same as for the plates in Minkowski bulk. In particular, in this limit the interaction 
forces are repulsive for Dirichlet boundary condition on one brane and non-Dirichlet boundary 
condition on the other, and are attractive for all other cases. For large distances between the 
branes, k^ih — a) ^ 1 (this limit is realized in the Randall-Sundrum model), by using the 
expressions for the modified Bessel functions for small values of the argument, one can see that 
Ptlt) ~ (^a/22;fe)^+^'' andpjl'^^) ~ {zal2zi,f'' . 

Now we consider the application of the results described in this section to the Randall- 
Sundrum braneworld models |138j based on the AdS geometry with one extra dimension. The 
fifth dimension y is compactified on an orbifold, /Z2 of length L, with —L ^ y ^ L. The 
orbifold fixed points at y = a = and y = b are the locations of two 3-branes. We will allow 



100 



these submanifolds to have an arbitrary dimension D. The metric in the Randall-Sundrum 
model has the form (|20.ip with cj(y) = fcclyl- The corresponding boundary conditions for an 
untwisted scalar field are in form (|2U.3p with 



Aj/Bj = -{n^^^cj + 4DCkD)/2, (20.39) 

and respectively 

Aj/Bj = \d{1 - 40 - n^^^Cj/kD] /2, (20.40) 

where Cj are the brane mass terms. For a twisted scalar Bj = 0, which corresponds to Dirichlet 
boundary conditions on both branes. Recently the EMT in the Randall-Sundrum braneworld 
for a bulk scalar with zero brane mass terms ci and C2 is considered in |141j . where a general 
formula is given for the unrenormalized VEV in terms of the differential operator acting on the 
Green function. In our approach the application of the GAPF allowed to extract manifestly 
the part due to the AdS bulk without boundaries and for the points away from the boundaries 
the renormalization procedure is the same as for the boundary-free parts. In addition, the 
boundary-induced parts are presented in terms of exponentially convergent integrals convenient 
for numerical calculations. 

From the point of view of embedding the Randall-Sundrum type braneworld models into 
a more fundamental theory, such as string/M theory, one may expect that a more complete 
version of this scenario must admit the presence of additional extra dimensions compactified on 
an internal manifold. The results discussed in this section can be generalized for the geometry 
of two parallel branes of codimension one on background of (D-l- 1) -dimensional spacetime with 
topology AdSoi+i x S and the line element 

^ e-^kny^^^^^f^^^a _ e-'^kny^.^^x'dX^ - dy^, (20.41) 

where ry^o- = diag(l, — 1, . . . , — 1) is the metric tensor for Di -dimensional Minkowski spacetime 
and the coordinates X*, i = 1, . . . , D2, cover the manifold S. The quantum effective 
potential and the problem of moduli stabilization in the orbifolded version of this model with 
zero mass parameters on the branes were discussed recently in [142j . In particular, it has been 
shown that one-loop effects induced by bulk scalar fields generate a suitable effective potential 
which can stabilize the hierarchy between the gravitational and electroweak scales. In [57] the 
Wightman function, the VEVs of the field square and the EMT for a scalar field brought about 
by the presence of two parallel infinite branes, located at y = a and y = b with boundary 
conditions (j20.3p . are investigated in a way similar to that used in this section. 



21 Casimir densities for spherical branes in Rindler-like space- 
times 

In previous section we have considered braneworld models on background of AdS spacetime. It 
seems interesting to generalize the study of quantum effects to other types of bulk spacetimes. In 
particular, bulk geometries generated by higher-dimensional black holes are of special interest. 
In these models, the tension and the position of the brane are tuned in terms of black hole 
mass and cosmological constant and brane gravity trapping occurs in just the same way as in 
the Randall-Sundrum model. Though, in the generic black hole background the investigation 
of brane-induced quantum effects is technically complicated, the exact analytical results can 
be obtained in the near horizon and large mass limit when the brane is close to the black 
hole horizon. In this limit the black hole geometry may be approximated by the Rindler-like 
manifold (for some investigations of quantum effects on background of Rindler-like spacetimes 
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see |143t 1144] and references therein). In this section, by using the GAPF, we investigate the 
Wightman function, the VEVs of the field square and the EMT for a scalar field with an arbitrary 
curvature coupling parameter for two spherical branes in the bulk Ri x S^~^, where Ri is a 
two-dimensional Rindler spacetime [58j . 



21.1 Wightman function 

Consider a scalar field (p{x) propagating on background of {D + l)-dimensional Rindler-like 
spacetime Ri x S^^^. The corresponding metric is described by the line element 

ds^ = ^'^dT^ -df -rj^dJ:l_^, (21.1) 

with the Rindler-like (r, ^) part and dT,'jj_^ is the line element for the space with positive constant 
curvature with the Ricci scalar R = n{n + l)/r|^, n = D — 2. Line element (121. ID describes the 
near horizon geometry of {D + l)-dimensional topological black hole with the line element 

ds^ = AH{r)dt^ - - r^d^l_,, (21.2) 

where Ajj{r) = k + r'^/P — /Pr^ and the parameter k classifies the horizon topology, with 
/c = 0, —1, 1 corresponding to flat, hyperbolic, and elliptic horizons, respectively. The parameter 
I is related to the bulk cosmological constant and the parameter tq depends on the mass of the 
black hole. In the non-extremal case the function Auir) has a simple zero aX, r = th, and in 
the near horizon limit, introducing new coordinates r and p in accordance with 

T = A'H{rH)t/2, T-rH = J^u{rH)ei^. (21.3) 

the line element is written in the form ()21.ip . Note that for a (D + l)-dimensional Schwarzschild 
black hole one has Auir^ = 1 — {rul'^Y^ and, hence, A'j^{r}{) = n/rn- 

We will assume that the field satisfies Robin boundary conditions (jl9.ip on the hypersurfaces 
^ = a and ^ = 6, a < 6. In accordance with (j21.3p . the hypersurface ^ = j corresponds to 
the spherical shell near the black hole horizon with the radius rj = + A'^{rH)j'^ In 
the corresponding braneworld scenario based on the orbifolded version of the model the region 
between the branes is employed only and the ratio Aj/Bj for untwisted bulk scalars is related 
to the brane mass parameters Cj of the field by the formula |144j 

Aj/B, = ^{cj-C/j), j = a,b. (21.4) 

For a twisted scalar Dirichlet boundary conditions are obtained on both branes. 

In hyperspherical coordinates (see Section [12] for the notations) the eigenfunctions in the 
region between the branes can be written in the form 

Mx) = P„GfJ{Xib,XiC)Yimk;^A)e-"^\ (21.5) 

where the function Gu\x, y) is defined by formula ()12.30p with the barred notations (j4.2p . where 
Aj = Aj, Bj = Bj/j, and 



Xl = —Jl(l+n)+Cn(n+l)+m'^r%. (21.6) 
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From the boundary conditions we find that the eigenvalues for uj are roots to the equation 
Ziuii^ia, Xib) = 0, where the function Zii^{u,v) is defined by formula (|5.14p . and, hence, uj = 
ilg = ^}s{Xia, Xib), s = 1, 2, . . .. For the normahzation coefficient one has 



2 



Pi 



r 



N{mk)lfJ{\ib)d^ZU\ia,\ih) 



(21.7) 



Substituting eigenfunctions (j21.5p into the mode-sum formula, for the Wightman function 
in the region between the branes one finds 

W{x,x') = ^ V(2/ + n)q"^^(cos e) V 

U i^i lTJ{kh)d^ZUMa,\ib) 

xGf^{\A\ii)Gf^{\A\ii')\ , (21.8) 

with Ar = T — t' . The application of summation formula ()5.18p to the sum over s leads to the 
result 



l-D oo 

W{x,x') = Wo(x,x') + (v5(2;)v9(x'))(^')--^y(2/ + n)C"/^(cos( 



imSn 

^ 1=0 



X / (kunjU>^ia,Xib)Gl^HXij,XiOG^j\Xij,\ie)coshiu;AT), (21.9) 



with j = b, where W(){x, x') is the Wightman function for the geometry without boundaries and 
the part 

(^(x)v.(x'))(^) = -^f:i2l + n)Cr^\cos9) Tdu^^^ 

xI^{XiC)lUXiC') cosh(u;Ar) (21.10) 

is induced in the region ,^ < 5 by the presence of the brane at ^ = 6. Note that the representation 
(|21.9p is vahd under the condition a^e'^^' < < t^el^"^'. As it has been shown in |144j . the 
Wightman function for the boundary-free geometry may be written in the form 

Woix,x') = VFo(x,x')-^^ V(2^ + n)Cf/'(cos0) 



POO 

/ du;e-'^^cos{u;AT)KiUWK^U^le), (21.11) 
Jo 



where Wo{x, x') is the Wightman function for the bulk geometry B? x S"^"^. Outside the horizon 
the divergences in the coincidence limit of the expression on the right of (j21.1ip are contained in 
the first term. It can be seen that the Wightman function in the region between the branes can 
also be presented in the form (j21.9p with j = a. In this representation {ip{x)Lp{x'))^°^ is induced 
in the region ^ > a by the presence of the brane at ^ = a and the corresponding expression is 
obtained from (j21.10p by the replacements 6 — > a, K^. 
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21.2 Vacuum densities 



In the coincidence limit, from the formulae for the Wightman function one obtains two equivalent 
forms for the VEV of the field square: 



VA / dconjU>^ia,Xib)Gl^^\Xij,XiO, (21.12) 

^0 



irSn 

^ 1=0 

corresponding to j = a and j = b, and |Oo) is the amplitude for the vacuum without boundaries. 
The coefficient Di in this formula is defined by relation (|12.14p . The VEV (Oo|9J^|Oo) is obtained 
from the corresponding Wightman function given by (j21.11|) . For the points outside the horizon, 
the renormalization procedure is needed for the first term on the right only, which corresponds 
to the VEV in the geometry x 5''^"^. This procedure is realized in ^144j on the base of the 
zeta function technique. 

In ()21.12p . the part ((^^ is induced by a single brane at ^ = j when the second brane is 
absent. For the geometry of a single brane at ^ = 6, from ()21.10p one has 

l-D oo oo K'^'')(\,h\ 

The expression for i^ip^)^"^ is obtained from (j21.13p by the replacements 6 — > o, I^; «^ K^^. The 
last term on the right of formula (jl9.13|) is induced by the presence of the second brane. It 
is finite on the brane at ^ = j and diverges for points on the other brane. By taking into 
account the relation g[^\u,u) = Bj/j, we see that for the Dirichlet boundary condition this 
term vanishes on the brane i = j. 

By using the formulae for the Wightman function and the VEV of the field square, one 
obtains two equivalent forms for the VEV of the EMT, corresponding to j = a and j = b (no 
summation over i): 



1=0 



X / cLunj^{Xia,Xib)F^'\G^J\Xij,XiO]. (21.14) 







The functions F^^^[g{z)] in this formula are defined by relations (|18.15p and (|18.16p with A — > A; 
and g{z) = G\i\Xii,z). In formula (l2LTiD . 

l-D oo „oo 

{Q^\T^\%) = 5^^y^DiXU du;sinh7ru;/«[i^,^(A;0], (21-15) 



1=0 



with f^'^\g{zy\ defined by Eq. (jlS.lSp . is the corresponding VEV for the vacuum without bound- 
aries. The term {T^)^^'> is induced by the presence of a single spherical brane located at ^ = j. 
For the brane at ^ = 6 and in the region ^ <b one has (no summation over i) 

„1-D oo „oo k'^^)(\,U\ 

^Sd ^ Jo I^iXib) 

For the geometry of a single brane at = a, the corresponding expression in the region ^ > a 
is obtained from (j21.16p by the replacements 6 — > a, <^ K^^. It can be easily seen that for a 
conformally coupled massless scalar the boundary induced part in the EMT is traceless. 
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Vacuum forces acting on the branes are presented in the form of the sum of self-action and 
interaction terms. The vacuum effective pressures corresponding to the interaction forces are 
obtained from the third summand on the right of Eq. (j21.14|) , taking i = k = l,^=j 



^' int) = ^-VE^' / dco[{Xff + u;')(3] + ACP,-l] %^(A;a,A^6), (21.17) 

with (3j = Bj/jAj. The interaction force acts on the surface ^ = a + O for the brane at ^ = a and 
on the surface ^ = 6 — for the brane at ^ = 6. In dependence of the values for the coefficients in 
the boundary conditions, effective pressures p9.23|) can be either positive or negative, leading 
to repulsive or attractive forces, respectively. For Dirichlet or Neumann boundary conditions on 
both branes the interaction forces are always attractive. For Dirichlet boundary condition on one 
brane and Neumann boundary condition on the other one has > ^ interaction forces 

are repulsive for all distances between the branes. The investigation of the vacuum densities in 
various asymptotic regions of the parameters can be forund in [58^ I144j . 

22 Radiation from a charge moving along a helical orbit inside 
a dielectric cylinder 

The radiation from a charged particle moving along a helical orbit in vacuum has been widely 
discussed in literature. This type of electron motion is used in helical undulators for generating 
electromagnetic radiation in a narrow spectral interval at frequencies ranging from radio or 
millimeter waves to X-rays (see, for instance, p3S])- The unique characteristics, such as high 
intensity and high collimation, have resulted in extensive applications of this radiation in a wide 
variety of experiments and in many disciplines. In this section we apply the GAPF for the 
investigation of the radiation on the lowest azimuthal mode by a charged particle moving along 
a helical orbit inside a dielectric cylinder [59]. 

Consider a dielectric cylinder of radius pi and dielectric permittivity Eq and a point charge q 
moving along the helical trajectory of radius po < pi. We assume that the system is immersed in 
a homogeneous medium with permittivity ei. The velocities of the charge along the axis of the 
cylinder and in the perpendicular plane we will denote by V|| and v±, respectively. In a properly 
chosen cylindrical coordinate system (p, (p, z) with the z-axis along the cylinder axis, the vector 
potential of the electromagnetic field is presented in the form of the Fourier expansion 

Ai{r,t)= e^'-('^--o*) / dk,e'''^^'-^ll'U^i{k,,p), (22.1) 

m=-oo ''-'^ 

where ujq = v±/pQ is the angular velocity of the charge. The radiation field inside the dielectric 
cylinder is presented as a sum over the eigenmodes of the cylinder. Unlike to the case of the 
waveguide with perfectly conducting walls, here the eigenmodes with m ^ are not decomposed 
into independent transverse electric (TE or M-type) and transverse magnetic (TM or E-type) 
parts. This decomposition takes place only for the m = mode. The radiation on this mode is 
present under the condition 

/3ij| < 1 < /3o||, All = ^llV^/c, i = l,2. (22.2) 
With these conditions the equation for the corresponding TM eigenmodes for k^ has the form 

so|Ai||^ + e,Ao|^^ = 0, X. = k.J^l. (22.3) 
Jo(Ao/Oi) Ao(|Ai|pi) V *ii 
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The corresponding eigenvalues for kz we will denote by ^kl , s = 1, 2, . . .. The eigenmodes of 
the M-type are solutions of the equation 

J^(Ao,.i) K',{\\M _ 
Jq{XqPi) Ao(|Ai|pi) 

and we will denote them by kz = it/ci^^ s = 1, 2, . . .. 

For the radiation intensities on the mode m = inside the dielectric waveguide we have 

Im=o = Imlo + -^^io' (22.5) 
where the contributions of the TM and TE modes are given by the formulae [59] 



(M) _ 2q^vl ^ JfjpoX?^^) , . 

Here we have introduced the notation 



Af) = fef)^/32, -1, F = E,M. (22i 



In formulae (I22.6P and (I22.7P the upper limit for the summation over s, which we will denote 
by Sm, is determined by the dispersion law of the dielectric permittivity = Eq^uj) through the 
condition eo{iL>'^^) > c'^/wy • The summands with a given s describe the radiation field with the 

frequency Ws = v^^kg . As the modes kg are not explicitly known as functions on s, formulae 
()22.6p . (j22.7p are not convenient for the evaluation of the corresponding radiation intensities. 
More convenient form may be obtained by making use of the GAPF. 

Let us derive a summation formula for the series over zeros of the function 

Caiv,z) = VMiz),Koir]z)}, (22.9) 

where and in what follows for given functions F{z) and G{z) we use the notation 

V{F{z), G{t]z)} = F{z)G'{r]z) + ar]F' {z)G{t]z), (22.10) 

with a ^ 1 and r] being real constants. In the GAPF we substitute 

«(^) = i/(^) '°'''°'t'^'°''"" ' 

For the combinations of the functions entering in the GAPF one has 

Let us denote positive zeros of the function G^i'q^z) by kg, s = 1,2,..., assuming that these 
zeros are arranged in the ascending order. Note that, for z ^ 1 we have CaiVj^) ~ ^oiv^) ~ 
ar]zKQ{riz)/2 < and, hence, kg > 1. By using the asymptotic formulae for the cylinder 
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functions for large values of the argument, it can be seen that for large values s one has kg ~ 
— arctan(l/a;r7) +7r/4 + 7rs. For the derivative of the function Ca{i], z) at the zeros kg one obtains 



Jo{z 



a-l)J^(z) 



k.e . 



(22.13) 



In particular, it follows from here that the zeros are simple. Assuming that the function f{z) is 
analytic in the right half-plane, for the residue term in the GAPF one finds 



Res g{z) 

z=ks 

where we have introduced the notation 

Po.{z, 



vr 



'Pa{ks)f(ks), 



2a/z 



a{l + ar]-^)J'^{z) + (a - l)4[z) ' 



(22.14) 



(22.15) 



Substituting the expressions for the separate terms into the GAPF, we obtain the following 
result [591 



lim 



y2PUks)fiks)- / dxfix) = — / 



dx 



Vc,{Ko{x),Hl^\r]x)} 



+f{-ix) 



14{/o(x),//i')(r?x)} 
Va{Ko{x),Hl,^\vx)} ] ^ 
Va{Io{x),Hl,^\vx)}\ ' 



22.16) 



where sq is defined by the relation kg^ < xq < A^sq+i. This formula is valid for functions f{z) 
obeying the condition 

\f{z)\ < e(x)e"l^l, z = x + iy, \z\ oo, (22.17) 

where c < 2 and e(x) — > for x oo. Formula ()22.16p is further simplified for functions 
satisfying the additional condition f{—ix) = —f{ix): 



where 



lim 



9air],x) 



so 



J2Pa{ks)f{ks 



XO 



dxfix) 



4ia 



vr^ Jo x'^ga{r],x)' 



(22.18) 



Ilix) [Jfirjx] + Y.^irjx)] + aWlfix) [JUtix) + Yi{r,x)] 
-2aT]Io{x)Ii{x) [Jo{T]x)Ji{r]x) + Yo{r]x)Yi{r]x)] . 



(22.19) 



Note that we have denoted ga{f],x) = \V{lQ{z),Hl^\riz)^\'^ and this function is always non- 
negative. 

Now, summation formulae for the series over m = TM and TE modes are obtained taking 



\ "Sii 



1 



(22.20) 



In formula (|22.16p we choose a = Eq/ei, f{z) = zJq{zpq/ pi) for the waves of the E-type, and 
a = 1, f{z) = zjf{zpo/ pi) for the waves of the M-type. For both types of modes one has 
f{—ix) = —f{ix) and we can use the version of the summation formula given by (j22.18p . In the 
intermediate step of the calculations, it is technically simpler instead of considering the dispersion 
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of the dielectric permittivity to assume that in formulae ()22.6p . (j22.7p a cutoff function ip^{Xs ) 
is introduced with // being the cutoff parameter and ipQ = 1 (for example, ip^{x) = exp(— /ix)), 
which will be removed after the summation. In this way, after the application of the summation 
formula we find the following results 



7-(E) 

'm=0 



(M) 
m=0 



dx—ip^{x)jQ{pQx) + -^T-^ 

£0 TT^Pl Jo 



dx 



dx 



dx 



ll{xpo/pi 



Ilixpo/pi) 
(/3o|| - ^)xgi{v,x) 



(22.21) 
(22.22) 



where the function (7q.(?7, x) is defined by formula (j22.19p . In the first terms in the square brack- 
ets replacing the integration variable by the frequency oo = v\\x / \ and introducing the 

physical cutoff through the condition /3o|| > 1 instead of the cutoff function, we see that these 
terms coincide with the radiation intensities on the harmonic m = for the waves of the E- and 
M-type in the homogeneous medium with dielectric permittivity £q. The second terms in the 
square brackets are induced by the inhomogeneity of the medium in the problem under consid- 
eration. Note that, unlike to the terms corresponding to a homogeneous medium, for po < Pi 
the terms induced by the inhomogeneity are finite also in the case when the dispersion is absent: 
for large values of x the integrands decay as exp[— 2x(l — po/ Pi)] (for this reason we have re- 
moved the cutoff function from these terms). In particular, from the last observation it follows 
that under the condition (1 — po/pi) ^ fy/wd, where ojd is the characteristic frequency for the 
dispersion of the dielectric permittivity, the influence of the dispersion on the inhomogeneity 
induced terms can be neglected. Note that in a homogeneous medium the corresponding radia- 
tion propagates under the Cherenkov angle 6c = arccos(l//3o|j) and has a continuous spectrum, 
whereas the radiation described by ()22.2ip . ()22.22p propagates inside the dielectric cylinder and 
has a discrete spectrum with frequencies uJs ■ As the function ga{ri^x) is always non-negative, 
from formulae (j22.2ip .( 122.22p we conclude that the presence of the cylinder amplifies the m = 
part of the radiation for the waves of the E-type and suppresses the the radiation for the waves 
of the M-type. In the helical undulators one has v±_ <^ and the contribution of the TM waves 
dominates. 



23 Summary 

The Abel-Plana summation formula is a powerful tool for the evaluation of the difference between 
a sum and the corresponding integral. This formula has found numerous physical applications 
including the Casimir effect for various bulk and boundary geometries. However, the applications 
of the APF in its standard form are restricted to the problems where the normal modes are 
explicitly known. In the present paper we have considered a generalization of the APF, proposed 
in [27j , which essentially enlarges the application range and allows to include problems where the 
eigenmodes are given implicitly as zeros of a given function. Well known examples of this kind are 
the boundary-value problems with spherical and cylindrical boundaries. The generalized version 
contains two meromorphic functions f{z) and g{z) and is formulated in the form of Theorem 
1. The special choice g{z) = — if (z) cot irz gives the APF with additional residue terms coming 
from the poles of the function f{z). We have shown that various generalizations of the APF 
previously discussed in literature are obtained from the GAPE by specifying the function g{z). 
Further we consider applications of the GAPE to cylinder functions. These applications include 
summation formulae for series over the zeros of these functions and their combinations (Sections 
I3|4|5p and formulae for integrals involving cylinder functions (Sections I6|7|8p . In the second part 
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of the paper, including Sections [U H221 we outline the applications of the summation formulae 
over the zeros of cylinder functions for the evaluation of the VEVs for local physical observables 
in the Casimir effect and the radiation intensity from a charge moving along a helical trajectory 
inside a dielectric cylinder. 

In Section [3l in the GAPF choosing the function g{z) in the form (j3.2p . we derive two 
types of summation formulae, (|3.18p . (|3.33p . for the series over the zeros Xi,^k of the function 
AJy{z) + BzJl{z). This type of series arises in a number of problems of mathematical physics 
with spherically and cylindrically symmetric boundaries. As a special case they include Fourier- 
Bessel and Dini series (see |64j). Using formula (13.18p . the difference between the sum over 
zeros Xy^k and the corresponding integral is presented in terms of an integral involving modified 
Bessel functions plus residue terms. For a large class of functions f{z) the last integral converges 
exponentially fast and, in particular, is useful for numerical calculations. The APF is a special 
case of formula (j3.18p with v = 1/2, A = 1^ B = Q and for an analytic function f{z). Choosing 
in ()3.18p z^=l/2,yl = l,i? = 2we obtain the APF in the form (I2.19P useful for fermionic field 
calculations. Note that formula (j3.18p may also be used for some functions having poles and 
branch points on the imaginary axis. The second type of summation formula, (13.33p . is valid for 
functions satisfying condition ()3.27p and presents the difference between the sum over zeros Xy^^ 
and the corresponding integral in terms of residues over poles of the function f{z) in the right 
half-plane (including purely imaginary ones). This formula may be used to summarize a large 
class of series of this type in finite terms. In particular, the examples we found in literature, 
when the corresponding sum may be presented in a closed form, are special cases of this formula. 
A number of new series summable by this formula and some classes of functions to which it can 
be applied is presented. 

As a next application of the GAPF, in Section [5] we consider series over zeros z = ^^^^ of the 
function Jy°'\z)Yu'^\riz) — Yp"'\z)jl)'\r]z), where the barred notations are defined by relation 
(14. 2p . The corresponding results are formulated in the form of Corollary 2 and Corollary 3. 
Using formula (j4.14p . the difference between the sum over 7,^^ and the corresponding integral 
can be expressed as an integral containing modified Bessel functions plus residue terms. For 
a large class of functions h{z) this integral converges exponentially fast. The formula of the 
second type, (I4.18p . allows to summarize a class of series over the zeros ju^k in closed form. 
To evaluate the corresponding integral, the formula can be used derived in section [8j This 
yields to another summation formula, (|4.2ip . containing residue terms only. The examples we 
have found in literature when the corresponding sum was evaluated in closed form are special 
cases of the formulae considered here. We present new examples and some classes of functions 
satisfying the corresponding conditions. In Section [Sj taking in the GAPF the functions in the 
form (j5.3p . we derive summation formula (j5.12p for series over the zeros z = uj^ oi the function 
AKiz{r]) + Br]drfKiz{r]). These zeros are eigenfrequencies for scalar and electromagnetic fields in 
the geometry of a uniformly accelerated plane boundary and the corresponding VEVs contain 
summation over them. For the geometry of two plane boundaries the eigenfrequencies in the 
region between the planes are zeros z = oi the function k\'^ {u)!^^ {v) — l\'^\u)Kf'^ {v). 
Summation formula ()5.18p for series over these zeros is obtained from ()2.1ip with the functions 
f{z) and g[z) given by (j5.15p . 

In Section [6] we consider relations for the integrals of type 

/•oo /'OO 

/ F{x)Ju{x)dx, I F{x)[Ju{x) COS 5 + Yi,{x)sm.S\dx, 
Jo Jo 

which are obtained from the GAPF. The corresponding formulae have the form ()6.4p . (16. 7p and 
(|6.18p . In particular, formula (j6.4p is useful to express the integrals containing Bessel functions 
with oscillating integrand through the integrals of modified Bessel functions with exponentially 
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fast convergence. Formulae (j6.7|) and (j6.18p allow to evaluate a large class of integrals involving 
the Bessel function in closed form. In particular, the results obtained in [68] are special cases 
of these formulae. Illustrating examples of applications of the formulae for integrals are given 
in Section [7] (see (|7.3p - (|7.7p and (j7.13p - (j7.17p ). Looking at the standard books (see, e.g., [13] . 
|64j-[72j) one will find many particular cases which follow from these formulae. A number of new 
integrals can be evaluated as well. We also consider two examples of functions having purely 
imaginary poles, (j7.19p and (|7.23p . with corresponding formulae (j7.20p and (j7.24p (for two special 
cases of these formulae see [64]). By choosing the functions f{z) and g{z) in accordance with 
(18. ip . formulae (18. 7p and (I8.16P for the integrals of the type 



can be derived from the GAPF. The corresponding results are formulated in the form of Theorem 
5 and Theorem 6 in Section |8l Several examples for integrals of this type we have been able 
to find in literature are particular cases of formula (j8.7p . New examples when the integral is 
evaluated in finite terms are presented. Some classes of functions are distinguished to which the 
corresponding formulae may be applied. 

In the second part of the paper, started from Section [U physical applications of the summa- 
tion formulae for series over zeros of cylinder functions obtained from the GAPF are reviewed. 
The general procedure for the evaluation of the VEVs of the field square and the EMT on 
manifolds with boundaries based on the point-splitting technique is described. This VEVs are 
obtained from the corresponding positive frequency Wightman function in the coincidence limit 
by using formulae ()9.9p . (jD.lOp . In these formulae instead of the Wightman function any other 
two-point function can be chosen. Our choice of the Wightman function is related to that it also 
determines the response of the Unruh-DeWitt type particle detectors. Mode-sum ()9.8p contains 
a summation over the eigenmodes of the problem under consideration. For curved boundaries 
these eigenmodes are given implicitly as zeros of the eigenfunctions for the corresponding bound- 
ary value problem or their combinations. The second difficulty for the direct evaluation of the 
mode-sum is related to that it diverges in the coincidence limit and exhibits very slow conver- 
gence characteristics when the arguments of the Wightman function are close to each other. In 
addition, the separate terms in the mode-sum are strongly oscillatory for higher modes. The 
application of summation formulae obtained from the GAPF enables to present the mode-sum in 
terms of integrals and, hence, in the corresponding procedure the explicit form of the eigenmodes 
is not necessary. These formulae explicitly extract from the VEVs of local physical observables 
the parts corresponding to the bulk without boundaries and the boundary-induced parts are 
presented in terms of integrals which are exponentially convergent for the points away from 
the boundaries. As a result, in the coincidence limit the renormalization is necessary for the 
boundary-free parts only and this procedure is the same as that in quantum field theory without 
boundaries. Note that, by using summation formulae obtained from the GAPF, the mode-sums 
for more general class of objects such as the heat and cylinder kernels can be evaluated as well. 

First, in Section [10] we consider examples where the eigenmodes are explicitly known and 
the application of the APF in standard form enables to obtain the renormalized VEVs. These 
examples include the problems for the evaluation of the VEVs in topologically non-trivial space 

X and for the geometry of two parallel plates with Dirichlet and Neumann boundary condi- 
tions on them. On these simple examples we have seen two important advantages of the applica- 
tion of the APF. First, this formula enables to extract explicitly from the VEVs the Minkowskian 
part and second, to present the parts induced by the non-trivial topology/boundaries in terms 
of rapidly convergent integrals. Already in the case of two parallel plates with Robin boundary 
conditions the eigenmodes are not explicitly known and the procedure for the evaluation of the 
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VEVs based on the APF needs a generalization. The corresponding problem is discussed in 
Section [TTJ In this problem the eigenmodes for the projection of the wave vector perpendic- 
ular to the plates are solutions of transcendental equation For the summation of the 
mode-sums over these eigenmodes we have derived formula (|11.7p by making use of the GAPF. 
The application of this summation formula enables to extract from the Wightman function the 
part corresponding to the geometry of a single plate (see formula (jll.9p ). Similar decomposed 
formulae are obtained for the VEVs of the field square and the EMT. 

In Section [12] we consider the scalar Casimir densities inside a spherical boundary and in the 
region between two concentric spherical boundaries on background of global monopole spacetime 
described by line element (jl2.ip . It is assumed that on the bounding surfaces the scalar field 
obeys Robin boundary conditions. For the region inside a single sphere the eigenmodes are zeros 
of the function Jy^{x)^ where the order of the Bessel function is defined by formula ()12.4p and 
the coefficients in the barred notation are related to the coefficients in the boundary condition 
by formula (112. 6p . For the evaluation of the corresponding mode-sum we have used summation 
formula (j3.18p . The term with the integral on the left of this formula corresponds to the 
Wightman function for the global monopole bulk without boundaries and the term with the 
integral on the right is induced by the spherical boundary. In the region between two spheres the 
eigenmodes are zeroes of the function C^^{r], z) defined by formula ()4.ip and for the summation 
over these eigenmodes formula (j4.14p is used. We have seen that the term with the integral 
on the left of this formula corresponds to the Wightman function in the region outside a single 
sphere and the term with the integral on the right is induced by the presence of the outer sphere. 

In section (I13p we have investigated the vacuum densities for a fermionic field induced by 
spherical boundaries in a global monopole background. It is assumed that the field obeys bag 
boundary condition. Both regions inside a single sphere and between two spheres are considered. 
For the first case the eigenmodes are zeros of the function Ju^{x) with the tilted notation defined 
by (jl3.7p . By special choice of the function g{z) in the GAPF we have derived formula (jl3.17p 
for the summation over these zeros. The application of this formula to the mode-sum for the 
VEVs enables to extract the parts corresponding to the boundary-free bulk without specifying 
the form of the cutoff function. In the region between two spheres the eigenmodes are zeros 
of the function Cl^{r),x) defined by formula (|13.29p . By special choice of the functions f{z) 
and g{z) in the GAPF given by ()13.35p . in subsection 113.21 we have derived summation formula 
()13.42p for the series over these zeros. This formula presents the VEV of the EMT in the form 
of the sum of a single sphere and second sphere induced parts. Scalar and fermionic Casimir 
densities for spherical boundaries in the Minkowski bulk are obtained from the corresponding 
results for the global monopole geometry as special cases with a = 1. In Section HM. by making 
use of the GAPF we have investigated the local properties of the electromagnetic vacuum inside a 
perfectly conducting spherical shell and in the region between two spherical shells on background 
of Minkowski spacetime. 

Vacuum polarization by cylindrical boundaries is considered in Sections 1151 1161 [T71 We start 
with the case of a scalar field on background of the cosmic string spacetime with line element 
(jlS.ip and with coaxial cylindrical boundary. On the boundary the field obeys Robin boundary 
condition. In the region inside the cylindrical shell the eigenmodes are zeros of the function 
Jq|„|(z), where the parameter q is determined by the planar angle deficit and the coefficients in 
the barred notation are defined by relation (jlS.Sp . The application of the GAPF to the mode- 
sum of the Wightman function extracts the part corresponding to the cosmic string geometry 
without boundaries and the boundary induced part is presented in the form which can be 
directly used for the evaluation of the VEVs for the field square and the EMT. In Section [16] 
similar consideration is done for the electromagnetic field with perfectly conducting boundary 
conditions. The corresponding eigenmodes are the zeros of the function Jq|„|(z) for the TM 
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waves and the zeros of the function J', Jz) for the TE waves. The formulae for the VEVs 

q\n\\ I 

induced by a cyhndrical boundary in the Minkowski bulk are obtained from the corresponding 
results for the cosmic string geometry in the special case q = \ with zero planar angle deficit. 
Vacuum densities for both scalar and electromagnetic fields in the region between two coaxial 
cylindrical shells in background of Minkowski spacetime are investigated in Section [T71 

The use of general coordinate transformations in quantum field theory in flat spacetime leads 
to an infinite number of unitary inequivalent representations of the commutation relations with 
different vacuum states. In particular, the vacuum state for a uniformly accelerated observer, 
the Fulling-Rindler vacuum, turns out to be inequivalent to that for an inertial observer, the 
Minkowski vacuum. In Section [THl the polarization of the Fulling-Rindler vacuum for a scalar 
field by a uniformly accelerated plate with Robin boundary conditions is investigated by using 
the GAPF. The corresponding eigenfrequencies are the zeros of the function Ki^{x), where 
the coefficients in the barred notation (jS.ip are related to the Robin coefficients by formula 
(118. 4p . On the base of the summation formula derived in Section [5] for the series over these 
zeros, we have extracted from the Wightman function the part corresponding to the Rindler 
wedge without boundaries. The part induced by the plate is given by formula (|18.10p and is 
directly used for the evaluation of the VEVs for the field square and the EMT. In the case 
of D = 3 electromagnetic field with perfect conductor boundary condition on the plate, the 
corresponding eigenfunctions for the vector potential are resolved into the TE and TM scalar 
modes with Dirichlet and Neumann boundary conditions respectively. The corresponding VEV 
of the EMT is given by formula (jl8.26p . The Casimir densities in the region between two infinite 
parallel plates moving by uniform proper accelerations are discussed in Section [19] for both scalar 
and electromagnetic fields. In the case of scalar field the eigenfrequencies are the zeros uj = of 
the function Zi^[u,v) defined by formula (j5.14p (see (|19.3p ). After the application of summation 
formula (I5.18|) . the Wightman function is presented in the form ()19.7p or equivalently (119. 12p . 
The first term on the right of these formulae is the Wightman function for the case of a single 
plate and the second term is induced by the presence of the second plate. Similar decomposition 
is obtained for the VEVs of the field square and the EMT, formulae (jl9.13p and (|19.19p . We 
have investigated the parts in the VEVs induced by the presence of the second brane in various 
asymptotic limits of the parameters. The vacuum interaction forces between the plates are 
considered as well. The VEV of the EMT for the electromagnetic field in the region between the 
plates is given by formula (|19.30p or equivalently by ()19.3ip . and the corresponding interaction 
forces are given by formulae (|19.32p . (119. 33p . 

Recent proposals of large extra dimensions use the concept of brane as a sub-manifold em- 
bedded in a higher dimensional spacetime, on which the Standard Model particles are con- 
fined. Braneworlds naturally appear in string/M-theory context and provide a novel setting 
for discussing phenomenological and cosmological issues related to extra dimensions. In Sec- 
tion [20] we consider the geometry of two parallel flat branes in the AdS bulk. In the re- 
gion between the branes the radial Kaluza-Klein masses m = rUn are zeros of the function 
jl"'\mZa)Yy^\mzij) — Yu"'\mza)jl^\mzi,), where the coefficients in barred notation (|4.2p are 
related to the Robin coefficients by formulae (120. 8p . The application of summation formula 
()4.14p allows to extract from the VEVs the part due to a single brane. Both single brane 
and second brane induced parts and vacuum interaction forces between the branes are investi- 
gated. An application to the Randall- Sundrum braneworld with arbitrary mass terms on the 
branes is discussed. Similar issues for two spherical branes in Rindler-like spacetimes are con- 
sidered in Section [2T] The corresponding eigenfrequencies to = Vln are zeros of the function 
^i!^ i^i^) ~ i^i^)^ where a and b are the branes radii and A; is defined 

by (|21.6p . In Section [22] we apply the GAPF for the investigation of the radiation intensity 
from a charge moving along a helical trajectory inside a dielectric cylinder. Summation formula 
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(|22.16p is derived for the series over the zeros of function (|2'2.9|) . This formula is used for the 
summation of the series over the TE and TM eigenmodes of the dielectric cylinder appearing in 
the expressions for the corresponding radiation intensities on the lowest azimuthal harmonic. 

Of course, the applications of the summation formulae obtained from GAPF are not restricted 
by the Casimir effect only. Similar types of series will arise in considerations of various physical 
phenomenon near the boundaries with spherical and cylindrical symmetries, for example in 
calculations of the electron self-energy and the electron anomalous magnetic moment (for similar 
problems in the plane boundary case see, e.g., [146j and references therein). The dependence of 
these quantities on boundaries originates from the modification of the photon propagator due 
to the boundary conditions imposed by the walls of the cavity. 
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